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CHAPTER    I 
INTRODUCTION 

The  wide  diversity  and  essential  nature  of  the  many  products 
obtained  from  wheat  makes  it,  perhaps,  the  most  important  commod- 
ity produced  by  man.  The  many  uses  for  wheat  byproducts  Is  a 
direct  result  of  the  complexity  of  the  composition  and  structure 
of  each  individual  wheat  kernel.  There  Is  currently  a  large 
number  of  different  varieties  of  wheat  grown,  some  of  which  have 
only  recently  been  developed.  It  is  likely  that  more  such  vari- 
eties will  be  developed  in  the  near  future.  Among  these,  dis- 
tinct classes  have  been  defined  within  which  wheat  varieties  with 
similar  physical  and  compositional  properties  are  placed.  In 
addition  to  variety,  a  number  of  other  factors  such  as  soil  and 
environment  may  also  influence  the  composition,  quality,  and 
characteristics  of  wheat  (1).  Accurate  categorization  of  wheat, 
on  the  basis  of  its  physical,  or  structural  properties,  has 
proven  to  be  Important  in  commercial  production,  where  certain 
properties  are  often  found  to  be  necessary  or  desirable  for  a 
given  application.  Since  the  demand  for  a  specific  type  of  wheat 
by  the  commercial  sector  largely  dictates  Its  ultimate  monetary 
value,  the  category  within  which  a  given  variety  is  placed  is 
also  of  importance  to  the  farmers  who  supply  the  wheat  and  plant 
breeders  who  develop  new  wheat  varieties  for  production.  It  is 
therefore  Important  that  unambiguous,  universally-acceptable 
standards    be    applied    to    wheat    classification    and    that    a    suffic- 


iently  objective  method  be   employed   for  achieving  accurate  place- 
ment  of   wheat    samples    within    these    classification    standards. 

1.1      Current   Grading    Standards 

Although  distinctions  on  the  basis  of  the  botany  of  the 
wheat  plant  itself  are  the  most  fundamental  means  of  classifica- 
tion, such  a  classification  scheme  is  not  practical  in  commercial 
grain  trade.  Kernel  characteristics  such  as  hard  or  soft,  either 
red  or  white,  and  spring  and  winter  habit  are  commonly  used 
throughout  the  world  for  wheat  classification  (2),  although  spe- 
cific market  classes  using  such  distinctions  have  a  tendency  to 
vary  somewhat  among  different  countries.  In  the  United  States, 
commercial  wheat  is  divided  into  the  following  seven  major  cate- 
gories: hard  red  spring,  durum,  red  durum,  hard  red  winter,  soft 
red  winter,  white,  and  mixed  wheat  (red  durum,  though  still 
recognized,  has  essentially  disappeared  as  a  market  class)  (3). 
Defined  beneath  most  of  these  classes  are  several  subclasses, 
into  which  individual  varieties  can  be  placed.  Further  distinc- 
tions are  made  on  the  basis  of  grain  quality,  with  grades  of  1-5 
and  sample  grade  given  on  the  basis  of  percent  damaged  kernels, 
amount  of  foreign  material  present,  etc.  Figure  1.1  gives  a  gen- 
eral outline  of  the  standards  currently  used  by  the  United  States 
in  classifying  and  grading  wheat  (4).  Traditionally,  the  grain 
classification  and  grading  process  has  involved  few  scientific 
procedures,  and  many  subjective  judgements.  The  ability  to  make 
such    judgements    has,    however,    deteriorated    greatly    in    recent 


FIGURE  1.1.   United  States  wheat  grading  system. 


years  due  to  the  introduction  by  plant  breeders  of  new  wheat 
varieties,  crossed  between  varieties  belonging  in  different  mar- 
ket   classes. 

The  criteria  needed  to  effectively  determine  the  suitability 
of  a  given  sample  of  wheat  for  a  particular  end-use  includes  more 
information  than  is  provided  by  the  classification  scheme  de- 
scribed above.  Users  of  wheat  often  require  further  knowledge  of 
variable  factors  such  as  moisture  and  protein  content,  relative 
kernel  hardness,  gluten  quality,  etc.  Much  attention  is  pres- 
ently being  given  to  more  inclusive  consideration  of  such  quality 
factors  in  the  grading  process  (5).  Current  grading  techniques, 
which  are  primarily  based  on  visual  inspection,  are  clearly 
inadequate  for  making  these  types  of  quality  determinations.  To 
facilitate  more  accurate  and  comprehensive  testing  of  grain  sam- 
ples, analytical  methods  must  be  developed  which  provide  quanti- 
tative   information   about    the   properties   of   interest. 

1.2      General    Structure    and    Composition    of   Wheat   Kernels 

A  wheat  kernel  is  most  accurately  described  as  dry,  one- 
seeded  fruit,  known  to  the  botanist  as  a  caryopsis  (6).  Although 
the  grains  of  different  varieties  of  wheat  may  vary  considerably 
in  size,  form,  and  color,  they  all  possess  a  similar  fundamental 
structure.  Figures  1.2  and  1.3  provide  longitudinal  and  cross- 
sectional  views  of  a  typical  wheat  kernel  (7).  These  figures 
reveal  the  general  shape  of  the  kernel  and  give  an  indication  of 
its    structural    complexity.       The    kernel    structure    may    be    broadly 
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FIGURE    1.2.       Longitudinal    section    of    a    wheat    kernel. 
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FIGURE  1.3.   Cross-sectional  view  of  a  wheat  kernel, 


divided  into  three  primary  parts,  these  being  the  germ  and  endo- 
sperm, which  are  enclosed  in  a  seed  coat  (bran).  Minor  portions 
of  the  kernel  not  included  in  these  regions  are  a  tuft  of  fine 
hairs  located  at  one  end  of  the  kernel  and  a  deeply  pigmented 
vascular  strand  located  at  the  bottom  of  the  crease.  The  rela- 
tive location  of  each  of  the  primary  kernel  parts  is  identified 
in  the  cross-sectional  diagram.  Approximate  proportions  of  the 
kernel  comprised  by  these  three  parts  are:  endosperm,  83%;  germ, 
2.5%;  and  bran,  14.0%.  Reported  compositional  values  tend  to 
vary  somewhat,  however,  depending  on  methods  of  separation  and 
analysis  (5).  Each  of  these  primary  regions  of  the  kernel  has  a 
relatively  complex  individual  structure  of  its  own.  Figure  1.4 
provides  a  detailed  outline  of  the  major,  separately-identifiable 
portions  of  a  wheat  kernel,  and  their  relation  to  one  another 
(8).  The  longitudinal  diagram  illustrates  the  respective  loca- 
tion  of   most   of    these    structural    regions. 

The  chemical  constituents  of  wheat  are  fairly  large  in 
number  and  are  distributed  non-homogeneous ly  throughout  the  ker- 
nel. These  include  carbohydrates,  proteins,  lipids,  vitamins, 
water,  enzymes,  and  various  minerals.  The  relative  amounts  of 
these  constituents  in  a  given  kernel  can  vary  widely,  due  to 
differences  in  variety  and  environmental  conditions.  The  major 
components,  and  the  percentage  of  each  of  the  primary  regions  of 
a  typical  kernel  comprised  by  these  components,  are  outlined  in 
the    following    table    (9). 
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Table    1.1.      Composition    of    the    Endosperm,    Germ,    and    Bran. 


Constituent 


Endosperm 


Germ 


Bran 


Moisture 

14.0% 

Protein 

9.6% 

Fat 

l.« 

Mineral 

0.7% 

Starch 

71.0% 

Hemicellulose 

1.8% 

Sugars 

1.1% 

Cellulose 

0.2% 

Carbohydrate    (total)      74.1% 


11.7% 
28.5% 
10.4% 

4.5% 
14.0% 

6.8% 
16.2% 

7.5% 

44.5% 


13.2% 

14.4% 
4.7% 
6.3% 
8.6% 

26.2% 
4.6% 

21.4% 

60.8% 


The  endosperm,  which  comprises  most  of  the  wheat  kernel,  is 
of  primary  commercial  interest,  since  most  products  produced  from 
wheat  utilize  material  from  this  region.  The  endosperm  mainly 
consists  of  starch  granules  embedded  in  a  matrix  which  is  largely 
protein.  Starch  is  by  far  the  most  abundant  carbohydrate  present 
in  the  wheat  kernel,  constituting  approximately  64%  of  the  dry 
matter  of  the  entire  wheat  grain  and  about  70%  of  the  endosperm 
(10).  The  starch  contained  in  the  endosperm  is  in  the  form  of 
lenticular  or  spherical  granules  which  are  tightly  packed  toge- 
ther, with  the  protein  filling  the  inter-granular  spaces.  The 
structure  and  composition  throughout  the  endosperm  is  not  com- 
pletely uniform,  however,  with  the  concentration  of  protein  de- 
creasing as  the  central  portion  of  the  endosperm  is  approached. 
The  concentration  of  protein  is  known  to  be  dependent  upon  envi- 
ronmental   factors    as   well    as    variety,    with  hard   red   winter  wheats 


typically    containing    a    higher    protein    content    than    soft    red 
winter    varieties. 

Although  the  germ  and  bran  contain  much  lower  percentages  of 
starch  than  does  the  endosperm,  the  concentration  of  protein  in 
these  regions  is  greater.  The  concentration  of  most  other  con- 
stituents of  the  kernel  are  also  greater  in  these  regions.  The 
composition  of  the  multi-layered  bran  varies  to  a  large  extent 
with  increasing  depth  beneath  the  kernel  surface,  as  the  various 
structural  portions  of  this  region  are  encountered.  Likewise, 
the  components  of  the  germ  are  non-homogeneous ly  distributed, 
with  their  relative  concentrations  varying  with  its  structural 
complexity.  A  summary  of  more  extensive  and  detailed  information 
regarding  the  structure  and  composition  of  a  wheat  kernel  has 
been   compiled    by    MacMasters,    Hinton,    and    Bradbury    (11). 

1.3  Kernel  Hardness 

Hardness  is  one  of  many  distinguishing  characteristics  of 
wheat  which  may  vary  from  sample  to  sample.  As  properties  of 
wheat  kernels,  hardness  and  softness  relate  to  the  manner  with 
which  the  endosperm  of  the  kernel  breaks  down  during  the  milling 
process.  In  hard  wheats,  fragmentation  of  the  endosperm  tends  to 
occur  along  the  lines  of  the  of  the  starch  and  protein  cell 
boundaries,  whereas  in  soft  wheats  the  endosperm  fractures  in  a 
random  way.  This  seems  to  suggest  a  pattern  of  areas  of  mechani- 
cal strength  and  weakness  in  hard  wheats  and  uniform  mechanical 
weakness  in  soft  wheats.   Upon  milling,  the  degree  of  mechanical 
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damage  to  the  starch  granules  of  the  endosperm  is  greater  for 
hard  wheats  than  for  soft,  and  hard  wheat  typically  requires  a 
greater  amount  of  grinding  force  in  the  milling  process.  Hard- 
ness also  affects  the  ease  of  detachment  of  the  endosperm  from 
the  bran,  as  the  bran  is  found  to  separate  more  cleanly  with  hard 
wheats  than  with  soft.  Although  it  is  not  completely  understood 
why  some  wheats  are  "harder"  than  others,  one  view  is  that  hard- 
ness is  related  to  the  degree  of  adhesion  between  the  starch 
particles  and  protein  matrix  which  comprise  the  endosperm  (12). 
This  hypothesis  is  supported  in  research  done  by  Hoseney  and  Seib 
(13),  where  it  was  shown  that  sliced  hard  wheat  kernels  typically 
dissociate  by  rupturing  of  the  starch  particles,  whereas  in  soft 
wheats  the  dissociation  occurs  at  the  starch-protein  bonds,  leav- 
ing the  starch  particles  intact.  Another  view  is  that  kernel 
hardness  is  dependent  upon  the  continuity  of  the  protein  matrix 
of  the  endosperm  within  which  the  starch  particles  are  embedded. 
Hardness  is  a  characteristic  known  to  be  inherited  in  the  Mendel- 
ian  fashion  (14),  which  implies  that  hardness  should  not  greatly 
vary  for  kernels  of  a  given  variety. 

Due  to  the  above-cited  differences,  the  type  of  flour  ob- 
tained from  hard  wheats  is  different  than  that  obtained  from 
soft,  tending  to  be  more  coarse,  gritty,  and  more  freely  flowing 
and  easily  sifted.  Soft  wheats  yield  a  very  fine  flour  consist- 
ing of  irregular-shaped  fragments  of  endosperm  cells  which  adhere 
together,  sift  with  difficulty,  and  tend  to  clog  the  apertures  of 
sieves  (12).   Flour  from  hard  wheats,  due  to  its  granular  nature, 
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as  well  as  Its  relatively  high  protein  content  and  desirable 
gluten  quality,  Is  well-suited  for  breadmaklng.  Soft  wheat 
flour,  however,  Is  found  to  be  preferable  for  production  of  cake 
and  pastry  products  (15).  Since  hard  and  soft  wheats  require 
different  processing  techniques  for  optimum  utilization,  It  Is 
important  that  knowledge  of  the  relative  hardness  of  wheat  to  be 
processed  be  obtained  prior  to  its  use. 

In  the  classification  of  wheat,  hardness  Is  usually  judged 
on  the  basis  of  visual  appearance  rather  than  actual  measurement. 
For  research  purposes,  objective  measurements  of  hardness  are 
often  made  by  determining  the  wheat's  "pearling  index",  which  is 
defined  as  percentage  of  material  pearled-off  from  a  sample  of 
wheat  of  prescribed  weight  in  a  laboratory  barley  pearler  oper- 
ated for  a  prescribed  period  of  time.  The  equipment  and  proced- 
ures for  this  technique  have,  however,  never  been  standardized. 
Furthermore,  in  studying  segregates  produced  from  crossing  hard 
and  soft  wheats,  it  has  been  found  that  the  pearling  technique 
produces  wide  differences  In  the  pearling  indexes,  and  visual 
classification  for  such  material  is  frequently  faulty  (16).  In 
addition  to  this  procedure,  a  number  of  other  tests  have  been 
designed  to  determine  kernel  hardness  (17),  with  varying  degrees 
of  success.  At  present,  no  single  one  of  these  techniques  is 
either  universally  used  or  accepted.  These  tests,  in  general, 
require  samples  consisting  of  a  relatively  large  number  of  ker- 
nels In  order  to  yield  meaningful  results.  A  satisfactory  method 
for  determining  the  hardness  of  a  single  kernel  has  not  yet  been 
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developed. 

It  is  suggested  in  the  above  discussion  that  the  hardness  of 
a  wheat  kernel  is  related  to  structural  properties  of  the  endo- 
sperm. It  is  therefore  quite  possible  that  the  bulk  thermal 
and/or  optical  properties  of  hard  wheats  differ  to  some  extent 
from  soft  wheats.  If  the  variation  between  hard  and  soft  wheats 
is  sufficiently  large  and  is  consistent,  it  might  be  feasible  to 
make  a  determination  of  the  relative  hardness  of  a  single  kernel 
on  the  basis  of  these  properties.  Since  hardness  is  genetically 
transmitted,  it  should  also  be  possible  to  make  a  determination 
of  the  hardness  of  a  kernel  by  positively  identifying  it  as  being 
of  a  particular  variety,  assuming  that  the  relative  hardness  of 
that  variety  is  understood.  A  method  for  making  distinctions  on 
the  basis  of  the  optical  or  thermal  properties  of  a  kernel  might, 
therefore,  prove  successful  in  determining  relative  kernel  hard- 
ness, either  by  providing  information  reflecting  properties  which 
are  directly  related  to  hardness,  or  Indirectly  by  identifying 
properties  which  are  uniquely  common  to  a  given  variety. 

1.4   Photoacous tic  Analysis  of  Wheat  Kernels 

A  sample  analysis  technique  utilizing  a  phenomenon  commonly 
known  as  the  pho toacous t ic  effect  was  studied  and  subsequently 
used  in  this  research  project  for  the  examination  of  Individual 
wheat  kernels.  The  pho toacous tic  effect  occurs  when  a  substance 
in  a  closed  cell  is  illuminated  by  chopped  light,  or  some  other 
intensity-modulated  energy  source.   The  result  of  the  illumina- 
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tion  is  an  acoustic  signal  with  frequency  equal  to  the  modulation 
or  chopping  frequency  of  the  incident  beam.  The  effect  is  more 
thoroughly  discussed  in  the  following  chapter.  The  photoacous tic 
response  of  a  material  is  known  to  be  dependent  upon  that  mater- 
ial's optical  and  thermal  properties.  Thus,  it  should  be  possi- 
ble to  differentiate  between  wheat  kernels  in  which  optical  or 
thermal  properties  vary  by  comparing  the  photoacous tic  responses 
obtained  from  them.  If  thermal  or  optical  properties  vary  be- 
tween varieties,  or  these  properties  differ  for  hard  and  soft 
kernels,  the  variations  should  be  reflected  in  their  respective 
photoacous tic    responses,    which   should    vary   as   well. 

As  a  method  for  making  determinations  on  the  basis  of  such 
factors  as  kernel  hardness,  a  technique  incorporating  photoacous- 
tic  sample  analysis  would  offer  several  advantages  over  analysis 
techniques  currently  used.  First,  the  test  would  provide  a  non- 
destructive examination  of  a  sample,  enabling  use  of  the  sample 
for  further  testing  or  continued  research.  Secondly,  photoacous- 
tic  analysis  accommodates  small  sample  sizes,  and  would  thus 
permit  evaluation  on  the  basis  of  a  sampling  as  small  as  a  single 
kernel.  Although  small  sample  size  is  not  imperative  in  most 
instances  requiring  the  grading  of  wheat,  such  as  milling  opera- 
tions and  International  grain  trade,  it  might  be  advantageous  for 
certain  grain  research  and  plant  breeding  applications.  Also, 
photoacoustic  sample  analysis  is  capable  of  providing  quantita- 
tive comparisons,  and  would  thus  allow  for  objective  determina- 
tions  of   kernel    properties,    relative    to    other   kernels. 
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1.5   Summarization  of  Thesis  Objectives 

Although  the  primary  objective  of  this  project  was  to  exam- 
ine the  photoacoustic  response  of  wheat  kernels,  and  attempt  to 
correlate  it  to  kernel  hardness,  a  number  of  other  preliminary 
procedures  and  experiments  needed  to  be  performed  as  well.  A 
general  outline  of  the  major  objectives  of  the  project  is  given 
below: 

1.  Develop  a  microcomputer-operated  data  acquisition  system 
to  work  with  the  photoacoustic  sampling  apparatus. 

2.  Test  the  apparatus  and  acquisition  system  by  checking 
acquired  data  with  theoretical  predictions,  and  make  any  neces- 
sary equipment  modifications. 

3.  Examine  the  effect  of  moisture  on  the  photoacoustic 
response,  and  determine  a  satisfactory  method  for  standardization 
of  moisture  content  for  wheat  samples  to  be  studied. 

4.  Obtain  the  photoacoustic  response  throughout  a  modula- 
tion frequency  range  of  1-1000  Hz  for  single  kernel  samples  of 
several  different  hard  and  soft  wheat  varieties,  and  compare 
results . 

A  theoretical  discussion  of  the  photoacoustic  effect  in 
solids,  a  description  of  the  apparatus  and  acquisition  system 
used,  and  a  presentation  and  discussion  of  the  results  obtained 
are  provided  in  the  subsequent  chapters  of  this  thesis. 
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CHAPTER    II 
PHOTOACOUSTIC    EFFECT    IN   SOLIDS;    THEORY 

2.1      Introduction 

In  order  to  understand  clearly  the  significance  of  the 
information  provided  by  pho toacous t Ic  examination  of  various 
wheat  varieties,  it  is  first  necessary  to  understand  the  theore- 
tical basis  for  the  photoacoustic  effect.  Since  all  photoacous- 
tic  studies  done  in  this  thesis  work  have  been  on  solid  samples, 
in  this  discussion  we  will  concentrate  on  the  photoacoustic 
effect    in    solids. 

The  photoacoustic  effect,  in  general,  refers  to  the  produc- 
tion of  an  acoustic  signal  by  a  sample,  either  solid,  liquid,  or 
gaseous,  when  that  sample  is  illuminated  by  an  Intensity  modulat- 
ed energy  beam,  usually  in  the  form  of  electromagnetic  radiation. 
It  was  initially  observed  in  the  early  1880's  when  Alexander 
Graham  Bell  detected  an  audible  sound  through  a  listening  tube 
connected  to  a  chamber  containing  a  solid  sample  surrounded  by  a 
gaseous  atmosphere,  through  which  the  sample  could  be  illuminated 
with  a  periodically  interrupted  sunlight  beam  (18).  Subsequent 
experimentation  shortly  following  this  discovery  revealed  that  an 
acoustic  signal  could  also  be  detected  when  an  intensity-modu- 
lated sunlight  beam  was  allowed  to  impinge  upon  various  liquids 
and  light-absorbing  gases.  Motivated  by  these  findings,  Tyndall 
and  Roentgen  found  that  the  effect  could  also  be  created  using 
chopped    artificial    light    (19).      At    the    time    of    its    discovery,     the 
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processes  responsible  for  production  of  the  acoustic  signal  were 
not  well  understood,  with  a  number  of  qualitative  theories  being 
proposed.  Interest  in  the  effect  was,  however,  soon  diminished 
as  it  was  generally  considered  merely  a  laboratory  curiosity  of 
no  practical  importance.  Although  approximately  fifty  years 
later  interest  in  the  photoacous tic  effect  in  gases  was  renewed, 
becoming  a  well-established  method  for  gas  analysis,  the  photoa- 
coustic  effect  in  solids  and  liquids  remained  essentially  forgot- 
ten until  the  early  1970's,  when  its  potential  for  analysis  of 
the  bulk  physical  properties  of  non-gaseous  materials  began  to  be 
realized . 

The  experimental  apparatus  used  by  Bell,  Tyndall,  and 
Roentgen  in  the  Initial  discovery  and  early  investigation  of  the 
photoacous tic  effect,  possessed  the  crude  essentials  necessary  to 
observe  the  effect.  This  apparatus  consisted  roughly  of  an 
intensity-modulated  light  source,  an  airtight  enclosure  through 
which  the  light  source  could  penetrate  containing  either  a  gase- 
ous sample  or  a  non-gaseous  sample  surrounded  by  a  gaseous  atmos- 
phere, and  an  acoustic  signal  detection  device  by  which  the 
photoacous tic  signal  could  be  monitored.  Figure  2.1  gives  a 
schematic  representation  of  such  a  photoacous tic  effect  detection 
arrangement  as  used  for  monitoring  the  photoacous tic  signal  emit- 
ted by  a  solid  sample.  Although  this  apparatus  alone  is  insuf- 
ficient for  performing  the  quantitative  studies  of  solid  samples 
as  done  in  this  thesis  work,  such  an  illustration  is  sufficient 
for  demonstrating  the  effect,  and  is  useful  in  describing  the 
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FIGURE    2.1.       Simple    arrangement    for    pho toacous t lc 
signal    detection. 


18 


basic  mechanisms  by  which  the  pho toacous tic  effect  arises. 

2.2   Fundamental  Description 

In  attempting  to  describe  the  underlying  physical  processes 
which  give  rise  to  the  pho toacous tic  effect  in  solids,  it  is 
perhaps  best  to  begin  by  giving  a  qualitative  discussion  of  these 
processes,  followed  by  a  more  rigorous  mathematical  treatment  of 
the  phenomenon. 

As  earlier  stated,  the  photoacoustic  effect  occurs  when  an 
intensity  modulated  energy  beam  is  allowed  to  impinge  upon  a 
sample  when  that  sample  is  confined  to  an  enclosed  chamber, 
generally  referred  to  as  the  photoacoustic  cell.  The  modulated 
energy  source  used  is  usually,  but  not  always,  monochromatic 
electromagnetic  radiation  in  or  near  the  visible  spectrum.  Un- 
less specified  otherwise,  incident  monochromatic  electromagnetic 
radiation  will  be  assumed  throughout  this  discussion.  The  photo- 
acoustic signal,  or  response,  is  in  the  form  of  an  acoustic 
signal  at  the  same  frequency  as  the  energy  beam  modulation  fre- 
quency. The  photoacoustic  signal  is  monitored  by  an  acoustic 
signal  detection  device  sufficiently  sensitive  to  measure  It, 
commonly  a  high  sensitivity  microphone  built  into  some  portion  of 
the  chamber.  In  the  case  of  photoacoustic  experimentation  with 
solids,  it  is  necessary  that  the  solid  sample  be  at  least  par- 
tially surrounded  by  a  gaseous  atmosphere,  preferably  using  a  gas 
that  will  absorb  a  minimal  amount  of  the  incident  energy. 

When  a  given  cycle  of  the  modulated  energy  beam  (e-m  radla- 
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tion)  becomes  incident  upon  the  sample,  the  sample  will  absorb 
all  or  part  of  this  energy  throughout  a  depth  beneath  the  illumi- 
nated surface  characterized  by  the  optical  absorption  depth  of 
the  sample  for  the  wavelength  of  the  e-m  radiation  used.  This 
energy  absorption  results  in  localized  excitation  of  atomic  ener- 
gy levels  within  the  sample  medium  at  the  absorption  depth.  When 
the  subsequent  deexcitation  of  these  energy  levels  occurs,  it 
generally  does  so  by  non-radiative  heat  producing  processes, 
which  results  in  a  localized  heating  of  the  sample  in  the  vicini- 
ty of  the  initial  absorption.  This  heat  then  diffuses  outward 
from  the  point  of  absorption  in  all  directions  as  Individual 
atoms  and/or  molecules  are  vibrational ly  excited  in  a  non-cooper- 
ative and  highly  dissipative  manner.  This  heat  diffusion  within 
the  sample  has  the  characteristics  of  a  highly  damped  wave,  being 
formally  equivalent  to  wave  propagation  in  a  very  lossy  medium, 
with  a  wavelength  approximately  equal  to  the  thermal  diffusion 
length  of  the  sample  in  the  region  through  which  it  propagates. 
Thus,  this  heat  propagation  through  the  sample  is  often  referred 
to   as    a    thermal    wave. 

If  the  initial  energy  absorption  occurs  at  a  shallow-enough 
depth  beneath  the  surface  so  that  the  thermal  wave  is  able  to 
propagate  to  the  surface  of  the  sample  before  being  entirely 
dissipated,  the  heat  which  reaches  the  surface  will  produce  an 
increase  in  temperature  at  the  sample  surface.  If  the  incident 
e-m  radiation  is  maintained  at  a  constant  modulation  frequency, 
this    will    result    in    a    periodic    heating    of    the    surface    of    the 

2  0 


sample,  with  the  heat  contribution  from  a  given  depth  beneath  the 
surface  of  the  sample  being  essentially  the  same  for  each  modula- 
tion cycle.  Because  a  finite  amount  of  time  is  required  for  the 
heat  produced  within  the  sample  to  propagate  to  its  outer  sur- 
face, a  thermal  wave  propagating  from  some  depth  beneath  the 
surface  of  the  sample  will  cause  a  heat  contribution  at  the 
sample  surface  which  lags,  to  some  extent,  in  phase  with  respect 
to  the  incident  e-m  source.  Since  a  gaseous  atmosphere  surrounds 
the  surface  of  the  sample,  this  heat  is  then  transferred  to  a 
thin  layer  of  gas  at  the  gas-sample  interface.  The  periodic 
heating  of  this  thin  gas  layer  causes  it  to  expand  and  contract 
at  a  frequency  equal  to  the  modulation  frequency  of  the  incident 
e-m  radiation,  creating  periodic  pressure  fluctuations  which 
occur  in  the  photoacoustic  cell  at  the  modulation  frequency,  and 
resulting  in  the  acoustic  signal  which  can  be  detected  by  a 
sufficiently   sensitive    microphone. 

From  the  above  discussion,  it  is  easily  seen  that  the  pho- 
toacoustic response  obtained  from  a  solid  sample  depends  on  the 
thermal  properties  of  the  sample,  which  influence  the  propagation 
of  the  thermal  wave  through  the  sample  to  the  sample  surface. 
Also,  the  response  is  dependent  upon  the  optical  properties  of 
the  sample,  which  dictate  the  distance  beneath  the  surface  at 
which  the  thermal  waves  will  be  generated.  Although  this  adds  to 
the  complexity  of  a  quantitative  treatment  of  the  problem,  both 
of  these  dependencies  are  individually  important,  for  each  can  be 
used    to   obtain   information   regarding    the   bulk   physical    properties 
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of  a  given  sample.  For  example,  If  the  modulation  frequency  of 
the  e-m  radiation  Is  held  constant  and  the  wavelength  of  the 
Incident  radiation  Is  allowed  to  vary  through  some  range,  Infor- 
mation similar  to  that  achieved  by  conventional  optical  absorp- 
tion spectroscopy  can  be  obtained;  as  strong  absorption  by  the 
sample  will  lead  to  thermal  waves  emerging  at  the  surface  before 
being  greatly  dispersed,  and  radiation  not  absorbed  will  not 
contribute  to  the  acoustic  signal.  Furthermore,  this  technique 
can  often  be  applied  In  Instances  where  the  opaqueness  of  the 
sample  or  problems  with  scattered  light  pose  serious  difficulties 
with  conventional  absorption  or  reflection  studies  (20).  Another 
important  application  of  the  pho toacous tic  effect  in  solid  sample 
analysis  takes  advantage  of  a  functional  dependence  of  the  ther- 
mal diffusion  length,  or  the  wavelength  of  the  thermal  waves 
which  contribute  to  the  measured  acoustic  signal  from  deepest 
within  the  sample,  on  the  modulation  frequency  of  the  incident 
electromagnetic  radiation.  This  technique,  commonly  referred  to 
as  pho toacous tic  depth  profiling,  is  unique  to  pho toacous tics , 
and  essentially  permits  one  to  obtain  a  depth  profile  of  the 
thermal  characteristics  of  a  sample  while  the  sample  itself 
remains  intact.  Both  the  phase  lag  of  the  acoustic  response  with 
respect  to  the  Incident  radiation  and  the  amplitude  of  the  signal 
provide  information  regarding  the  thermal  properties  encountered 
in  the  portion  of  the  solid  through  which  the  thermal  wave  propa- 
gates prior  to  reaching  the  surface.  Variation  of  the  thermal 
diffusion    length  as    the   modulation    frequency    is    changed    allows 

22 


information  from  different  sample  depths  to  be  acquired.  It  is 
this  aspect  of  photoacoustic  sample  analysis  that  was  primarily 
utilized    in    the    studies    done    throughout    this    thesis    work. 

In  the  above  discussion,  only  a  qualitative  description  of 
the  photoacoustic  effect  in  solids  has  been  given.  Although 
dependencies  of  the  effect  on  the  optical  and  thermal  properties 
of  the  sample  were  alluded  to,  the  exact  nature  of  these  depend- 
encies was  not  stated.  If  meaningful  quantitative  sample  analy- 
sis is  to  be  performed  utilizing  this  effect,  a  more  exact  theo- 
retical model  needs  to  be  employed.  Rosencwaig  and  Gersho,  in 
1976,  published  such  a  theoretical  model,  presently  referred  to 
as  the  Rosencwaig-Gersho,  or  R-G  theory  (19).  Although  this 
theoretical  description  Is  one-dimensional,  and  assumes  a  homo- 
geneous solid  sample,  it  yields  valuable  information  concerning 
many  aspects  of  the  photoacoustic  effect  in  solids.  Hence,  the 
main  points  of  this  theory,  and  some  of  its  predictions  in  spe- 
cial   cases   of   interest,    are   given  below. 

2.3      Rosencwaig-Gersho    Theory 

A.    Heat   Flow   Equations 

To  begin  this  theoretical  treatment  of  the  problem,  some 
basic  assumptions  first  need  to  be  made.  Let  us  assume  a  simple, 
closed,  cylindrical  photoacoustic  cell  with  length  L  and  diameter 
D,  as  illustrated  in  Figure  2.2.  This  cell  contains  a  disk- 
shaped  sample  of  length  l  which  has  the  same  diameter  D  as  the 
cell.      The    sample    is    positioned    in    the    cell    so    that    its    back 
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surface  Is  placed  directly  against  a  backing  material  of  low 
thermal  conductivity,  and  Its  front  surface  faces  the  gas  filled 
volume  of  the  cell.  Thegas  and  the  backing  materials  are  as- 
sumed non-absorbing  for  the  wavelength  X  of  the  incident  e-m 
radiation,  modulated  at  frequency  in  (rad/sec),  impinging  upon  the 
sample's  front  surface.  The  sample  itself  is  assumed  to  have  an 
absorption   coefficient    B   (1/cm). 

We   now   define    the    following    parameters    for    the   gas,    sample, 
and   backing   material;    where    the   subscript   i    takes    on    the    value    g 
for    the    gas,    s    for    the    solid,    and    b    for    the   backing   material: 
k»:    thermal   conductivity  (cal/cm-sec-deg)   of   material    i 
P^:    density   (g/cm   )   of   material    i 
Cji    specific    heat    (cal/g-deg)    of   material    i 
a1=k1/P1Ci:    thermal    dlffusivity    (cm    /sec)    of   material    i 
aj  =  (  i)/2c(1)1  '2:     thermal    diffusion    coef.    (1/cm)    of    material    i 
u^=l/aj:    thermal   diffusion   length  of  material    i 
Let   the   incident  e-m   radiation   be   sinusoidally   modulated 
with    Intensity    I,    given   by: 

I=-il0(l+cosiot), 

where  IQ  is  the  incident  monochromatic  light  flux  (H/cm2). 

Using  these  parameters,  the  heat  produced  at  a  given  point  x 
within  the  sample  due  to  the  light  absorbed  at  that  point  is 
given  by: 


il0exp(Bx)(l+coswt), 
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with  x  taking  on  negative  values  from  x=0,  where  the  light  is 
incident  upon  the  sample,  to  x=-l,  at  the  back  end  of  the  sample. 
The  thermal  diffusion  equation  for  the  interior  of  the  solid 
(-Kx<0),  taking  into  account  the  distributed  heat  source,  can  be 
written    as: 

fs^-r-H"  A  exp(Bx)[l+exp(jcot)] 

for   -Kx<0,  (1) 

with 

A-gl^i  /2k,, 

where  $  is  the  temperature  and  n  is  the  efficiency  at  which  the 
absorbed  light  is  converted  to  heat  by  the  non-radiative,  heat 
producing  deexcitation  processes.  Here,  we  will  assume  that  n=l, 
which  is  reasonable  for  most  solids  at  room  temperature.  For  the 
backing  material  and  the  gas,  the  respective  heat  diffusion 
equations  are  simplified,  since  the  optical  absorption  coeffic- 
ient   for    these    is    zero.      They   are: 

32(j>  1    3ifr 

3JF"       a.     3t  -9.-lh<x<-l  (2) 

q 
The    real    part   of    the   complex-valued    solutions    for    the    temp- 
erature   4>    for    the    above    thermal    diffusion    equations    is    the    solu- 
tion of   physical    interest,    representing   the    temperature   in   the 
cell    relative    to    the    ambient    temperature    as    a    function    of    posi- 

2(> 


tion  and  time.  Therefore,  the  actual  temperature  field  in  the 
cell  is  given  by: 

T(x,t)=Re[*(x, t)]+* 

with  "Re"  denoting  the  real  part  of  the  solution,  and  *  being  the 
ambient   (room)    temperature. 

In  order  to  completely  specify  the  solution  to  the  above 
heat  diffusion  equations,  the  appropriate  boundary  conditions 
need  to  be  applied.  These  boundary  conditions  are  obtained  from 
the  requirement  of  temperature  and  heat-flux  continuity  at  the 
boundaries  of  the  sample  (x=0  and  x=-2),  and  from  assuming  that 
the  temperature  at  the  cell  walls  x=  L  and  x—l-l.  is  ambient. 
Finally,  it  is  assumed  that  the  dimensions  of  the  cell  are  small 
enough  that  convective  heat  flow  in  the  gas  can  be  ignored  at 
steady    state    conditions. 

B.    Temperature    Distribution    in    the    Cell 

The  general  solution  for  cj>  (x,  t)  In  the  cell,  neglecting 
transients,     can   be   written   as: 

<f>(x,t)  =  (l/l>b)(x+S.+J>b)W0+Wexp[ab(x+n)+j(Dt], 

-i-*b<x<-« 

=e1+e2X+dexp(Bx)  +  [Uexp(  asx)+Vexp(-asx) 

-Eexp(Bx)]exp(jut),        -£<x<0 

=  (l-x/ Jlg)e0+9exp(-a  x+jut), 

0<x<£ 

8  (4) 
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where  U,  U,   V,   E,  and    8   are  complex-valued  constants;  e,,   e->,   d, 
U0,    and  80   are   real    valued   constants,    and   aji'd+jjaj,    with 
a1  =  (M/2o1)1/2. 

It  should  be  noted  that  6  and  W  represent  the  complex  ampli- 
tudes of  the  periodic  temperature  at  the  sample-gas  Interface 
(x=0)  and  the  sample-backing  material  boundary  (x=-£),  respec- 
tively. The  quantities  8Q  and  WQ  represent  the  real,  dc  compo- 
nent of  the  temperature,  relative  to  ambient,  at  these  boundaries 
(x=0  and  x=-i),  respectively.  The  dc  solution  in  the  backing 
material  and  gas  already  make  use  of  the  assumption  that  the 
temperature,  relative  to  ambient,  is  zero  at  the  ends  of  the 
cell.  The  quantities  E  and  d,  which  are  determined  by  the  forc- 
ing   function    in    Equation    (1),    are    given    by: 

d=-A/B2  (5a) 

and 

F  A  Big 

"   (B'-o*)   "  ?ks(B2-aa)  (5b) 

In  the  general  solution,  Equation  (4),  the  growing  exponen- 
tial component  of  the  solutions  to  the  gas  and  the  backing  mater- 
ial has  been  omitted,  since  for  all  frequencies  u  of  interest, 
the  thermal  diffusion  length  Is  small  compared  to  the  length  of 
both  the  gas  and  backing  materials.  Thus  the  sinusoidal  compo- 
nents of  these  solutions  are  sufficiently  damped  so  as  to  be 
essentially  zero  at  the  cell  walls,  hence  the  growing  exponential 
components    of    the    solutions    must   have   coefficients    essentially 


28 


equal  to  zero  In  order  to  satisfy  the  temperature  constraint  at 
the  cell  wal  Is. 

The  temperature  and  flux  continuity  conditions  at  the  sam- 
ple's front  and  back  surfaces  are  explicitly  given  by: 


and 


*g(0,t)-+s(0,t), 

*b(-e,t)  =  4,s(-i,t), 
kg  f^o.t)=ksfii(o,t) 

kbJ^-t.t).^(-i.t) 


(6a) 
(6b) 
(6c) 

(6d) 


where  the  subscripts  s,  b,  and  g  Identify  the  solution  to  equa- 
tion (4)  for  the  solid,  backing,  and  gas.  The  above-stated 
constraints  apply  separately  to  the  dc  component  and  the  sinusoi- 
dal component  of  the  solution.  From  Equation  (6),  we  obtain  for 
the    dc    component   of    the    solution: 


e0=e1+d, 

U0=e1-e2il+d'exp(-S2 ) , 


-(WVkse2+ks8d 


(7a) 
(7b) 

(7c) 
(7d) 


(kb/  lib)U=kse2+ksB  dexp(-6«  ) . 

Equations    (7)    determine    the   coefficients    e,,    e,,   W    ,    and    6      for 
the   time-independent  (dc)   component  of    the    solution.      Applying 
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Equations  (6)  to  the  sinusoidal  component  of  the  solution  yields: 
e=u+v-E,  (8a) 

W=Uexp(-as£ )+Vexp( asH)-Eexp(-BH ) ,  (8b) 

-kgage=ksosa-ksasV-ksBE,  (8c) 

and 

kbabW=Uexp(-osJ.)-ksosVexp(  asl  )-ks  BEexp(-gji  ) .  (8d) 

These  equations  together  with  the  expression  for  E  in  Equation 
(5b)  determine  the  coefficients  0,  V,  W,  and  6.  Hence  the  solu- 
tions to  Equations  (7)  and  (8)  allow  for  evaluation  of  the  temp- 
erature distribution,  Equation  (4),  in  the  cell  in  terras  of  the 
optical,  thermal,  and  geometric  parameters  of  the  system.  The 
explicit  solution  of  the  complex  amplitude  of  the  periodic  tem- 
perature  at    the   solid-gas   boundary   (x=0)    is    given   by: 

n-.      Ha     ..,(r-lHb+l)exo(gs£)-(r+l)(b-l)exp(-q,.z)+2{b-r)exp(-e)L) 
8~2ks(sy-a2)  (q+l)(b+l)exp(as*)-(g-l)(b-I)exp(-os&) 

(9) 
where 

b=Ml 

Ms  do) 

ksas  (11) 


r=(l-j)_B_> 

2as  (12) 

and 

cr.-d-JK 
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Thus,  Equation  (9)  can  be  evaluated  for  specific  parameter  values 
yielding  a  complex  number  whose  real  and  imaginary  parts  9j  and 
62  respectively,  give  the  ln-phase  and  quadrature  components  of 
the  periodic  temperature  variation  at  the  sample's  front  surface. 
Specifically,  the  actual  temperature  at  the  sample  surface  (x»0) 
is  given  by: 

T(0,  t)=«>+9o+e1coswt-82sin(i)t, 

where  $  is  the  ambient  temperature  (room  temperature)  at  the  cell 
walls  and  8Q  is  the  increase  in  temperature  from  the  steady  state 
component    of    the    absorbed    heat. 

C.    Production    of    the    Acoustic    Signal 

As  previously  stated,  the  acoustic  signal  detected  in  the 
photoacous tic  cell  arises  from  the  periodic  heat  flow  from  the 
sample  to  a  layer  of  the  gas  at  the  gas-sample  interface.  This 
periodic  heat  diffusion  causes  a  periodic  temperature  variation 
in  the  gas  given  by  the  sinusoidal  (ac)  component  of  Equation 
(4),    which   is: 


Sac(x, t)=9exp(-o  x+ju  t) . 


(13) 


Taking    the   real    part   of   Equation   (13),    it   can  be   seen   that 
the   actual    temperature    variation    in    the   gas    Is: 


Tac(x,t)=exp(-agx)[9  jcos  (u)  t-a   x)-e2sin(u  t-a_x)  ]  , 


(14) 


where  8j  and  82  are  the  real  and  imaginary  parts  of  6,  as  given 
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by  Equation  (9).  The  time-dependent  component  of  the  temperature 
in  the  gas  attenuates  rapidly  to  zero  with  increasing  distance 
from  the  surface  of  the  solid  sample,  and,  at  a  distance  of  only 
2irug  from  the  surface,  is  effectively  fully  damped  out.  He  can 
therefore  define  a  boundary  layer  of  thickness  2iru  ,  and  make  the 
approximation  that  only  this  gas  layer  is  capable  of  responding 
thermally  to  the  periodic  temperature  at  the  sample  surface. 

The  spatially  averaged  temperature  within  this  gas  boundary 
layer  can  be  determined,  as  a  function  of  time,  by  evaluating: 

From  Equation  (13)  this  gives: 


(15) 

using  the  approximation  that  exp(-2T)«l. 

Since  the  gas  boundary  layer  expands  and  contracts  periodi- 
cally due  to  the  periodic  diffusion  of  heat  into  this  layer  from 
the  sample  surface,  it  can  be  treated  as  an  acoustic  piston 
acting  on  the  rest  of  the  gas  column,  producing  an  acoustic 
pressure  signal  that  is  transmitted  through  the  entire  column. 
The  displacement  of  the  gas  piston  due  to  its  periodic  heating 
can  be  estimated  by  using  the  ideal  gas  law, 


(5x( 


t)  =  2,Uq41]-=>rexp[j(.t-V)]  (16) 

0     c  ■  o 
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where  we  have  set  the  average  dc  temperature  of  this  gas  boundary 
layer  equal  to  the  dc  temperature  at  the  sample  surface  (T  =»  + 
80).  Equation  (16)  Is  a  reasonable  approximation  to  the  actual 
displacement  of  the  layer  since  2ttu  is  only  approximately  0.1  cm 
for   <u/2tt  =100    Hz    and    even    smaller    for  higher   frequencies. 

If  we  assume  that  the  rest  of  the  gas  responds  to  the  action 
of  this  piston  adiabatical ly,  then  the  acoustic  pressure  in  the 
cell  due  to  the  displacement  of  the  gas  piston  is  derived  from 
the   adiabatic   gas    law: 

PVY=constant, 

where  P  is  the  pressure,  V  is  the  gas  volume,  and  y  is  the  ratio 
of    the    specific    heats.       Thus    the    incremental    pressure    is: 

6P(t)    =-^6V  =-*%x(t) 

0  *g 

where  PQ  and  VQ  are  the  ambient  pressure  and  volume,  and  -617  is 
the   incremental    volume.      Then   from   Equation   (16)   we   have: 

SP(t)-Qexp[j(wt-%ir)],  (17) 

where 


Q=i^ 


VgTo  (18) 


Therefore,  the  actual  physical  pressure  variation,  Ap(t),  is 
given  by  the  real  part  of  <5p(t),  and  is: 

Ap(t)»Q1cos(a)t-%ir)-Q2sin(u>t-%ir),  (19) 
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AP(  t)=qcos(n)  t-4  -H*)  ,  (20) 

where  Qj  and  Q2  are,  respectively,  the  real  and  Imaginary  parts 
of  Q,  and  q  and  -i|i  are  the  magnitude  and  phase  of  Q,  thus, 

Q=Ql+jQ2=?exP(-:N'  )• 

The  value  Q  specifies  the  complex  envelope  of  the  sinusoidal 
pressure  variation.  Combining  Equations  (9)  and  (18)  we  get  the 
explicit  formula: 

n=  B'n  Tpn (r-l)(b+l)exp(ga)-(r+l)(b-l)exp(-o.t)+2(b-r)exp(-BZ) 

^VqagT°(s  -°s' (9+l)(b+iyexp(asl)-(g-l)(b-l)exp(-asn) 

(21) 

where  b  =  kbab/ksas,  g  =  kgag/ksas,  r  =  (l- j  )B/2as ,  and  as  =  (l  +  j)as  as 
previously  defined.  At  ordinary  temperatures  TQ  Is  essentially 
equal  to  *  ,  thus  the  dc  components  of  the  temperature  distribu- 
tion need  not  be  evaluated.  Hence,  Equation  (21)  may  be  evalua- 
ted for  the  magnitude  and  phase  of  the  pho toacous tic  response, 
represented  by  the  acoustic  pressure  wave  produced  In  the  cell. 

D.  Special  Cases 

Although  the  full  expression  for  <5P(t)  is  somewhat  difficult 
to  interpret  due  to  the  complicated  expression  for  Q  as  given  by 
Equation  (21),  physical  insight  may  be  gained  by  examining  spec- 
ial cases  of  Interest,  where  the  expression  for  Q  can  be  simpli- 
fied. Some  such  special  cases  for  optically  opaque  and  optically 
transparent  samples  are  considered  here,  with  opaqueness  being 
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determined  by  the  relation  of  the  optical  absorption  length, 

vi/e 

to  the  thickness  I  of  the  solid.  For  each  category  of  optical 
opaqueness,  three  cases  are  considered  according  to  the  relative 
magnitude  of  the  thermal  diffusion  length  ug,  as  compared  to  the 
physical  length  I  and  the  optical  absorption  length  yg.  For  all 
of  these  cases,  we  make  the  reasonable  assumptions  that  g<b  and  b 
is    approximately    one.      Also,    it    is    convenient    to    define 


rajo 

>/F8,gT0  (22) 


which  appears  in  the  expression  for  Q  in  each  case.  The  six 
cases  are  illustrated  in  Figure  2.3  (19)  and  are  individually 
considered    in    the    following    discussion. 

Case    1:    Optically   Transparent    Solids    (yg><0 

In  these  cases,  absorption  of  light  occurs  throughout  the 
entire    length   of    the   sample. 

Case    1(a):      Thermally    thin    solids    (us>>H;   Ps>ug) 

Here  we  set  exp(-6fc)  -  1-W,  exp(±o  <.(>)  =  1,  and  |r|>l  in 
Equation    (21).       We    then    obtain: 


■"gVb  s  Zag  D  (23) 

The  acoustic  signal  is  thus  proportional  to  S£  and,  since  vb/ag 

is  proportional  to  1/u,  the  acoustic  signal  has  a  u"   dependence. 
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Since  the  thermal  diffusion  length  of  the  sample  is  greater  then 
the  sample  length,  the  thermal  properties  of  the  backing  material 
affect  the  photoacous tic  response,  and  the  expression  for  Q  is 
thus  a  function  of  these  properties. 

Case  1(b):  Thermally  thin  solids  (u„>8,;  p„<u„) 

S  S     g 

Here  we   set   exp(-g{,)   -    1-61,    exp(±os«)    =    (l±as«),    and    |r|<l 
in    Equation    (21).       We    then    obtain: 

5  -  4kYa  a   3br(B2+2as2H(g2.2a$2)] 
s  g  s 

(24) 

Here,  the  acoustic  signal  is  again  proportional  to  BH,  varies  as 
w  ,  and  depends  on  the  thermal  properties  of  the  backing  mater- 
ial as  well.  The  expression  for  the  acoustic  signal  in  this  case 
is  identical  to  Equation  (23). 

Case  1(c):  Thermally  thick  solids  (y„<£  :  u„«u„) 

&      s    g 

Here,  the  approximations  that  are  made  are  exp(-£S£)  =  l-gst, 
exp(-as«)  =  0,  and  |r|<<l.  The  acoustic  signal  then  becomes: 

9  (25) 

The  acoustic  signal  here  is  now  proportional  to  By g  rather  than 
St.  That  is,  only  the  light  absorbed  within  the  first  thermal 
diffusion  length  of  the  sample  contributes  to  the  signal.  Also, 
since  Ps<  ?.,  the  thermal  properties  of  the  backing  material  pres- 
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ent  in  Equation  (24)  are  replaced  by  those  of  the  solid.  In  this 
expression,  the  dependence  of  Q  on  the  modulation  frequency 
varies    as   u)~3/2_ 

Case    2:    Optically   Opaque    Solids    (\i0«l) 

p 

In  these  cases,  most  of  the  light  being  absorbed  by  the 
sample  is  within  a  distance  small  compared  to  the  length  I  of  the 
sample. 

Case  2(a):  Thermally  thin  solids  (us>>£;  M3»y.) 
In  Equation  (21),  we  set  exp(-6«  )  -  0,  exp(±o^  )   =  1,  and 
lr|>>l.   We  then  obtain: 


9  (26) 


In  this  case,  we  have  pho toacous tic  "opaqueness"  as  well  as 
optical  opaqueness,  in  that  the  acoustic  signal  is  independent  of 
8,  This  would  be  the  case  of  a  very  black  absorber  such  as  carbon 
black.  The  acoustic  signal  Is  quite  strong,  and  depends  on  the 
thermal  properties  of  the  backing  material.  For  this  case,  the 
dependence  of  the  signal  on  frequency  Is  given  as  a  function  of 
oi-l. 

Case   2(b):      Thermally   thick   solids    (us<£;    ^S>P8) 
In    Equation    (21),    we    set    exp(-B8.)    =    0,    exp(-osH)    =    0,     and 
|rl>l.       We    obtain: 


9  (27) 
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Equation  (27)  is  analogous  to  Equation  (26),  but  the  thermal 
parameters  of  the  backing  are  now  replaced  by  those  of  the  solid. 
Again,  the  acoustic  signal  is  independent  of  the  optical  absorp- 
tion coefficient   and    varies    as  a "*, 

Case   2(c):      Thermally   thick   solids    (v  «{,;  u    <p    ) 
Here,    again,    we    set    exp(-BJt)    =    0,     exp(-as«)    =    0,     and    let 
|r|<l    in    Equation    (21).       We    obtain: 

(28) 

In  this  case,  although  the  solid  is  optically  opaque,  it  is  not 
photoacoustically  opaque,  since  only  the  light  absorbed  within 
the  first  thermal  diffusion  length  will  contribute  to  the  acous- 
tic signal.  Even  though  the  solid  is  opaque,  the  photoacous tic 
signal  will  be  proportional  to  the  optical  absorption  coeffic- 
ient.  As  in  case  1(c),  the  signal  varies  as  id"3/2. 

2.4  Summary 

The  Rosencwaig-Gersho  theory  presented  above  provides  a 
derivation  of  the  exact  solutions  for  the  acoustic  pressure 
produced  in  an  enclosed  cell  by  the  photoacous  tic  effect.  These 
equations  are  then  evaluated  to  give  specific  solutions  for 
special  cases  of  physical  significance.  Three  characteristic 
lengths  are  important  in  determining  which  of  the  specific  solu- 
tions are  applicable  for  a  given  sample,  these  being  the  thermal 
diffusion   length,   optical   absorption   length,   and    the   physical 
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thickness  of  the  sample.  Since  our  experiments  primarily  in- 
volved wheat  kernels,  which  are  optically  opaque,  only  the  theo- 
retical cases  for  which  the  optical  absorption  length  of  the 
sample  is  less  than  the  length  of  the  sample  itself  are  relevant 
with  respect  to  our  results.  It  should  be  noted  that  the  theory 
assumes  one-dimensional  pho toacous tic  apparatus,  no  reflection 
from  the  sample  surface,  100%  non-radiative  deexcltation  within 
the  sample,  a  non-absorbing  gas,  and  a  gas  thermal  diffusion 
length  less  than  the  length  of  the  column  of  gas  contained  in  the 
cell.  These  assumptions  will  generally  not  all  be  completely 
valid  in  practice,  and  may  affect  the  agreement  of  experimental 
results  with  theoretical  predictions.  The  following  chapter 
gives  a  detailed  description  and  discussion  of  the  actual  photoa- 
coustic  apparatus  which  we  used  in  the  experiments  performed 
throughout    this    thesis    work. 
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CHAPTER    III 
EXPERIMENTAL    APPARATUS    FOR    PHOTOACOUSTIC   ANALYSIS 

All  sample  analysis  studies  utilizing  the  pho toacous tic 
effect  are  similar  in  that  they  require  the  basic  components 
needed  to  produce  the  effect;  these  being  an  enclosed  sample 
cell,  some  intensity  modulated  energy  source  made  incident  upon 
the  enclosed  sample,  a  suitable  acoustic  signal  detection  device, 
and  a  method  by  which  this  signal  can  be  recorded  and  analyzed. 
Even  so,  the  particular  components  used  and  the  methods  by  which 
they  are  employed  will  vary  widely  with  respect  to  the  material 
being  studied  and  the  information  sought  regarding  this  material. 
In  this  section,  the  actual  experimental  apparatus  and  analysis 
techniques    used    in    these    studies    will    be    described    and    discussed. 

Shown  In  Figure  3.1  is  a  schematic  diagram  of  the  photoa- 
coustic  sampling  apparatus  used  in  these  studies.  The  primary 
components  of  this  setup  include:  1,  a  photoacoustic  cell  with 
built  in  LED,  light  sensitive  PIN  diode,  and  microphone;  2,  a 
Hewlett-Packard  model  3311A  function  generator;  3,  a  current 
amplifier;  4,  an  Ithaco  model  393  lock-in  amplifier;  and  5,  a 
Zenith  Z-150  personal  computer  with  installed  Qua-Tech  D-A  and  A- 
D  converters,  connected  to  a  Qua-Tech  input/output  terminal  con- 
nector board.  Additional  supporting  components  include  two  dc 
power  supplies  needed  to  power  the  current  amplifier  and  the 
microphone,  and  a  Keithley  Instruments  model  155  microvol tmeter 
used    to    monitor    the    light    intensity    in    the    photoacoustic   cell    as 
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FIGURE  3.1.   Schematic  diagram  of  the  pho toacous t ic 
sampling  apparatus. 
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detected  by  the  PIN  diode.  A  more  detailed  discussion  of  these 
components,  and  how  they  are  incorporated  into  the  actual  ar- 
rangement used  in  these  photoacoustic  sample  analysis  studies,  is 
given   below. 

3.1      The  Photoacoustic  Cell 

There  is  probably  no  single  photoacoustic  cell  design  which 
is  best  for  all  types  of  photoacoustic  sample  analysis.  The  cell 
must  be  able  to  accommodate  the  sample  being  analyzed,  the  modu- 
lated energy  beam  used,  the  acoustic  signal  detector,  and  any 
other  experimental  requirements.  Designs,  therefore,  will  vary 
according  to  these  constraints.  There  are,  however,  certain 
criteria  that  must  be  considered  In  the  design  in  order  to  maxi- 
mize the  acoustic  signal  obtained  from  the  sample  and  minimize 
all  extraneous  signals  and  noise.  These  include  the  size  of  the 
sample  chamber,  the  acoustic  isolation  of  the  cell,  the  material 
of  which  the  cell  is  constructed,  and  the  relative  placement  of 
the   microphone    In    the    cell. 

A.    Description   of    Photoacoustic    Cell 

A  diagram  of  the  cell  used  In  these  studies  is  given  in 
Figure  3.2.  The  main  structure  of  the  cell  was  constructed  from 
a  rectangular  block  of  aluminum,  having  a  height  of  7.6  cm  a 
length  of  10.2  cm  and  a  width  of  6.4  cm.  Each  of  the  two  6.4x7.6 
cm  sides  of  the  cell  had  rectangular-shaped  cavities  machined 
into  the  main  cell  structure.  Installed  into  the  cavity  on  the 
left  side  was   a  Hewlett-Packard    5082-4200    light-sensitive   PIN 
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diode,  and  Installed  into  the  right-side  cavity  was  a  Radio  Shack 
Archer  electret  microphone  element.  The  volume  in  these  cavities 
not  occupied  by  the  diode  and  microphone  was  packed  with  an 
insulating  filler  material.  The  cavities  were  sealed  by  1.25-cm- 
thick  lead  plates,  which  were  fastened  to  each  side  by  two  tight- 
ening screws.  A  thin  sheet  of  compressible  plastic  foam  was 
positioned  between  each  lead  plate  and  the  main  structure  of  the 
cell  to  improve  the  seal  at  the  aluminum-lead  interface. 

Centered  at  the  top  of  the  main  structure  of  the  cell  was  a 
cylindrical  cavity  having  a  3.8  cm  diameter,  which  extended  from 
the  top  surface  to  a  depth  of  3.0  cm.  Centered  at  the  base  of 
this  cavity  was  a  smaller  cylindrical  cavity  with  1.0  cm  diameter 
and  adjustable  depth.  This  portion  of  the  cell  was  the  actual 
sample  chamber.  The  base  of  this  chamber  consisted  of  a  separate 
cylinder  held  in  place  by  a  pair  of  0-rings  which  were  installed 
in  two  separate  circular  depressions  formed  about  the  circumfer- 
ence of  the  cylinder,  permitting  a  tight  fit  between  the  sample 
chamber  cavity  and  the  chamber  base.  By  applying  sufficient 
pressure,  the  cylindrical  base  could  be  forced  to  slide  up  or 
down  within  the  chamber  cavity,  thus  giving  the  sample  chamber  an 
adjustable  depth  and  variable  volume,  with  a  maximum  chamber 
depth  of  approximately  1.7  cm.  Surrounding  the  circumference  of 
this  chamber,  In  the  base  of  the  larger  cylindrical  cavity,  was  a 
circular  notch  Into  which  was  placed  a  rubber  0-ring.  On  the  top 
of  this  0-ring,  which  protruded  slightly  above  the  surface  of  the 
cavity  base,  was  centered  a  square,  1.57  millimeter  thick,  Supra- 
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sil  fused  quartz  window.  Gas  ports  were  drilled  in  the  front  and 
rear  sides  of  the  photoacoustic  cell,  extending  from  the  sample 
chamber  to  the  outside  of  the  cell.  These  ports  could  be  opened 
or  closed  using  needle  valves  installed  into  these  openings  on 
the  outside  of  the  cell's  main  structure.  Also,  a  port  extended 
from  the  sample  chamber  to  the  right  side  cavity  of  the  main 
structure,  where  it  met  with  the  microphone  installed  into  this 
cavity. 

Into  the  larger  cylindrical  cavity  extending  downward  from 
the  top  of  the  cell  was  fitted  a  removable  light  source  contained 
in  a  cylindrical  plug.  The  plug  was  constructed  of  aluminum, 
with  a  diameter  equal  to  that  of  the  cavity  and  length  slightly 
greater  then  the  cavity  depth.  Once  inserted  completely  into  the 
cavity,  this  cylindrical  plug  could  be  fastened  in  place  by 
tightening  a  plastic  set  screw  installed  into  the  front  side  of 
the  cell's  main  structure.  The  bottom  end  of  the  plug  (the  end 
placed  into  the  cavity)  had  a  shallow,  disk-shaped  recess.  The 
light  source  built  Into  this  plug,  which  protruded  slightly  from 
a  small  hole  at  the  center  of  this  recess,  was  a  Texas  Instru- 
ments TIL-31  high-intensity,  infrared-emitting  LED,  positioned  so 
that  it  irradiated  downward,  through  the  quartz  window  and  into 
the  sample  chamber.  The  recess  in  the  base  of  the  plug  was  small 
enough  in  diameter  so  that  when  the  plug  was  inserted  into  the 
cavity,  the  corners  of  the  quartz  window  encountered  the  outer, 
non-recessed  portion  of  the  plug  base,  thus  allowing  the  quartz 
window    to    be    pressed    up    against    the    rubber    O-ring,     with    pressure 
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applied  to  the  inserted  plug  being  transmitted  to  the  O-rlng 
seal.  In  the  side  of  the  plug,  near  Its  base,  extending  through 
to  the  recess  in  the  plug  base,  was  a  small  hole  allowing  for  a 
small  amount  of  the  scattered  radiation  emitted  by  the  LED  to 
escape.  Also,  In  the  left  side  of  the  main  structure  of  the 
photoacoustic  cell  was  a  hole  which  extended  to  the  left  side 
cavity  in  the  main  structure,  into  which  was  installed  the  PIN 
diode.  The  relative  position  of  these  two  holes  was  such  that, 
when  the  cylindrical  plug  was  properly  aligned  the  cavity,  the 
holes  would  be  aligned  as  well,  and  an  opening  would  completely 
extend  from  the  recess  in  the  plug  to  the  light  sensitive  PIN 
diode,  allowing  for  a  measure  of  the  relative  intensity  of  e-m 
radiation    produced    by    the    LED. 

In  order  to  monitor  the  output  of  the  microphone  and  PIN 
diode,  BNC  connectors  were  installed  on  the  left  and  right  sides 
of  the  top  surface  of  photoacoustic  cell's  main  structure  with 
electrical  connections  extending  to  the  PIN  diode  and  the  micro- 
phone respectively.  Also,  toward  the  left  side  of  the  rear  of 
the  main  structure,  a  BNC  connector  also  was  installed,  extending 
to  the  microphone,  so  as  to  provide  it  with  a  power  input.  A  BNC 
connector  was  also  positioned  on  the  top  surface  of  the  cylindri- 
cal plug,  with  electrical  connections  extending  to  the  LED,  in 
order    to    provide    the   LED   with   a    power    source. 

B.    Discussion   of    the   Photoacoustic    Cell    Design 
The  photoacoustic  cell   used   in   these   experiments   was   con- 
structed  of    aluminum    so    as    to    minimize    the   photoacoustic   response 
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produced  by  the  walls  of  the  sample  chamber.  Although  It  is  not 
feasible  to  completely  eliminate  the  chamber  response,  use  of  a 
material,  such  as  aluminum,  with  relatively  low  absorption  (high 
reflectivity),  high  density,  and  high  thermal  conductivity, 
should  give  rise  to  a  chamber  response  small  in  comparison  with 
that  of  the  sample  being  studied,  assuming  that  the  sample  does 
not  possess  similar  characteristics.  The  microphone  in  the  cell 
was  sufficiently  remote  from  the  illuminated  portion  of  the 
chamber  so  that  any  scattered  light  incident  upon  it  should  be 
Insignificant  in  comparison  to  the  direct  illumination  of  the 
sample;  thus  any  contribution  to  the  photoacoustic  signal  arising 
from  the  microphone's  interaction  with  this  stray  light  should  be 
minimal.  The  fused  quartz  window  used  in  the  cell  has  a  high 
degree  of  transparency,  allowing  the  sample  to  receive  as  much  of 
the  emitted  radiation  as  possible  from  the  LED,  and  absorbing 
very  little  so  as  to  minimize  the  response  of  the  window  itself. 

In  photoacoustic  examination  of  solids,  the  purpose  of  the 
gas  in  the  cell  is  to  undergo  pressure  fluctuations  in  response 
to  temperature  increases  at  the  surface  of  the  solid  sample,  and 
to  transmit  these  pressure  increases  to  the  microphone.  Thus,  it 
is  important  that  the  gas  used  in  the  cell  be  as  non-absorbing  as 
possible,  as  e-m  radiation  absorbed  in  the  gas  will  lead  to  a 
strong  undesirable  photoacoustic  response  from  the  gas  itself. 
The  gas  used  in  the  sample  chamber  was  helium,  which  was  intro- 
duced through  the  needle-valve-controlled  openings  in  the  front 
and  rear  sides  of  the  cell.   Since  helium  Is  a  non-absorbing  gas, 
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problems  with  absorption  and  subsequent  photoacoustic  responses 
from  the  gas  were  minimized.  Furthermore,  the  Rosencwaig-Gersho 
theory  states  that  the  amplitude  of  the  photoacoustic  response 
for  a  given  sample  increases  as  the  thermal  conductivity  of  the 
gas  used  in  the  chamber  is  increased,  varying  as  k1^2.  Since  it 
is  desirable  to  obtain  as  high  of  a  signal  amplitude  as  possible, 
it  is  beneficial  to  use  a  gas,  such  as  helium,  with  a  relatively 
high  thermal  conductivity. 

Another  concern  in  photoacoustic  cell  design  is  minimization 
of  the  effects  of  noise  from  the  outside  environment.  Since  a 
lock-in  amplifier  was  used  to  amplify  the  photoacoustic  signal, 
locking-in-on  and  amplifying  only  signals  having  the  same  fre- 
quency and  a  constant  phase  relationship  as  compared  to  the 
modulation  frequency  of  the  incident  e-m  radiation,  a  certain 
amount  of  noise  could  be  tolerated,  as  long  as  the  noise  was  not 
periodic  with  the  same  frequency  as  the  modulation  frequency 
being  used.  For  the  most  part,  external  noise  did  not  present 
any  serious  problem,  which  is  at  least  partially  due  to  the 
design  of  our  photoacoustic  cell.  Although  the  cell  was  primar- 
ily constructed  of  aluminum,  a  material  of  fairly  low  density, 
the  walls  of  the  cell  were  sufficiently  thick  for  adequate  acous- 
tic isolation  of  the  sample  chamber  from  outside  noise,  as  re- 
flected in  the  dimensions  of  the  cell  previously  given.  Also 
the  lead  plates  used  to  enclose  the  ends  of  the  cell  helped  to 
provide  better  acoustic  isolation  for  the  chamber,  and  the  entire 
photoacoustic  cell  was  placed  on  a  thick  lead  block  throughout 
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all  the  experiments  to  help  reduce  the  effects  of  external  vibra- 
tions. 

The  physical  size  of  the  sample  chamber  Is  also  an  Important 
consideration  In  the  design  of  a  photoacoustlc  cell.  As  previ- 
ously stated,  the  design  of  the  sample  chamber  was  such  that  the 
volume  of  the  chamber  In  the  photoacoustlc  cell  could  be  varied 
by  changing  the  the  chamber  depth.  The  volume  of  the  chamber  Is 
Important,  since  the  acoustic  signal  Is  Inversely  proportional  to 
the  length  of  the  gas  column  occupying  the  cell,  as  stated  previ- 
ously In  the  theoretical  section.  This  dependency  suggests  using 
as  small  a  volume  as  possible  In  order  to  maximize  the  amplitude 
of  the  response.  However,  according  to  the  Rosencwalg-Gersho 
theory,  It  Is  also  Important  that  the  dimensions  of  the  sample 
chamber  be  greater  than  one  thermal  diffusion  length  of  the  gas 
occupying  the  cell,  since  the  gas  layer  of  one  diffusion  length 
thickness  at  the  gas-sample  Interface  expands  to  create  the 
pressure  fluctuation  in  the  cell,  or  acoustic  signal,  which  is 
ultimately  transmitted  to  the  microphone.  In  addition,  the  cham- 
ber size  should  be  large  enough  so  that  dissipation  of  the  signal 
to  the  walls  and  window  of  the  cell  prior  to  the  detection  of  the 
signal  by  the  microphone  is  not  excessive. 

The  Rosencwaig-Gersho  theory  is  strictly  valid  only  for 
cases  where  the  length  of  the  gas  column  is  substantially  greater 
than  the  thermal  diffusion  length  of  the  gas  used  in  the  chamber. 
For  helium,  the  density,  specific  heat,  and  thermal  conductivity 
are  0.180  kg/cm3,  5232  J/kg-deg  C,  and  0.142  J/m-s-deg  C  respec- 

50 


tively  (21).  The  values  of  these  parameters  result  In  a  thermal 
diffusion  length  for  helium  of  1.73  cm,  at  a  modulation  frequency 
of  1  Hz,  which  corresponds  approximately  to  the  maximum  chamber 
depth  obtainable  for  the  photoacoustlc  cell  used.  Of  course,  as 
the  modulation  frequency  Is  Increased,  the  thermal  diffusion 
length,  which  varies  as«""2,  will  decrease,  becoming  0.055  cm 
at  a  modulation  frequency  of  1000  Hz.  It  Is  therefore  to  be 
expected  that  the  Rosencwalg-Gersho  theory  will  not  be  completely 
applicable  to  the  results  obtained  at  low  modulation  frequencies, 
and  becoming  even  less  so  as  the  chamber  depth  Is  decreased. 

The  size  effects  of  a  photoacoustlc  chamber  having  dimen- 
sions comparable  to  or  less  than  the  thermal  diffusion  length  of 
the  gas  used  In  the  sample  chamber  have  been  Investigated,  both 
theoretically  and  experimentally,  by  Aamodt,  Murphy,  and  Parker 
(22).  In  a  theoretical  derivation,  they  obtain  a  relationship 
giving  the  photoacoustlc  signal  amplitude  as  a  product  of  two 
separately  defined  parameters  A  and  g;  with  A  being  dependent 
upon  the  optical  absorption  coefficient,  thermal  properties,  and 
length  of  the  sample,  as  well  as  the  thermal  properties  of  the 
backing  material  for  optically  and  thermally  thin  samples;  and  g 
being  essentially  a  sensitivity  factor,  determined  by  such  para- 
meters as  the  cell  length  and  gas  thermal  properties,  being 
relatively  insensitive  to  the  sample's  thermal  properties  when 
the  cell  is  long  compared  to  the  thermal  diffusion  length  of  the 
gas.  In  the  case  where  the  the  length  of  the  gas  column  is  much 
greater  than  the  gas  thermal  diffusion  length,  their  results 
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parallel  the  conclusions  of  the  Rosencwalg-Gersho  theory,  with 
the  amplitude  of  the  photoacoustic  response  varying  inversely  as 
the  length  of  the  sample  chamber.  As  the  length,  L,  of  the  sam- 
ple chamber  is  decreased  and  approaches  the  thermal  diffusion 
length,  Ug,  of  the  gas,  the  amplitude  vs  chamber  length  curve 
levels  off  and  reaches  a  maximum  at  a  chamber  length  value  of 
approximately  2  =  1^.  For  H<ug,  the  amplitude  vs  chamber  length 
curve  generally  takes  on  a  positive  slope,  with  signal  amplitude 
decreasing  as  chamber  length  is  decreased.  In  this  region  the 
exponential  dependence  of  the  amplitude  of  the  signal  on  the 
sample  chamber  length  is  dependent  on  the  thermal  properties  of 
both  the  sample  and  chamber  window.  Figure  3.3  shows  some  of 
their  published  experimental  data,  for  both  the  amplitude  and  the 
phase  of  the  signal,  illustrating  this  effect.  The  sample  used 
in  their  experiment  was  carbon  black  and  the  gas  used  was  helium. 
The  maximum  amplitude  for  a  given  modulation  frequency  occurs  at 
approximately  the  point  at  which  the  chamber  length  is  equal  to 
the   thermal   diffusion   length   for   helium. 

3.2      The   Modulated   Electromagnetic  Radiation  Source 

The  source  of  e-m  radiation  used  throughout  all  experiments 
done  in  this  thesis  work  was  a  Texas  Instruments  TIL-31  light 
(infrared)  emitting  diode.  This  LED  was  built  into  the  photoa- 
coustic cell,  as  previously  described,  and  was  oriented  so  that 
e-m  radiation  emitted  from  the  diode  was  directed  downward, 
illuminating    the   sample    located   directly   below   It.      The   diode   had 
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FIGURE  3.3.  A  plot  of  the  variation  in  the  pho toacous t ic 
response  of  carbon  black  with  cell  length  at 
several  different  modulation  frequencies. 
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a  fixed  lens  positioned  on  Its  surface  to  help  focus  the  beam. 
The  half-intensity  beam  angle  of  the  diode  was  only  10  degrees, 
thus  a  minimal  amount  of  beam  spreading  was  experienced,  enabling 
the  majority  of  the  e-m  radiation  output  of  the  diode  to  be 
directed  to  the  sample  surface.  The  output  radiation  of  the  LED 
was  essentially  monochromatic,  having  a  peak  emission  wavelength 
of  9400  angstroms  with  a  typical  total  bandwidth  of  500  ang- 
stroms. The  on-off  and  off-on  switching  times  of  the  LED  were 
350  ns  and  600  ns  respectively,  which  were  more  than  adequate  for 
the   modulation    frequencies    used    In    these    experiments. 

Intensity  modulation  of  the  LED's  output  radiation  was 
achieved  by  modulation  of  Its  power  supply.  The  modulation  was 
in  the  form  of  an  on/off  square  wave,  with  the  signal  being 
supplied  by  a  Hewlett-Packard  model  3311A  variable  frequency 
function  generator.  This  function  generator  allowed  the  user  to 
select  a  frequency  range  via  a  front  panel  selector  switch  and 
then  vary  the  output  frequency  of  the  generator  throughout  this 
range  either  by  adjusting  a  frequency  calibrated  front  panel 
potentiometer,  or  by  sending  a  0  to  -15  volt  dc  signal  into  a 
voltage-controlled-oscillator  (VC0)  input  located  in  the  back  of 
the  instrument  (0  volts  correspond-ing  to  lowest  frequency  of 
selected  range).  The  frequency  output  of  the  instrument  with 
respect  to  a  given  voltage  input  was  designed  to  be  linear 
throughout  the  selected  frequency  range.  A  factor  of  10  frequen- 
cy increase  from  the  lowest  frequency  of  a  given  range  was  per- 
mitted,   using   either    of    these    two   methods,    before    an    increase    in 
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the  frequency  range  of  the  instrument  needed  to  be  made.  A 
current  amplifier  was  used  to  amplify  the  signal  supplied  to  the 
LED  by  the  frequency  generator,  drawing  power  from  a  Lambda  model 
LP  410  FM  dc  power  supply. 

A  measure  of  the  relative  LED  intensity  as  detected  by  the 
PIN  diode  installed  in  the  pho toacous tic  cell  was  provided  by 
connection  of  the  PIN  diode  output  to  a  Keithley  Instruments 
model  155  mlcrovol tme ter.  A  front  panel  meter  with  variable 
range  selector  was  provided  by  the  Instrument  for  monitoring  of 
the  input  voltage.  In  addition,  the  microvol tmeter  also  provided 
a  voltage  output,  ranging  from  0  to  1  volt,  directly  proportional 
to  the  input  voltage,  with  a  proportionality  constant  equal  to 
the  selected  range  setting.  Since  at  low  modulation  frequencies, 
a  stable  voltage  reading  Initially  could  not  be  obtained  due  to 
the  periodicity  of  the  light  source  being  detected,  a  resistor- 
capacitor  circuit  was  connected  in  parallel  with  the  microvoltme- 
ter  to  rectify  the  output  voltage  from  the  diode. 

3.3   The  Acoustic  Signal  Detection  Apparatus 

In  photoacoustic  sample  analyses  which  utilize  a  gas  trans- 
fer medium  for  the  acoustic  signal,  the  most  commonly  used  device 
to  monitor  the  acoustic  signal  is  a  sensitive  microphone,  either 
flat  or  cylindrical-surfaced.  Although  cylindrical  microphones 
are  generally  more  sensitive  then  flat  microphones  due  to  their 
larger  surface  area,  the  acoustic  response  of  a  cylindrical 
microphone  typically  varies  with  the  frequency  of  the  signal  that 
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it  monitors,  while  a  flat  microphone  typically  maintains  a  flat 
response  over  a  wide  frequency  range.  A  flat  frequency  response 
is  more  desirable  for  depth  profiling  studies  where  the  modula- 
tion frequency  is  varied  over  a  wide  frequency  range.  Conse- 
quently, since  a  frequency-independent  response  was  desired  in 
these  experiments,  a  flat-surfaced  microphone  was  used. 

The  microphone  incorporated  into  the  photoacous tic  cell  was 
a  Radio  Shack  Archer  electret  microphone  element  with  a  built-in 
preamplifier.  It  had  a  flat  frequency  response  throughout  a  20- 
3000  Hz  frequency  range,  with  the  frequency  response  dropping  off 
somewhat  at  frequencies  lower  than  20  Hz.  Throughout  this  exper- 
imentation, the  microphone  proved  to  be  sufficiently  sensitive  to 
detect  the  acoustic  signals.  The  microphone  preamplifier  was 
powered  by  a  Lambda  model  LL-903  variable  voltage  power  supply, 
maintained  at  a  constant  voltage  of  8  volts.  The  preamplifier 
output  was  sent  directly  to  an  Ithaco  model  393  lock-In  ampli- 
fier. 

This  lock-in  amplifier  was  capable  of  accepting  a  periodic 
external  reference  signal  within  a  limited  frequency  range,  and 
amplifying  from  an  input  signal  only  the  ac  components  with  fre- 
quencies Identical  to  the  reference.  Modification  of  the  range 
of  reference  frequencies  which  could  be  accepted  by  the  lock-in 
amplifier  was  possible,  being  achieved  by  replacement  of  several 
removable  Internal  component  boards.  In  these  experiments,  two 
sets  of  these  component  boards  were  used;  one  set  enabling  ampli- 
fication of  input  signals  within  a  frequency  range  of  1-100  Hz, 
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and  the  other  within  a  10-1000  Hz  frequency  range.  For  a  given 
set  of  installed  boards,  amplification  of  signals  at  frequencies 
higher  than  a  factor  of  ten  increment  above  their  lowest  accepta- 
ble frequency  required  changing  the  setting  of  a  toggle  switch 
located  on  the  front  panel  of  the  lock-in  amplifier.  Also  lo- 
cated on  the  front  panel  was  a  sensitivity  selector,  which  needed 
to  be  positioned  to  an  appropriate  sensitivity  setting,  with 
respect  to  the  signal  being  monitored.  The  amplitude,  in-phase 
component,  and  out-of-phase  component  of  the  amplified  signal 
could  be  monitored  either  by  viewing  a  selectable  front  panel 
meter  built  into  the  instrument,  or  by  measuring  voltages  sup- 
plied by  three  BNC  outputs  located  in  the  back  of  the  instrument. 
These  voltage  outputs  were  proportional  to  the  amplified  signal, 
with  the  proportionality  constant  being  essentially  equivalent  to 
the   sensitivity   setting   of    the   amplifier. 

In  order  to  achieve  amplification  of  the  true  photoacoustlc 
response,  the  signal  used  to  modulate  the  LED  was  also  used  as  a 
reference  input  for  the  lock-in  amplifier.  Since  the  photoacous- 
tlc signal  is  of  the  same  frequency  as  that  of  the  modulated 
light  source,  this  allowed  the  lock-in  amplifier  to  amplify  only 
the  photoacoustlc  signal,  and  and  essentially  disregard  all  am- 
bient noise  In  the  chamber  detected  by  the  microphone.  One 
source  of  noise  with  distinct  periodicity  which  needed  to  be 
minimized  was  inductive  pickup  from  the  60  Hz  ac  line  voltage  and 
ripple    in    the    power    supply    for    the    microphone    preamplifier. 
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Modulation   frequencies   close    to   60   Hz,    and   multiples    thereof, 
were    therefore    avoided. 

3.4      The    Data   Acquisition    System 

In  acquiring  the  data  with  the  above-described  apparatus,  a 
software  driven  data  acquisition  system  was  used.  The  components 
of  the  acquisition  system  consisted  of  a  Zenith  Z-150  microcompu- 
ter, used  in  conjunction  with  a  Qua  Tech  PXB-721  parallel  expan- 
sion board,  with  added  ADM8-10  analog-to-digital  (A-D)  and  DM8-10 
digital-to-analog  (D-A)  converter  modules.  The  expansion  board, 
with  converter  modules,  was  installed  into  one  of  the  computer's 
expansion  slots,  with  ribbon  cables  extending  from  these  to  the 
exterior  of  the  computer.  The  ribbon  cables  were  attached  to  a 
Qua  Tech  UIO-10  universal  I/O  board  which  permitted  connecting  of 
the  installed  expansion  board  and  converter  modules  to  the  ex- 
perimental apparatus.  A  capacity  for  eight  inputs  to  the  A-D 
converter  and  eight  outputs  from  the  D-A  converter  were  provided. 
The  analog  voltage  output  from  the  D-A  converter  could  be  varied 
from  -5  to  5  volts  in  increments  of  0.03922  volts.  The  analog  to 
digital  converter  was  configured  to  accept  analog  voltage  inputs 
from  0  to  0.5  volts,  with  a  resolution  of  8  bits,  or  1.953 
millivolts. 

The  software  used  to  control  the  operation  of  the  A-D  and  D- 
A  converters  was  written  in  Microsoft  BASICA,  and  used  in  con- 
junction with  a  software  package  provided  by  Qua  Tech,  enabling 
access    to    the   converter  modules   within    the   context   of   a   BASIC 
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program.   This  software  package,  entitled  LABSTAR,  defined  a 
series  of  BASIC  commands,  allowing  various  modes  of  access  to  the 
D-A  and  A-D  modules  via  a  machine  language  program,  LABSTAR.EXE; 
which  needed  to  be  separately  loaded  into  computer  memory  prior 
to  loading  the  Microsoft  BASIC  interpreter.   The  actual  BASIC 
application  program  utilizing  these  specialized  commands  required 
Initial  execution  of  a  small  BASIC  program  segment,  provided  by 
Qua-Tech  under  the  filename  LABSTAR. BAS,  designed  to  compose  the 
first  99  lines  of  the  application  program.  The  purpose  of  this 
segment  of  the  application  program  was  to  define  commands  to  be 
used  within  the  program  to  invoke  the  machine  language  subrou- 
tines.  Prior  to  execution  of  one  of  these  specialized  commands, 
a  CALL  statement  specifying  memory  address  information  regarding 
the  appropriate  machine  language  subroutine  first  needed  to  be 
executed,  and  any  user-assigned  variables  associated  with  the 
command  needed  to  be  properly  defined.   Some  programming  capabil- 
ities for  A-D  conversion  provided  by  these  commands  included  the 
ability  to  rapidly  acquire  blocks  of  data,  at  a  maximum  sampling 
rate  of  one  data  point  every  33  microseconds,  either  from  a 
single  channel  or  from  a  specified  range  of  channels.   Also 
possible  was  the  acquiring  of  either  a  single  point  or  block  of 
data  upon  the  request  of  an  external  trigger  signal  or  an  exter- 
nal clock.   Similar  commands  associated  with  the  D-A  converter 
were  capable  of  specifying  either  output  of  single-valued  volt- 
ages from  a  given  channel,  maintained  until  a  different  voltage 
at  that  channel  is  specified,  or  block  voltage  outputs  specified 
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within  an  array.  For  block  voltage  outputs,  voltages  could  be 
rapidly  varied  from  a  single  channel,  or  sequentially  from  a 
range  of  channels,  at  a  maximum  output  rate  of  42  microseconds. 
Also  single  point  or  block  data  output  could  be  initiated  by  way 
of  request  from  an  external  trigger.  A  more  detailed  listing  and 
description  of  allowed  commands  is  provided  in  a  user's  manual 
supplied    by    Qua-Tech. 

The  primary  tasks  of  the  data  acquisition  program  were  to 
scan  through  a  range  of  modulation  frequencies,  achieved  by 
incrementing  the  frequency  output  of  the  waveform  generator,  and 
obtain  and  record  the  pho toacous tic  response  from  the  lock-in 
amplifier  for  each  of  the  modulation  frequencies  produced.  In- 
crementation of  the  modulation  frequency  was  achieved  by  incre- 
mentally varying  a  dc  voltage  sent  to  the  VCO  inputs  of  the 
waveform  generator,  provided  by  a  channel  from  the  D-A  converter. 
Since  the  D-A  converter  was  capable  of  providing  negative  voltage 
outputs  within  a  range  of  only  0  to  -5  volts,  the  waveform  gener- 
ator was  modified  so  that  the  entire  frequency  range  for  a  given 
front  panel  range  setting  of  the  instrument  could  be  scanned 
within  this  voltage  range,  instead  of  its  normal  0  to  -15  volt 
range.  This  modification  affected  somewhat  the  linearity  of  the 
output  frequency  vs  input  voltage  curve,  therefore,  a  best-fit 
polynomial  function  was  found  to  relate  the  modulation  frequency 
to  a  given  input  voltage.  The  frequency  values  provided  by  this 
function  were  accurate  to  approximately  +/-  1%,  with  the  devia- 
tion being   consistent   from   run    to   run.      The   amplitude   and    the   in- 
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phase  component  of  the  photoacoustic  response,  as  provided  by  the 
lock-in  amplifier,  were  accessed  by  the  computer  via  two  separate 
channels  of  the  installed  A-D  converter.  A  default  number  of 
fifty  consecutive  samplings  were  taken  from  the  lock-in  amplifier 
at  evenly  spaced  intervals,  over  a  length  of  time  equal  to  the 
time  constant  setting  of  the  lock  in  amplifier,  for  both  of  these 
parameters  at  a  given  modulation  frequency.  The  fifty  values  for 
each  parameter  were  subsequently  averaged  to  give  the  value  of 
that  parameter  corresponding  to  the  current  modulation  frequency. 
Since  the  A-D  converter  was  configured  only  to  accept  input 
voltages  from  0  to  0.5  volts,  the  lock-in  amplifier  was  used  so 
as    to   provide   output   voltages   below   this    0.5    volt   constraint. 

Since  the  amplitude  of  the  acoustic  signal  decreases  as  the 
modulation  frequency  of  the  incident  radiation  increases,  manual 
changes  in  the  sensitivity  setting  of  the  lock-in  amplifier 
needed  to  be  made  periodically  in  the  course  of  analysis  of  a 
given  sample.  When  the  amplitude  of  the  photoacoustic  response 
became  low  enough  that  an  increase  in  the  sensitivity  setting  of 
the  lock-in  amplifier  could  be  achieved  without  producing  an 
output  signal  to  the  A-D  converter  exceeding  0.5  volts,  the 
execution  of  the  program  was  temporarily  stopped  to  allow  for 
such  a  sensitivity  increase.  Once  the  setting  change  was  made, 
the  program  could  be  allowed  to  continue.  Another  point  in  the 
sampling  procedure  that  required  manual  setting  changes  was  the 
reaching  of  the  upper  end  of  the  frequency  range  for  a  given 
range    setting    of    the    waveform    generator.       At    this    point    the 
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execution  of  the  program  needed  to  be  temporarily  terminated  to 
allow  for  increasing  this  frequency  range.  Since  the  upper  end 
of  the  input  frequency  range  for  the  lock-in  amplifier  was 
reached  at  this  point  as  well,  a  manual  setting  adjustment  to 
increase  the  acceptable  frequency  range  needed  to  be  made  on  the 
lock-in  amplifier  at  this  point  also.  Due  to  the  power  law 
dependence  of  the  signal  amplitude  on  the  modulation  frequency, 
fewer  data  points  were  needed  at  the  higher  modulation  frequen- 
cies. Therefore,  periodically  (at  multiples  of  250  Hz),  the  size 
of  the  modulation  frequency  increment  was  allowed  to  increase  to 
help   increase    the   speed   of    the   sampling   process. 

Another  experimental  parameter  monitored  by  the  A-D  conver- 
ter was  the  PIN  diode  voltage,  as  detected  by  the  microvol tmeter 
to  which  it  was  connected.  In  this  instance,  the  0  to  1  volt 
voltage  output  range  of  the  instrument  was  reduced  to  a  0  to  0.5 
voltage  range  using  a  2-1  voltage  divider.  The  output  voltage 
from  the  raicrovoltmeter  within  this  range  was  sent  to  a  channel 
of  the  A-D  converter  from  which  a  series  of  fifty  readings  were 
taken  over  a  one  second  interval,  and  subsequently  averaged,  for 
each  modulation  frequency.  Due  to  the  fact  that  the  intensity  of 
the  LED  did  not  vary  greatly  throughout  the  entire  range  of 
modulation  frequencies  experimentally  used,  modification  of  the 
range  selection  of  the  microvol  tmeter  was  not  required  in  the 
course  of  a  given  run,  provided  that  a  sufficiently  low  range 
selection   had    initially   been   chosen. 
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When  scanning  through  the  desired  range  of  modulation  fre- 
quencies was  completed  for  a  given  sample  (1-50  Hz  for  the  low 
frequency,  and  10-1000  Hz  for  the  high  frequency  lock-in  ampli- 
fier component  boards),  sampling  was  terminated  and  the  acquired 
data  was  filed  on  5.25  inch  floppy  disks  for  future  analysis.  A 
hard  copy  of  the  data  was  also  printed  as  acquired  during  the 
actual  experimental  run.  Given  below  is  a  general  outline  sum- 
marizing the  major  objectives  used  in  programming  the  Zenith  Z- 
150  microcomputer  for  the  task  of  data  acquisition.  Two  slightly 
different  versions  of  the  program  were  created,  one  for  data 
acquisition  in  the  1-50  Hz  frequency  range,  and  one  for  the  10- 
1000  Hz  frequency  range,  since  both  ranges  could  not  be  scanned 
at  once  due  to  the  frequency  range  limitations  of  the  lock-in 
amplifier.  However,  both  programs  followed  this  general  outline 
(with  the  exception  of  step  10  and  part  of  step  5  being  omitted 
in  the  lower  frequency  range  acquisition  program).  A  generalized 
flowchart  of  the  data  acquisition  routine  utilized  in  these 
programs  (steps  5  through  14)  is  provided  in  Appendix  I. 

Data  Acquisition  Programming  Objectives 

1.  Initially  define  Labstar  commands  to  be  used  with  the  D- 
A  and  A-D  converters  (this  portion  of  the  program,  lines  1-99, 
provided  by  Qua-Tech). 

2.  Initialize  0  volts  to  VCO  input  of  waveform  generator. 

3.  Calibrate  waveform  generator  frequency  output  with  fre- 
quency counter  using  front  panel  potentiometer  so  that  0  volts 
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input    results    in    lowest    frequency    output    for    selected    frequency 
range. 

4.  Request  and/or  initialize  required  parameters  for  data 
acquisition    process    (increment    size,     sampling    rate,     etc.). 

5.  Compute  modulation  frequency  from  input  voltage  to  wave- 
form generator,  and  check  to  see  if  it  is  a  multiple  of  60  Hz;  If 
so,  then  increase  modulation  frequency  by  appropriate  increment 
and  continue  execution  of  program,  otherwise,  increase  data 
counter  by   one   and   store   value   in   data   matrix. 

6.  Execute  a  delay  loop  of  length  approximately  equal  to  4 
times  the  time  constant  setting  on  the  lock-in  amplifier,  to 
allow   lock-in   amplifier    to    reach   equilibrium. 

7.  Acquire  fifty  samplings  of  the  relative  LED  signal 
intensity  as  provided  from  the  micro  vo  1  tme  ter's  measurement  of 
the  PIN  diode  voltage,  average  these  values,  and  store  the  aver- 
age   value    in    the    data    matrix. 

8.  Acquire  from  the  lock-in  amplifier  numerous  samplings  of 
the  photoacoustic  signal  amplitude  and  the  in-phase  component  of 
the  amplitude  via  the  A-D  converter,  average  these  values,  and 
store   the  averaged   values    in    the   data   matrix. 

9.  Compute  the  phase  angle  from  the  amplitude  and  in-phase 
component  of  the  photoacoustic  response,  and  store  this  value  in 
the    data    matrix. 

10.  Check  to  see  if  the  size  of  the  frequency  increment 
should    be    increased;     If    so,     increase    it    by    the    prescribed    amount. 

11.  Check    to    see    if    the    sensitivity    setting    of    the    lock-in 
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amplifier   should   be    increased;    if   so,     temporarily  halt   execution 
of    the   program   to   allow   for    the   setting   adjustment. 

12.  Check  to  see  if  the  maximum  frequency  has  been  reached 
for  the  current  range  setting  of  the  waveform  generator;  if  it 
has  not,  then  increase  the  modulation  frequency  via  the  D-A 
converter   and    continue    execution    of    program   at    step    4. 

13.  Check  to  see  if  the  current  modulation  frequency  has 
reached  the  desired  termination  point;  if  it  has  not,  then  tempo- 
rarily halt  execution  of  the  program  to  allow  for  increasing  the 
frequency  range  setting  of  the  waveform  generator,  reinitialize 
appropriate  parameters  for  continued  program  execution  at  this 
new   setting,    and   continue   program   execution   at   step    4. 

14.  Print  out  the  important  parameter  values  used  in  acqui- 
ring the  data,  provide  the  option  for  saving  the  data  matrix  as  a 
data  file  on  floppy  disk  for  future  analysis,  and  terminate 
program   execution. 

Complete  listings  of  the  two  acquisition  programs  actually 
used  in  these  experiments,  both  for  the  low  modulation  frequency 
(1-50  Hz)  and  high  modulation  frequency  (10-1000  Hz)  ranges,  are 
provided    in   Appendices    II    and    III. 

Once  the  data  were  acquired  for  a  given  run  and  saved  as  a 
data  file  on  floppy  disk,  the  Lotus  1-2-3  (23)  spreadsheet  and 
graphics  software  was  used  manipulate  the  data  and  display  the 
results.  Data  files  saved  by  the  acquisition  software  were 
organized   so   that   they  could  be  easily  read   into  Lotus    1-2-3  by 


65 


using  its  FILE  IMPORT  command.  Once  the  data  were  read  into  a 
Lotus  spreadsheet  file,  procedures  such  as  normalization  and  data 
averaging  could  be  experimented  with  and  performed,  and  graphic 
displays  of  the  data  could  be  easily  generated  for  comparisons 
between    kernels. 
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CHAPTER  IV 


RESULTS  AND  DISCUSSION 


The  results  of  the  various  experiments  performed  throughout 
this  thesis  work  will  be  presented  and  discussed  in  this  chapter. 
Although  the  ultimate  goal  of  these  experiments  was  to  study  the 
photoacous tic  response  of  various  varieties  of  wheat  kernels,  a 
number  of  other  studies  were  also  done  to  determine  the  validity 
and  reproducibility  of  the  results  and  to  examine  the  nature  of 
the  photoacoustic  response  for  the  empty  sample  chamber  and  for 
various  chamber  sizes. 

At  a  given  modulation  frequency,  theory  states  that  the 
photoacoustic  response  of  a  sample  is  dependent  upon  the  sample's 
thermal  and  optical  properties.  However,  it  also  predicts  a 
dependence  upon  such  parameters  as  the  incident  e-m  radiation 
intensity,  the  ambient  pressure,  thermal  diffuslvity,  and  length 
of  the  gas  column  contained  In  the  cell,  and  the  ambient  tempera- 
ture of  the  sample  and  its  surroundings.  In  order  to  attribute 
variations  observed  In  the  photoacoustic  response  to  actual  vari- 
ations in  the  samples  tested,  it  Is  imperative  that  parameters 
not  associated  with  the  sample  itself,  whose  alteration  could 
affect  the  response,  be  maintained  at  constant  values  between 
runs  for  a  given  experiment.  The  variable  parameters  capable  of 
influencing  the  photoacoustic  response  in  these  experiments,  and 
the  values  used  and/or  procedures  required  to  initialize  them, 
are  as  follows: 
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1.  The  waveform  generator  was  operated  at  maximum  ampli- 
tude, with  0  volts  dc  offset  throughout  all  experiments,  and  was 
checked  for  proper  frequency  calibration  with  0  volt  VCO  input 
prior  to  each  run, 

2.  Settings  on  the  current  amplifier  and  power  supply 
powering  the  amplifier  were  maintained  so  that  at  a  modulation 
frequency  of  10  Hz  the  voltage  output  of  the  power  supply  was  3.1 
volts  and  the  average  current  output  was  0.19  amps. 

3.  The  microphone  voltage  input  was  maintained  at  a  con- 
stant value  of  approximately  8  volts. 

4.  The  lock-in  amplifier  was  operated  with  zero  suppress  in 
off  position,  signal  mode  in  maximum  noise  rejection  position, 
and  phase  adjust  set  at  0  degrees.  Also,  the  ground/ iso 1  a te 
toggle  was  set  for  isolated  signal  input  and  an  external  refer- 
ence input  was  provided  for  all  experiments. 

5.  The  sample  chamber  was  provided  with  a  helium  gas  atmos- 
phere at  a  pressure  approximately  equal  to  atmospheric  pressure. 
In  order  to  achieve  this,  the  following  procedure  was  followed. 
For  each  run,  the  input  and  output  needle  valves  were  opened 
approximately  two  full  turns,  and  a  mild  flow  of  helium  gas  was 
allowed  to  pass  through  the  chamber  for  five  minutes.  After  five 
minutes  had  elapsed,  both  valves  were  closed  nearly  simultaneous- 
ly, with  the  output  valve  being  closed  approximately  one  second 
after  the  input  valve  to  allow  for  pressure  equalization  with  the 
outside  environment. 

6.  The  ambient  temperature  of  the  chamber  and  the  samples 
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tested   was    room    temperature,    which   was   essentially   constant. 

7.  For  each  run,  the  sample  being  examined  was  centered  as 
closely  as  possible  in  the  chamber  to  allow  for  direct  illumina- 
tion by   the   LED. 

8.  For  experiments  concerned  with  detection  of  variations 
between  wheat  varieties,  the  cell  volume  was  left  unaltered 
between  runs. 

9.  The  sample  chamber  was  kept  closed  between  runs,  and  was 
periodically  cleaned  with  ethyl  alcohol,  as  was  the  cell  window, 
to   prevent   accumulation   of   contaminants   within   the   cell. 

The  experiments  were  performed  using  the  apparatus  described 
in  the  previous  chapter,  with  the  data  acquisition  process  being 
driven  by  the  BASIC  programs  listed  in  the  appendix.  The  quanti- 
ties measured  were  the  amplitude  and  in-phase  component  of  the 
amplitude,  from  which  the  phase  lag  of  the  photoacoustic  response 
was  also  derived.  Also  measured  for  many  of  the  experiments  was 
the  relative  intensity  of  the  LED.  These  quantities  were  all 
monitored  as  a  function  of  the  modulation  frequency  of  the  LED, 
throughout  overlapping  frequency  ranges  of  1-50  Hz  and  10-1000 
Hz,  with  separate  runs  required  to  scan  each  of  these  frequency 
ranges  for  a  given  sample,  due  to  the  earlier-mentioned  input 
frequency    limitations    of    the    lock-in  amplifier. 

4.1      Sample    Chamber    Studies 
A.    Black    Cloth    Study 
As    a    test    of    the    photoacoustic    cell    and    supporting    appara- 
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tus,  the  pho toacous tic  response  of  a  black  cloth  sample  was 
measured  as  a  function  of  modulation  frequency,  throughout  a 
frequency  range  of  1-1000  Hz.  The  cloth  sample  was  cut  in  the 
shape  of  a  circular  disk,  and  sized  so  as  to  completely  cover  the 
bottom  of  the  sample  chamber.  The  depth  of  the  sample  chamber 
for  this  experiment  was  4.8  mm.  The  log  of  the  amplitude  and  the 
phase  lag  of  the  response  were  then  plotted  as  a  function  of  the 
log  of  the  modulation  frequency,  and  are  shown  in  Figures  4.1  and 
4.2  respectively.  The  slight  vertical  shift  observed  in  the 
data,  which  is  most  evident  in  the  overlapping  10-50  Hz  frequency 
range,  is  a  result  of  a  variation  in  the  signal  processing  elec- 
tronics of  the  lock-in  amplifier  caused  by  replacing  its  set  of 
removable   internal   component  boards. 

The  black  cloth  sample  used  in  this  experiment  was  optically 
opaque,  and  it  would  be  reasonable  to  expect  that  its  optical 
absorption  coefficient  was  relatively  large,  resulting  in  a  rela- 
tively short  optical  absorption  length.  If  the  optical  absorp- 
tion length  is  much  less  then  the  thermal  diffusion  length  of  the 
sample  for  the  range  of  modulation  frequencies  used  in  this 
experiment,  which  would  be  quite  likely  in  this  instance,  the 
Rosencwaig-Gersho  theory  predicts  that  the  amplitude  of  the  pho- 
toacoustic  response  should  vary  as  u-1.  Upon  evaluating  the  data 
for  this  experiment,  such  a  dependence  Is  clearly  observed,  with 
the  exception  of  the  response  at  frequencies  below  approximately 
15  Hz,  where  the  slope  of  the  log-log  plot  becomes  less  negative. 
At  a  modulation   frequency   of    approximately    13   Hz,    the    thermal 
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Figure  4.1.   A  log-log  plot  of  the  amplitude  of  the 
photoacous tic  response  vs  modulation  frequency 
for  the  black  cloth  sample.   Sample  chamber 
depth  =  4.8  mm. 
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Fig  4.2.   A  plot  of  the  phase  of  the  photoacoustlc 
response  vs  log  modulation  frequency  for  the 
black  cloth  sample.   Sample  chamber  depth  =  4.8  mm. 
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diffusion  length  for  helium,  which  is  frequency  dependent,  be- 
comes equivalent  to  the  sample  chamber  length  used  for  this 
experiment.  A  plausible  explanation  for  this  behavior  is  that  at 
these  low  modulation  frequencies,  the  Rosencwaig-Gersho  theory 
becomes  invalid,  as  it  assumes  a  gas  column  of  length  much  great- 
er than  the  thermal  diffusion  length  of  the  gas.  The  behavior 
is,  however,  in  agreement  with  the  previously  discussed  findings 
of  Aamodt,  Murphy,  and  Parker  (22),  who  predict  a  decrease  in  the 
amplitude  of  the  pho toacous tic  signal  for  cases  where  the  size  of 
the  sample  chamber  becomes  comparable  to,  or  less  than,  the 
thermal  diffusion  length  of  the  gas.  Although,  in  our  case,  the 
relative  amplitude  of  the  signal  still  increases  as  the  thermal 
diffusion  length  of  the  gas  is  increased  at  low  modulation  fre- 
quencies, the  amplitude  of  the  signal  at  a  given  modulation 
frequency  is  decreased  from  what  it  would  potentially  be  if  a 
larger   chamber   volume   were    used. 

B.    Empty   Chamber    Studies 

Although  the  sample  chamber  was  designed  to  minimize  Its  own 
empty  chamber  signal,  it  is  reasonable  to  expect  that  the  chamber 
walls  will  exhibit  a  measurable  pho toacous tic  response,  as  a 
relatively  large  surface  area  of  solid  material  is  being  illumi- 
nated. Even  with  a  wheat  kernel  sample  present  in  the  chamber, 
the  cell  walls  may  contribute  significantly  to  the  resultant 
photoacoustic  response,  due  to  the  limited  size  of  the  kernel  and 
scattered  radiation  reflected  from  its  surface.  Thus,  it  Is 
Important    to    have    a    knowledge    of    the    photoacoustic    response    from 
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the  empty  chamber  itself. 

Figures  4.3  and  4.4  show  the  amplitude  and  phase  of  the 
empty  chamber  response  as  compared  to  the  signals  measured  with 
the  black  cloth  sample  present.  These  results  show  that  the 
amplitude  of  the  empty  chamber  signal,  although  certainly  measur- 
able, is  more  than  one  order  of  magnitude  lower  then  the  black 
cloth  signal,  and  that  the  phase  lag  vs  modulation  frequency 
curves  are  distinctly  different  for  the  empty  chamber  and  black 
cloth  sample.  This  result  is  not  surprising  in  that  the  black 
cloth,  being  a  very  strong  absorber  and  having  a  fairly  large 
surface  area,  should  be  expected  to  produce  a  relatively  strong 
photoacoustic  response;  whereas  the  empty  chamber,  whose  design 
was  aimed  at  minimizing  its  own  response,  should  produce  a  signal 
small  in  comparison.  It  is  also  seen  that  the  phase  of  the  black 
cloth  signal  is  very  much  different  empty  chamber  response  at 
frequencies  greater  than  approximately  5  Hz.  In  the  empty  cham- 
ber phase  angle  vs  frequency  curve,  there  is  seen  a  fairly  large 
difference  between  values  obtained  for  the  low-frequency  and 
high-frequency  runs.  It  is  our  contention  that  this  discrepancy 
was  due  to  a  variation  in  the  manner  with  which  the  complex  wave 
forms  detected  by  the  microphone  were  interpreted  by  the  lock-In 
amplifier,  caused  by  the  exchanging  of  the  internal  component 
boards  of  the  amplifier  between  these  runs. 

To  make  a  determination  of  the  empty  chamber  response  as 
compared  to  the  response  obtained  from  wheat  kernel  samples,  a 
series  of  tests  were  performed  in  which  the  photoacoustic  re- 
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Figure  4.3.   A  log-log  comparison  of  the  amplitude  of 

the  photoacoustlc  response  vs  modulation  frequency 
obtained  with  the  black  cloth  sample  present  to 
that  obtained  from  the  empty  chamber.   Sample  cham- 
ber depth  =  4.8  mm. 

Legend : 

square  -  black  cloth  sample  present 
plus  -  no  sample  present  (empty  chamber) 
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Figure  4.4.   A  comparison  of  the  phase  of  the  photo- 
acoustic  response  vs  log  modulation  frequency  ob- 
tained with  the  black  cloth  sample  present  to  that 
obtained  from  the  empty  chamber.   Sample  chamber 
depth  -  4.8  mm. 

Legend: 

square  -  black  cloth  sample  present 
plus  -  no  sample  present  (empty  chamber) 
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sponse  was  monitored  as  a  function  of  modulation  frequency  for 
the  empty  chamber,  a  single  kernel,  and  five  kernels  (five  being 
the  maximum  number  of  normal-sized  kernels  which  could  be  placed 
on  the  floor  of  the  chamber).  In  the  case  of  the  single  kernel, 
the  kernel  was  placed  crease-up  at  the  center  of  the  floor  of  the 
sample  chamber.  For  five  kernels,  one  kernel  was  centered  in  the 
chamber  with  the  other  four  kernels  placed  around  it.  The  re- 
sults of  these  studies,  for  both  the  amplitude  and  phase  of  the 
photoacoustic  response,  are  shown  in  Figures  4.5  and  4.6  respec- 
tively. It  can  be  seen  from  these  results  that  the  measured 
amplitude  of  the  photoacoustic  response  when  a  sample  was  present 
in  the  chamber  was  slightly  greater  then  the  empty  chamber  re- 
sponse throughout  most  of  the  frequency  range,  with  the  response 
being  greatest  for  five  kernels.  The  phase  angle  measurements 
for  these  runs  reflect  a  similar  trend,  with  the  phase  lag  being 
smaller  for  the  empty  chamber,  except  for  within  approximately 
the    100-500    Hz    frequency    range. 

Upon  inspection  of  the  sample  chamber,  we  decided  that  the 
protruding  portion  of  the  0-ring  sealing  the  window  to  the  sample 
chamber  was  possibly  receiving  a  sufficient  amount  of  e-m  radia- 
tion to  cause  It  to  contribute  significantly  to  the  measured 
photoacoustic  response  of  the  chamber.  To  reduce  this  problem, 
and  to  increase  the  percentage  of  the  incident  radiation  striking 
the  sample,  we  reduced  the  0-ring  thickness  so  that  It  barely 
protruded  above  the  plug  cavity  base,  and  Installed  a  3  mm  diame- 
ter  aperture    at    the    base    of    the   cylindrical    plug   containing    the 
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Figure  4.5.   A  log-log  comparison  of  the  amplitude  of 

the  photoacous tic  response  vs  modulation  frequency 
obtained  with  a  single  kernel  and  five  kernels  of 
wheat  present  to  that  obtained  from  the  empty  cham- 
ber.  Sample  chamber  depth  =  4.8  mm. 

Legend: 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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Figure  4.6.   A  comparison  of  the  phase  of  the  photo- 
acoustic  response  vs  log  modulation  frequency  ob- 
tained with  a  single  kernel  and  five  kernels  of 
wheat  present  to  that  obtained  from  the  empty  cham- 
ber.  Sample  chamber  depth  =  4.8  mm. 

Legend: 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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LED.  Assuming  the  LED  to  be  a  point  light  source,  this  aperture 
provided  an  absolute  restriction  of  approximately  13  degrees  on 
the  half-angle  divergence  of  the  cone  of  e-m  radiation  entering 
the  sample  chamber,  and  prevented  the  0-ring  from  being  directly 
illuminated. 

After  making  these  modifications,  we  repeated  the  measure- 
ments for  the  empty  chamber,  a  single  kernel,  and  five  kernels. 
Figures  4.7  and  4.8  show  the  results  of  these  studies.  Since  the 
total  amount  of  e-m  radiation  entering  the  chamber  was  reduced  by 
the  installed  aperture,  the  amplitude  of  the  pho toacous tic  re- 
sponse is  reduced  in  all  cases,  as  is  theoretically  expected. 
However,  the  separation  between  the  amplitude  of  the  empty  cham- 
ber response  and  the  response  obtained  with  a  wheat  sample  pre- 
sent was  improved  substantially,  except  for  at  frequencies  above 
approximately  200  Hz,  where  the  curves  are  found  to  converge. 
Also,  at  frequencies  above  approximately  15  Hz,  the  phase  of  the 
single  kernel  response  become  distinctly  different  from  the  empty 
chamber  response,  although  converging  once  again  at  frequencies 
approaching  1000  Hz.  There  are  also  some  significant  differences 
which  occur  in  both  the  amplitude  and  phase  angle  vs  modulation 
frequency  curves  for  a  single  kernel  and  five  kernels.  Because 
of  the  increased  separation  from  the  empty  chamber  response 
achieved,  the  modifications  were  considered  an  improvement,  and 
were  maintained  in  all  experiments  subsequently  performed  with 
the  cell. 

Given  the  appropriate  parameters  for  the  thermal  properties 
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Figure  4.7.   A  log-log  comparison  of  the  amplitude  of 

the  photoacous tic  response  vs  modulation  frequency 
obtained  with  a  single  kernel  and  five  kernels  of 
wheat  present  to  that  obtained  from  the  empty  cham- 
ber.  An  aperture  (diameter  =  3  mm)  was  installed 
to  restrict  spot  size.   Sample  chamber  depth  =  4.7  mm. 

Legend : 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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Figure  4.8.   A  comparison  of  the  phase  of  the  photo- 
acoustic  response  vs  log  modulation  frequency  ob- 
tained with  a  single  kernel  and  five  kernels  of 
wheat  present  to  that  obtained  from  the  empty  cham- 
ber.  An  aperture  (diameter  =  3  mm)  was  installed 
to  restrict  spot  size.   Sample  chamber  depth  =  4.7  mm. 

Legend: 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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of  helium,  aluminum,  and  wheat,  the  Rosencwaig-Gersho  theory 
should  permit  a  theoretical  determination  of  the  relative  magni- 
tudes of  the  amplitude  of  the  pho toacous t lc  response  for  these 
two  substances,  as  well  as  the  manner  with  which  the  amplitude 
will  change  with  variation  of  the  modulation  frequency.  Known 
values  for  the  density,  thermal  conductivity,  and  specific  heat 
for  the  materials  associated  with  these  studies  are  tabulated 
below  (21,  24).  In  addition,  a  plot  showing  the  thermal  diffu- 
sion length  for  these  materials  as  a  function  of  modulation 
frequency,  as  derived  from  these  parameters,  is  given  in  Figure 
4.9. 

TABLE   4.1.      Thermal    Parameters. 


Therm.    Conduct. 
J/m-sec-K 


0.1423 

237 

0.140 

0.117 

0.128 


Special  case  2(b)  of  the  Rosencwaig-Gersho  theory  should  be 
most  applicable  in  the  case  of  the  empty  chamber  response,  since 
the  thermal  diffusion  length  of  aluminum,  which  Is  at  a  maximum 
of  1.4  cm  at  1  Hz,  Is  less  than  the  thickness  of  the  chamber 
walls,  and  never  becomes  less  than  the  optical  absorption  length 
throughout    the    frequency    range    used.       Theory    thus    predicts     that 


Material 

Density 
Kg/m3 

Helium  (1  atm) 

0.180 

Aluminum 

2699 

Wheat  (hard) 

1300 

Wheat  (soft) 

1320 

Wheat  (avg) 

1310 

Specific 

Heat 

J/Kg-K 

5232 

899 

1548 

1397 

1473 
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Figure  4.9.   A  plot  of  the  thermal  diffusion  length 

as  a  function  of  modulation  frequency  for  various 
materials  of  interest  in  these  experiments. 

Legend: 

square  -  helium 
plus  -  aluminum 
diamond  -  wheat  (hard) 
triangle  -  wheat  (soft) 
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the  signal  amplitude  should  vary  as  to-1,  which  is  experimentally 
observed  at  frequencies  above  approximately  300  Hz.  At  lower 
frequencies,  however,  the  log  amplitude  vs  log  modulation  fre- 
quency curve  is  found  to  deviate  from  this  dependence,  with  the 
signal  amplitude  not  increasing  as  much  as  predicted  with  de- 
creasing frequency.  The  reason  for  this  behavior  is  most  likely 
that,  at  low  modulation  frequencies,  the  thermal  diffusion  length 
becomes  significantly  large  with  respect  to  the  spot  size  being 
illuminated,  making  a  one-dimensional  theory  no  longer  applica- 
ble, as  three-dimensional  diffusion  occurs  from  the  illuminated 
portion    of    the    chamber. 

In  the  case  of  a  wheat  kernel,  which  is  optically  opaque  as 
well,  the  thermal  diffusion  length,  which  is  approximately  0.37 
mm  at  a  modulation  frequency  of  1  Hz,  is  also  less  than  the 
sample  thickness  throughout  the  entire  frequency  range  used.  If 
It  is  assumed  that  the  optical  absorption  length  Is  less  than  the 
thermal  diffusion  length  of  the  sample,  then  In  this  Instance,  as 
was  the  case  with  aluminum,  special  case  2(b)  of  the  Rosencwaig- 
Gersho  theory  should  most  accurately  describe  the  experimental 
observations.  However,  when  the  ratio  of  the  amplitude  for  a 
wheat  kernel  to  the  amplitude  for  the  empty  chamber  is  computed 
making  this  assumption,  theoretically  a  signal  with  amplitude 
approximately  50  times  greater  should  be  observed  with  a  wheat 
kernel  sample  present,  throughout  the  entire  frequency  range. 
Such  a  large  separation  is  clearly  not  seen  In  the  experimental 
data,    where    the    greatest    separation    is    no    more    than   approximately 
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a  factor  of  two.  It  is,  however,  quite  possible  that  the  optical 
absorption  length  is  greater  than  the  thermal  diffusion  length, 
as  the  thermal  diffusion  length  for  wheat  is  fairly  small,  even 
for  low  modulation  frequencies.  If  this  is  correct,  then,  the 
amplitude    of    the    signal    would    be    dependent    upon    the    optical 

absorption   coefficient   of   wheat,    and   would    theoretically    vary   as 

-3/5 

to  '  .  The  fact  that  the  slope  of  the  log  amplitude  vs  log  modu- 
lation frequency  curve  experimentally  obtained  Is  more  negative 
with  both  one  and  five  wheat  samples  present  In  the  chamber 
substantiates  this  hypothesis.  Although  the  slope  of  this  curve 
deviates  from  the  expected  -3/2  slope,  except  for  at  modulation 
frequencies  approaching  1000  Hz,  where  the  slope  appears  to 
converge  to  this  value,  the  degree  to  which  it  deviates  from 
theory  is  similar  to  the  deviation  of  the  empty  chamber  curve 
from  Its  expected  slope  of  -1.  This  suggests  that  these  observed 
deviations  may  be  due  to  the  same  causes,  namely,  that  the  three- 
dimensional  nature  of  the  sample  becomes  significant  as  the 
thermal  diffusion  length  becomes  sufficiently  large.  Assuming 
that  the  Rosencwaig-Gersho  theory  is  valid  at  1000  Hz,  and  that 
the  optical  absorption  length  is  greater  than  the  thermal  diffu- 
sion length  for  wheat  at  this  frequency;  in  order  for  the  re- 
sponse of  a  wheat  kernel  sample  to  have  an  amplitude  equal  to  the 
amplitude  obtained  for  the  empty  chamber,  as  is  experimentally 
observed,  the  optical  absorption  length  for  a  wheat  kernel  would 
have  to  be  approximately  0.56  mm.  This  number  seems  to  be  rea- 
sonable,   and    implies    that    the   optical    absorption    length    is    great- 
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er  than  the  thermal  diffusion  length,  even  at  modulation  frequen- 
cies as  low  as  1  Hz,  where  the  thermal  diffusion  length  was 
calculated  to  be  approximately  0.36  mm. 

Another  result  obtained  from  this  experiment  was  that,  at 
the  lower  modulation  frequencies  (less  than  and  approximately 
equal  to  40  Hz),  the  amplitude  of  the  pho  toacous  t  ic  signal  ob- 
tained for  a  single  kernel  was  greater  than  that  obtained  for 
five  kernels.  However,  at  higher  frequencies,  there  was  essen- 
tially no  difference  between  the  signals.  This  variation  is  most 
likely  due  to  a  decrease  in  the  volume  of  gas  occupying  the  cell, 
as  more  volume  was  occupied  by  the  sample  with  five  kernels 
present,  causing  the  thermal  diffusion  length  of  helium  to  become 
comparable  to  the  effective  chamber  length  at  a  higher  modulation 
frequency  than  in  the  case  of  a  single  kernel.  For  this  experi- 
ment (with  the  aperture  installed),  the  chamber  was  maintained  at 
a  depth  of  4.7  mm.  If  an  average  kernel  thickness  of  2.4  mm  is 
subtracted  from  this  depth  with  5  kernels  present,  the  resultant 
effective  length  of  the  gas  column  becomes  approximately  2.3  mm. 
The  modulation  frequency  at  which  the  thermal  diffusion  length 
for  helium  is  2.3  mm  is  approximately  50  Hz,  which  closely  corre- 
sponds to  the  40  Hz  modulation  frequency  at  which  the  variation 
between  the  single  kernel  and  five  kernel  response  becomes  appar- 
ent. The  fact  that  the  entire  chamber  floor  was  not  completely 
covered  when  five  kernels  were  present,  due  to  gaps  between  the 
kernels,  is  possibly  responsible  for  the  fact  that  the  above 
approximation  gives  a  theoretical  frequency  at  which  separation 
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between  the  curves  should  begin  slightly  higher  than  that  actual- 
ly observed,  as  a  larger  gas  volume  than  that  allowed  for  In  the 
above   approximation   was    present. 

C.    Chamber   Volume   Studies 

The  effects  of  the  size  of  the  sample  chamber  In  relation  to 
the  thermal  diffusion  length  of  the  gas  used,  In  this  case  he- 
lium, have  been  treated  in  the  discussion  above.  However,  all  of 
the  experiments  previously  discussed  were  performed  at  essential- 
ly the  same  chamber  depth.  The  purpose  of  the  following  experi- 
ment was  to  study  specifically  how  an  actual  variation  in  the 
depth  of  the  sample  chamber  will  affect  the  pho toacous tic  re- 
sponse, throughout  the  entire  modulation  frequency  range  of  1- 
1000   Hz. 

For  this  experiment,  the  photoacoustic  cell  was  adjusted  so 
that  its  chamber  depth  was  1.66  cm,  which  was,  by  design,  the 
maximum  chamber  depth  possible.  At  this  chamber  depth,  runs  were 
made  for  the  empty  chamber,  a  single  kernel,  five  kernels,  and 
the  black  cloth  sample.  Figures  4.10  and  4.11  show  the  amplitude 
and  phase  variation  with  modulation  frequency  of  the  black  cloth 
photoacoustic  response  at  this  new  chamber  depth,  as  compared  to 
the  response  obtained  at  a  chamber  depth  of  0.47  cm.  It  is  clear- 
ly seen  that  the  amplitude  of  the  response  Is  decreased  with  the 
increased  chamber  depth,  as  is  theoretically  predicted.  However, 
since  the  volume  of  the  chamber  was  increased  by  more  than  a 
factor  of  three,  theory  predicts  that  the  amplitude  of  the  signal 
at    this   new  chamber   depth   should  be  approximately  one-third  as 
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Figure  4.10.   A  log-log  comparison  of  the  amplitude  of 
the  photoacous tic  response  vs  modulation  frequency 
obtained  for  the  black  cloth  sample  at  two  different 
chamber  depths. 

Legend: 

square  -  chamber  depth  =  4.7  mm 
plus  -  chamber  depth  =  17  mm 
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Figure  4.11.   A  comparison  of  the  phase  of  the  photo- 
acoustic  response  vs  log  modulation  frequency  ob- 
tained for  the  black  cloth  sample  at  two  different 
chamber  depths. 

Legend: 

square  -  chamber  depth  =  4.7  mm 
plus  -  chamber  depth  ■  17  mm 
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great  as  the  amplitude  at  the  more  shallow  depth.  Experimental- 
lyf  we  see  an  amplitude  approximately  0.7  times  as  great.  A 
possible  reason  for  this  deviation  from  theory  is  that,  at  the 
larger  chamber  depth,  a  greater  amount  of  the  black  cloth  surface 
is  being  illuminated  due  to  the  increased  beam  spreading  of  the 
incident  e-m  radiation.  This  would  cause  an  increase  in  the 
volume  of  gas  actually  being  heated  in  the  case  of  the  larger 
chamber  depth,  and  may  result  in  complicating  a  strictly  theore- 
tical comparison  of  the  two  cases.  Another  factor  which  may  be 
responsible  for  the  deviation  of  our  experimental  results  from 
theoretical  predictions  is  the  assumption  of  one-dimensional  heat 
propagation  in  the  Rosencwaig-Gersho  theory,  whereas,  in  reality, 
propagation  of  heat  parallel  to  the  surface  occurs  within  the 
sample  as  well.  Since  as  the  chamber  depth  is  made  larger,  the 
amount  of  material  illuminated  increases,  it  is  reasonable  to 
expect  that  the  effect  of  three-dimensional  propagation  of  heat 
at  larger  chamber  depths  is  less  significant  than  at  shallow 
depths.  Thus,  in  our  experiment,  it  is  possible  that  at  the 
larger  chamber  depth  a  one-dimensional  approximation  is  valid, 
whereas  at  the  shallower  depth,  the  three-dimensional  nature  of 
the  heat  propagation  process  becomes  significant,  leading  to  a 
higher    percentage    of   heat   being    dissipated    within    the    sample. 

Since  the  thermal  diffusion  length  for  helium  at  1  Hz  is 
approximately  equal  to  the  length  of  the  sample  chamber  at  this 
new  chamber  depth,  as  the  modulation  frequency  approaches  1  Hz, 
the    point   where    the    thermal    diffusion    length    of    the    gas    becomes 
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equal  to  the  length  of  the  gas  column  Is  approached.  It  Is  seen 
that  at  these  low  modulation  frequencies,  the  slope  of  the  log 
amplitude  vs  log  frequency  curve  becomes  less  negative  for  the 
1.66  cm  chamber  depth,  although  more  gradually  than  at  the  0.47 
cm  depth.  This  behavior  Is  reasonable  since  for  the  deeper 
chamber  depth,  the  thermal  diffusion  length  of  helium  only  be- 
comes comparable  to  the  gas  column  length;  whereas  for  the  shal- 
lower chamber  depth,  the  thermal  diffusion  length  of  helium 
actually  becomes  significantly  greater  than  the  gas  column  depth, 
at    low  modulation    frequencies. 

Figures  4.12  and  4.13  show  the  amplitude  and  phase  of  the 
pho toacous t ic  response  at  a  chamber  depth  of  1.66  cm  for  the 
empty  chamber,  a  single  kernel,  and  5  kernels.  The  change  In  the 
slope  of  the  amplitude  curves  at  low  modulation  frequencies  is 
much  more  gradual  than  that  seen  In  Figure  4.7  for  the  smaller 
chamber  depth,  as  is  expected.  The  relative  amplitude  of  the 
pho  toacous tic  response  for  five  kernels  as  compared  to  a  single 
kernel  at  this  depth  seems  to  be  shifted  up  somewhat  with  respect 
to  the  results  obtained  at  the  lower  chamber  depth  throughout  the 
entire  frequency  range.  Here  we  find  that  at  the  lower  modula- 
tion frequencies,  the  single  kernel  and  five  kernel  amplitudes 
are  nearly  identical,  whereas,  at  the  higher  modulation  frequen- 
cies, the  amplitude  for  five  kernels  becomes  greater  than  that 
obtained  for  a  single  kernel.  This  is  possibly  due  to  an  in- 
creased contribution  from  those  kernels  surrounding  the  central 
kernel,    as    these  kernels   receive   a   greater   amount   of    Incident 
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Figure  4.12   A  log-log  plot  of  the  amplitude  of  the 

photoacous tic  response  vs  modulation  frequency  for 
the  empty  chamber,  a  single  kernel,  and  five  ker- 
nels.  Sample  chamber  depth  ■  17  mm. 

Legend : 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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Figure  4.13.   A  plot  of  the  phase  of  the  photoacoustlc 
response  vs  log  modulation  frequency  for  the  empty 
chamber,  a  single  kernel,  and  five  kernels.   Sam- 
ple chamber  depth  ■  17  mm. 

Legend : 

square  -  no  sample  present  (empty  chamber) 
plus  -  one  kernel  present 
diamond  -  five  kernels  present 
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radiation  due  to  the  increased  spot  size  at  this  chamber  depth. 
This  reasoning  parallels  that  used  above,  in  the  case  of  the 
black  cloth  signal  comparisons,  to  help  explain  the  discrepancy 
between  the  theoretical  and  experimental  drop  in  the  signal 
amplitude.  From  these  results,  it  seems  that  increasing  the 
surface  area,  while  keeping  the  total  radiation  flux  entering  the 
sample  chamber  constant,  results  in  an  increase  in  the  amplitude 
of  the  response.  Similar  reasoning  can  be  used  to  explain  the 
fact  that  the  separation  between  the  empty  chamber  and  single 
kernel  amplitudes  seems  to  decrease  at  this  larger  chamber  depth. 
Since  the  empty  chamber  has  a  much  larger  surface  area  than  a 
single  kernel,  an  increase  in  the  spot  size,  assuming  that  it 
becomes  greater  than  the  size  of  the  kernel,  will  result  in  a 
larger  increase  in  the  amount  of  the  empty  chamber  illuminated  as 
compared  to  the  increased  amount  of  wheat  illuminated.  The  phase 
relationship  of  the  single  kernel  pho toacous tic  response  with 
modulation  frequency  appears  to  become  less  distinctly  separate 
from  the  empty  chamber  response  as  the  chamber  depth  is  in- 
creased. This  may  be  due  to  the  increased  contribution  of  illum- 
inated portions  of  the  empty  chamber  when  the  single  kernel  run 
was  performed  because  of  the  Increased  spot  size.  The  phase 
relationship  with  five  kernels  present  appears  to  be  essentially 
the  same  for  both  chamber  depths,  which  is  reasonable  since  as 
the  spot  size  increases,  kernels  surrounding  the  central  kernel 
are  illuminated  instead  of  the  chamber  itself,  causing  the  type 
of  material  being  illuminated  to  be  essentially  unchanged. 
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4.2      Reproducibility    Studies 

A  series  of  runs  were  made  to  find  out  the  extent  to  which 
reproducible  data  could  be  obtained.  Although  these  studies  were 
done  using  an  aperture  of  somewhat  larger  diameter  than  that 
previously  described,  they  still  should  give  a  reasonably  valid 
representation  of  the  variability  in  our  results  due  to  uncer- 
tainties in  the  signal  measurement  by  the  electronics  systems, 
along  with  variations  in  the  percentage  of  helium  and  presence  of 
contaminants    in    the    chamber    between    runs. 

Figures  4.14  and  4.15  show  the  amplitude  and  phase  results 
of  five  runs  for  which  the  empty  chamber  response  was  monitored 
as  a  function  of  modulation  frequency,  throughout  the  entire  1- 
1000  Hz  frequency  range.  The  runs  1-5  were  sequentially  per- 
formed, without  reflushing  of  the  sample  chamber  between  runs.  A 
single  reflushing  of  the  chamber  was  performed  only  when  the 
internal  component  cards  of  the  lock-in  amplifier  were  changed, 
and  the  higher-frequency  range  portion  of  the  runs  were  begun. 
These  results  should  represent  the  limit  to  which  variations 
between  samples  can  be  detected  with  the  apparatus  used  in  these 
experiments.  In  the  case  of  the  amplitude,  it  is  seen  that  the 
scatter  between  runs  is  not  large.  The  percentage  of  uncertainty 
represented  by  the  standard  deviation  of  the  values  obtained  at  a 
given  frequency,  averaged  over  all  frequencies  at  which  measure- 
ments were  made,  gives  an  experimental  uncertainty  of  1.70%  for 
the    amplitude    and    1.78%    for    the    phase    measurement. 

Some   of    the   most  probable   sources   of   uncertainty   in   these 
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Figure  4.14.   A  log-log  plot  of  the  amplitude  of  the 
photoacoustlc  response  vs  modulation  frequency 
for  the  empty  chamber.   The  results  of  five  separ- 
ate runs  which  were  sequentially  performed  without 
ref lushing  the  sample  chamber  are  plotted. 

Legend: 

square  -  run  #1 
plus  -  run  #2 
diamond  -  run  #3 
triangle  -  run  #4 
cross  -  run  #5 
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Figure  4.15.   A  plot  of  the  phase  of  the  pho toacous tic 
response  vs  log  modulation  frequency  for  the  empty 
chamber.   The  results  of  five  separate  runs  which 
were  sequentially  performed  without  reflushing  the 
sample  chamber  are  plotted. 

Legend: 

square  -  run  #1 
plus  -  run  #2 
diamond  -  run  If 3 
triangle  -  run  M 
cross  -  run  #5 
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measurements  include  variations  in  the  intensity  of  the  e-m 
radiation  intensity  and  a  decrease  with  time  in  the  percentage  of 
helium  actually  present  in  the  chamber  due  to  leakage,  as  well  as 
possible  fluctuations  in  ambient  temperature,  variations  in  out- 
side noise  levels,  and  drift  or  uncertainty  in  the  microphone's 
signal  detection  capability.  Studies  in  which  the  LED  intensity, 
as  detected  by  the  PIN  diode,  was  monitored  over  a  period  of  time 
show  that  the  intensity  of  the  incident  radiation  does  decrease 
very  slightly  over  a  period  of  continuous  usage,  with  a  decrease 
of  no  more  than  one  percent  over  several  hours.  Looking  closely 
at  the  amplitude  data,  it  is  seen  that  the  slope  of  the  log 
amplitude  vs  log  frequency  curve  becomes  slightly  more  negative 
with  successive  runs.  It  is  likely  that  this  change  in  the 
frequency  dependence  of  the  amplitude  is  due  to  a  gradual  diffu- 
sion of  air  into,  and  helium  out  of,  the  sample  chamber.  This 
hypothesis  is  backed  by  the  fact  that  the  empty  chamber  amplitude 
without  helium  flushing  is  greater  than  with  helium  flushing  at 
low  modulation  frequencies,  and  less  than  the  helium  flushed 
chamber  signal  at  higher  modulation  frequencies.  This  implies 
that  it  is  advantageous  to  perform  the  run  soon  after  the  chamber 
is  flushed.  Another  source  of  uncertainty,  which  ultimately  dic- 
tates the  minimum  accuracy  of  the  acquired  data,  is  the  limited 
resolution  of  the  A-D  converter  used  in  the  data  acquisition 
process.  Although  the  consistency  of  the  analog  to  digital 
conversion  process  was  made  was  very  good,  as  tested  by  inputting 
constant  voltages   and  monitoring    the    resultant    digital    values, 
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the  accuracy  to  which  the  converter  could  detect  changes  in 
voltage  was  limited  to  1/255  of  its  0.5  volt  working  range,  or 
approximately  1.5%  at  half  scale  output  from  the  lock-in  ampli- 
fier. 

Figures  4.16  and  4.17  show  the  results  of  an  experiment 
similar  to  that  described  above,  in  which  five  runs  were  made 
monitoring  the  photoacous tic  response  of  a  single  wheat  kernel  as 
a  function  of  modulation  frequency.  For  each  run  in  this  experi- 
ment, however,  the  sample  was  removed,  and  then  centered  once 
again  in  the  sample  chamber,  and  the  chamber  was  reflushed. 
Despite  the  potential  increased  uncertainties  associated  with  the 
repositioning  of  the  kernel,  presence  of  unwanted  contaminants, 
and  the  percentage  of  helium  successfully  introduced  into  the 
chamber  for  each  run,  the  scatter  in  this  data  is  even  slightly 
less  than  the  scatter  observed  in  the  empty  chamber  runs  de- 
scribed above,  with  an  average  percent  standard  deviation  from 
average  of  1.25%  for  the  amplitude  and  0.96%  for  the  phase.  The 
decrease  in  scatter  in  this  case  is  probably  due  to  the  fact  that 
the  chamber  was  reflushed  prior  to  each  run,  and  an  additive 
reduction  in  the  percentage  of  helium  in  the  chamber  does  not 
occur  with  successive  runs.  The  scatter  that  is  observed  is 
likely  caused  by  other  factors  cited  in  the  above  discussion  for 
the  empty  chamber  experiment,  which  are  relevant  to  this  case. 
The  experimental  procedures  used  in  this  particular  study  closely 
corresponded  to  those  used  in  the  single  kernel  comparisons  later 
discussed.       These    results    therefore    demonstrate    that    a    satisfac- 
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Figure  4.16.   A  log-log  plot  of  the  amplitude  of  the 
photoacous tic  response  vs  modulation  frequency 
for  a  single  wheat  kernel  (variety:  Scout-66). 
For  each  of  the  five  runs  plotted,  the  kernel  was 
removed  and  placed  in  the  sample  chamber  again, 
and  the  chamber  was  reflushed  with  helium. 

Legend: 

square  -  run  #1 
plus  -  run  #2 
diamond  -  run  #3 
triangle  -  run  #4 
cross  -  run  #5 
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Figure  4.17.   A  plot  of  the  phase  of  the  photoacous tic 
response  vs  log  modulation  frequency  for  a  single 
wheat  kernel  (variety:  Scout-66).   For  each  of  the 
five  runs  plotted,  the  kernel  was  removed  and  placed 
In  the  sample  chamber  again,  and  the  chamber  was 
reflushed  with  helium. 

Legend: 

square  -  run  #1 
plus  -  run  #2 
diamond  -  run  #3 
triangle  -  run  #4 
cross  -  run  #5 
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tory    degree    of    reproducibility    is    obtainable    with    the    apparatus 
and    procedures    used    in    these    studies. 

4.3      Kernel    Moisture    Evaporation    Studies 

In  earlier  photoacoustic  studies  of  wheat  kernels  (25),  it 
has  been  seen  that  moisture  content  has  a  very  pronounced  effect 
on  the  amplitude  of  the  photoacoustic  response.  Due  to  the  high 
specific  heat  of  water,  its  presence  in  a  kernel  results  in  an 
reduction  In  the  amplitude  of  the  photoacoustic  response,  as, 
theoretically,  the  amplitude  varies  inversely  as  the  specific 
heat  for  an  opaque,  thermally  thick  solid  with  a  thermal  diffu- 
sion length  less  than  Its  optical  absorption  length.  Further- 
more, moisture  content  is  a  kernel  property  which  has  the  poten- 
tial to  vary  substantially  between  samples.  Therefore,  In  order 
to  make  valid  comparisons  of  the  Inherent  photoacoustic  signals 
from  different  kernels,  it  is  first  necessary  to  bring  those 
kernels  to  be  compared  to  the  same  state,  with  respect  to  moist- 
ure content.  The  purpose  of  this  study  is  to  see  how  the  photoa- 
coustic response  varies  as  moisture  is  removed  from  a  wheat 
kernel    by    drying. 

These  studies  were  performed  relatively  early  In  the  course 
of  our  experimentation,  prior  to  the  earlier  described  O-ring  and 
aperture  modifications  made  to  the  photoacoustic  cell,  and  before 
we  had  the  capability  to  extend  our  studies  to  modulation  fre- 
quencies below  10  Hz  (before  we  acquired  the  1-100  Hz  component 
boards    for    the    lock-in    amplifier).       Thus,     the    photoacoustic    data 


120 


is  not  directly  comparable  to  data  later  obtained,  and  the  re- 
sults may  be  somewhat  inconclusive.  Although  a  number  of  differ- 
ent varieties  of  kernels  were  examined  in  this  study,  they  all 
demonstrated  a  similar  variation  in  their  photoacous tic  response 
with  Increased  drying.  Thus,  only  the  results  for  two  represent- 
ative kernels  with  similar  initial  kernel  mass  are  considered 
here,  these  varieties  being  Arkan  (hard  red  winter)  and  Pike 
(soft  red  winter). 

Since  one  of  the  purposes  of  this  study  was  to  determine  the 
length  of  drying  time  necessary  for  preparation  of  wheat  samples 
for  our  future  experiments,  no  special  sample  preparation  was 
done  prior  to  the  study,  and  the  kernels  were  studied  with  an 
initial  moisture  content  as  received.  Before  the  drying  process 
was  begun,  each  kernel  was  weighed  and  its  photoacous tic  response 
was  obtained  throughout  the  10-1000  Hz  frequency  range.  The 
samples  were  then  dried  in  a  conventional  air  oven  at  a  tempera- 
ture of  103  degrees  Centigrade  (+/-2  degrees)  for  five  hours,  re- 
moved, weighed,  and  the  pho toacous t ic  response  was  again  mea- 
sured. This  process  was  continued,  with  increased  drying  times 
being  allotted  between  samplings  as  the  mass  and  signal  varia- 
tions became  more  gradual,  until  a  total  drying  time  of  95  hours 
was  reached.  Figures  4.18-4.21  display  the  results  obtained  for 
the  Pike  and  Arkan  kernels  analyzed,  throughout  the  entire  fre- 
quency range  scanned,  for  both  the  amplitude  and  phase  of  the 
response.  Figure  4.22  shows  the  variation  in  mass  as  a  function 
of  drying  time  for  the  both  kernels.   Assuming  that  mass  loss 
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represents  moisture  evaporation  from  the  kernel,  this  plot  should 
give  a  relative  indication  of  the  variation  in  total  kernel 
moisture  content.  Figures  4.23  and  4.24  more  clearly  demonstrate 
the  change  in  photoacous tic  signal  amplitude  with  drying  time  for 
modulation    frequencies    of    10    and    100   Hz    respectively. 

These    results    show    that   as    the    drying    time    is    increased,    and 
the    total    kernel    moisture    content    is    decreased,     as    indicated    by 
the  decreased  kernel   mass,    the   amplitude   of    the   photoacous tic   re- 
sponse   increases.      The    change    in    the   amplitude    is    rapid    when    the 
drying   process    begins,    but   becomes    more    gradual    as    the    drying 
time   increases,   appearing   to  eventually   level   off,   as   does    the 
kernel    moisture    content.       The    fact    that    the    signal    amplitude 
seems    to   vary   inversely  with  moisture  content   is    in  compliance 
with    theory,   as    the   Rosencwaig-Gersho    theory   predicts    that   ampli- 
tude   should    vary   as    the    inverse   of   specific   heat   for   a   solid   with 
the    thermal    and   optical    properties    of   wheat,    and    the   specific 
heat  of  water   is   much  higher   than  any  other  substances  of  which 
the   kernel    is    composed.      Upon   examining    the    amplitude    of    the 
photoacoustic    response    shown    in    Figures    4.18    and    4.20,     it    is    seen 
that   the  amplitude  does   not   increase  uniformly  over    the   entire 
frequency   range.      Due    to   the   fact   that   the  presence  of  moisture 
in    a    kernel    has    such    a    large    effect    on    the    photoacoustic    re- 
sponse,   and    variation  of    the   modulation   frequency   enables    the 
profiling  of  a   sample  as   a   function   of   depth,    it    is    reasonable    to 
expect   that   this   behavior  reflects   a  non-uniform  extraction  of 
moisture    from    a    given    depth    within    the    kernel.       If    the    change    in 
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Figure  4.18.   A  log-log  plot  showing  the  variation  in 
the  amplitude  of  the  photoacoustic  response  vs 
modulation  frequency  for  a  single  kernel  sample  of 
Pike  (soft  red  winter),  as  the  kernel  was  dried 
over  a  period  of  75  hours.   Oven  temperature  = 
103  deg.  C. 

Legend: 

square  -  before  drying  process  begun 

plus  -  after  5  hours  drying  time 

diamond  -  after  10  hours  drying  time 

cross  -  after  35  hours  drying  time 

inverted  triangle  -  after  75  hours  drying  time 
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Figure  4.19.   A  plot  showing  the  variation  In  the  phase 
of  the  photoacoustlc  response  vs  log  modulation 
frequency  for  a  single  kernel  sample  of  Pike  (soft 
red  winter),  as  the  kernel  was  dried  over  a  period 
of  75  hours.   Oven  temperature  =  103  deg.  C. 

Legend: 

square  -  before  drying  process  begun 

plus  -  after  5  hours  drying  time 

diamond  -  after  10  hours  drying  time 

cross  -  after  35  hours  drying  time 

inverted  triangle  -  after  75  hours  drying  time 
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Figure  4.20.   A  log-log  plot  showing  the  variation  in 
the  amplitude  of  the  photoacous tic  response  vs 
modulation  frequency  for  a  single  kernel  sample  of 
Arkan  (hard  red  winter),  as  the  kernel  was  dried 
over  a  period  of  75  hours.   Oven  temperature  ■ 
103  deg.  C. 

Legend: 

square  -  before  drying  process  begun 

plus  -  after  5  hours  drying  time 

diamond  -  after  10  hours  drying  time 

cross  -  after  35  hours  drying  time 

inverted  triangle  -  after  75  hours  drying  time 
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Figure  4.21.   A  plot  showing  the  variation  in  the  phase 
of  the  photoacoustic  response  vs  log  modulation 
frequency  for  a  single  kernel  sample  of  Arkan  (hard 
red  winter),  as  the  kernel  was  dried  over  a  period 
of  75  hours.   Oven  temperature  =  103  deg.  C. 

Legend: 

square  -  before  drying  process  begun 

plus  -  after  5  hours  drying  time 

diamond  -  after  10  hours  drying  time 

cross  -  after  35  hours  drying  time 

inverted  triangle  -  after  75  hours  drying  time 
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Figure  4.22.   A  plot  showing  the  variation  in  the 

total  kernel  mass  with  increased  drying  time  (at 
103  deg.  C.)  for  the  varieties  Pike  and  Arkan. 

Legend: 

square  -  Pike  (soft  red  winter) 
diamond  -  Arkan  (hard  red  winter) 
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Figure  4.23.   A  plot  showing  the  variation  In  the  amp- 
litude of  the  photoacoustic  response  with  Increased 
drying  time  (at  103  deg.  C),  at  a  modulation  fre- 
quency of  10  Hz,  for  the  varieties  Pike  and  Arkan. 

Legend: 

square  -  Pike  (soft  red  winter) 
diamond  -  Arkan  (hard  red  winter) 
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Figure  4.24.   A  plot  showing  the  variation  in  the  amp- 
litude of  the  photoacoustic  response  with  increased 
drying  time  (at  103  deg.  C),  at  a  modulation  fre- 
quency of  100  Hz,  for  the  varieties  Pike  and  Arkan. 

Legend: 

square  -  Pike  (soft  red  winter) 
diamond  -  Arkan  (hard  red  winter) 
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the  shape  of  these  curves  are  interpreted  in  this  manner,  it  is 
seen  that  as  the  drying  is  begun,  a  larger  change  in  kernel 
moisture  occurs  more  deeply  within  the  kernel  (lower  modulation 
frequencies)  than  at  the  kernel  surface  (higher  modulation  fre- 
quencies).  After  a  period  of  drying  time,  however,  the  amount  of 
moisture  at  the  kernel  surface  greatly  decreases,  and  afterward 
the  variation  seems  to  become  relatively  uniform  throughout  the 
entire  range  of  depths  from  which  information  was  obtained.   This 
seems  to  parallel  the  manner  by  which  moisture  is  presently 
understood  to  be  driven  from  a  wheat  kernel  upon  drying.   If  a 
kernel  has  not  been  recently  damped  (bran  drier  than  endosperm), 
moisture  will  first  be  drawn  from  the  interior  of  the  grain  until 
a  moisture  distribution  in  which  the  bran  is  damper  than  the 
endosperm  is  established  (26).   It  appears  from  our  experiment 
that  after  much  of  the  moisture  has  been  drawn  from  the  interior 
of  the  kernel,  a  rather  extreme  reduction  in  the  kernel's  surface 
moisture  occurs.   Although  the  phase  of  the  signal  did  vary  with 
drying  time,  as  seen  in  Figures  4.19  and  4.21,  the  behavior  was  a 
great  deal  more  complex  than  the  amplitude  variation.   No  con- 
sistent trend  of  either  decrease  or  increase  in  phase  angle  with 
drying  time  was  observed,  and  the  manner  with  which  it  did  change 
was  not  the  same  for  each  kernel  tested.   It  is  possible  that  the 
phase  lag  Is  more  sensitive  to  the  distribution  of  moisture 
throughout  the  portion  of  the  kernel  In  which  the  thermal  wave 
must  propagate  to  reach  the  kernel  surface.   Thus,  even  though 
the  phase  of  the  signal  seems  to  be  dependent  on  the  kernel 
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moisture  content,  the  exact  nature  of  its  dependency  is  not 
clear. 

The  temperature,  103  degrees  Centigrade,  used  in  this  study 
was  chosen  since  a  USDA  standard  wheat  drying  technique  specifies 
heating  at  this  temperature  for  fifty  hours.  Since  our  results 
indicate  that  very  little  change  in  the  moisture  content  occurs 
after  fifty  hours  drying  time,  we  decided  to  use  these  standard 
temperature  and  drying  time  parameters  in  our  initial  preparation 
of  samples  tested  in  subsequent  studies. 

4.4   Kernel  Variety  Studies 

The  objectives  of  these  experiments  were  to  examine  the 
photoacoustic  response  of  a  number  of  different  wheat  varieties, 
and  attempt  to  pinpoint  differences  in  the  photoacoustic  response 
characteristic  of  a  given  variety.  Twelve  different  varieties  of 
wheat  were  chosen  to  be  used  in  this  experiment,  eight  hard  red 
winter  and  four  soft  red  winter  varieties.  The  specific  varieties 
studied  are  listed  below. 

TABLE  4.2.   Listing  of  Varieties  Analyzed. 

Hard  Red  Winter  Soft  Red  Winter 

Arkan  Arthur 

Bounty-310  Caldwell 

Hawk  Coker-916 

Mustang  plke 

Newton 

Scout-66 

Tam-105 

Vona 
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To  verify  our  observations,  runs  were  made  for  five  differ- 
ent kernels  of  each  variety.  In  order  to  bring  all  kernels 
examined  to  the  same  state  with  respect  to  moisture  content, 
moisture  was  removed  from  each  of  the  kernels  prior  to  the  study 
by  the  procedure  described  in  the  previous  section,  with  all 
kernels  studied  being  dried  at  the  same  time.  After  the  drying 
process,  the  kernels  were  stored  in  a  desiccator  up  until  and 
after  they  were  individually  examined.  The  analysis  was  per- 
formed for  all  kernels  first  at  the  low  frequency  range  (1-50 
Hz),  after  which  the  internal  component  boards  of  the  lock-in 
amplifier  were  exchanged  and  the  high  frequency  range  (10-1000 
Hz)  analysis  was  done,  completing  the  experiment.  As  in  previous 
studies,  the  chamber  was  flushed  for  five  minutes  with  helium 
before  each  run.  The  chamber  depth  used  for  this  experiment  was 
0.47    cm. 

Although  variations  were  observed  in  the  photoacous tic  re- 
sponse for  the  different  kernels  tested,  even  kernels  of  the  same 
variety  were  found  to  vary  somewhat  from  kernel  to  kernel.  Fig- 
ures 4.25-4.28  show  the  extent  of  the  variation  in  the  amplitude 
and  phase  of  the  pho toacous t ic  signal  for  the  varieties  that 
seemed  to  exhibit  the  most  and  least  variance  between  kernels. 
Figures  4.29-4.34  show  the  pho toacous tic  responses  for  the  twelve 
varieties  tested,  as  averaged  over  the  five  kernels  tested  within 
each  variety.  The  average  of  the  amplitude  and  phase  angle  vs 
modulation  frequency  curves  for  a  given  variety  were  seen  to 
differ    somewhat    between   varieties.      However,    the   kernel    to  kernel 
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Figure  4.25.   A  log-log  plot  of  the  amplitude  of  the 
photoacoustic  response  vs  modulation  frequency 
for  five  different  kernels  (variety:  Hawk),  each 
analyzed  separately.   This  plot  represents  the 
least  amount  of  kernel  to  kernel  scatter  In  the 
signal  amplitude  observed  within  a  given  variety 
for  the  twelve  varieties  tested. 

Legend : 

square  -  kernel  #1 
plus  -  kernel  #2 
diamond  -  kernel  #3 
triangle  -  kernel  #4 
cross  -  kernel  #5 
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Figure  4.26.   A  log-log  plot  of  the  amplitude  of  the 
photoacoustic  response  vs  modulation  frequency 
for  five  different  kernels  (variety:  Vona),  each 
analyzed  separately.   This  plot  represents  the 
most  kernel  to  kernel  scatter  in  signal  amplitude 
observed  within  a  given  variety  for  the  twelve 
varieties  tested. 

Legend: 

square  -  kernel  #1 
plus  -  kernel  #2 
diamond  -  kernel  #3 
triangle  -  kernel  #4 
cross  -  kernel  #5 
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Figure  4.27.   A  plot  of  the  phase  of  the  photoacous tic 
response  vs  log  modulation  frequency  for  five  dif- 
ferent kernels  (variety:  Mustang),  each  analyzed 
separately.   This  plot  represents  the  least  amount 
of  kernel  to  kernel  scatter  in  the  phase  angle  ob- 
served within  a  given  variety  for  the  twelve  varie- 
ties tested. 

Legend: 

square  -  kernel  #1 
plus  -  kernel  #2 
diamond  -  kernel  $3 
triangle  -  kernel  #4 
cross  -  kernel  #5 
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?igure  4.28.   A  plot  of  the  phase  of  the  photoacous tic 
response  vs  log  modulation  frequency  for  five  dif- 
ferent kernels  (variety:  Scout-66),  each  analyzed 
separately.   This  plot  represents  the  most  kernel 
to  kernel  scatter  in  the  phase  angle  observed  within 
a  given  variety  for  the  twelve  varieties  tested. 

Legend: 

square  -  kernel  #1 
plus  -  kernel  #2 
diamond  -  kernel  $3 
triangle  -  kernel  #4 
cross  -  kernel  #5 
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Figure  4.29.   A  log-log  plot  showing  the  average  amp- 
litude of  the  photoacoustlc  response  vs  modulation 
frequency  for  each  of  the  four  soft  red  winter 
wheat  varieties  tested.   The  photoacoustlc  responses 
from  a  total  of  five  separate  kernels  were  used  to 
obtain  the  average  results  for  each  variety. 

Legend: 

square  -  Arthur  (AT) 
plus  -  Caldwell  (Ctf) 
diamond  -  Coker-916  (CK) 
triangle  -  Pike  (PK) 
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Figure  4.30.   A  log-log  plot  showing  the  average  amp- 
litude of  the  photoacoustlc  response  vs  modulation 
frequency  for  four  of  the  eight  hard  red  winter 
wheat  varieties  tested.   The  photoacoustlc  responses 
from  a  total  of  five  separate  kernels  were  used  to 
obtain  the  average  results  for  each  variety. 

Legend: 

square  -  Arkan  (AK) 
plus  -  Bounty-310  (BY) 
diamond  -  Hawk  (HK) 
triangle  -  Mustang  (MG) 
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Figure  4.31.   A  log-log  plot  showing  the  average  amp- 
litude of  the  photoacoustic  response  vs  modulation 
frequency  for  the  remaining  four  hard  red  winter 
wheat  varieties  tested.   The  photoacoustic  responses 
from  a  total  of  five  separate  kernels  were  used  to 
obtain  the  average  results  for  each  variety. 

Legend : 

square  -  Newton  (NW) 
plus  -  Seout-66  (SC) 
diamond  -  Tam-105  (TM) 
triangle  -  Vona  (VN) 
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Figure  4.32.   A  plot  showing  the  average  phase  of  the 

photoacous tic  response  vs  log  modulation  frequency 
for  each  of  the  four  soft  red  winter  wheat  varie- 
ties tested.   The  photoacous tic  responses  from  a 
total  of  five  separate  kernels  were  used  to  obtain 
the  average  results  for  each  variety. 

Legend : 

square  -  Arthur  (AT) 
plus  -  Caldwell  (CW) 
diamond  -  Coker-916  (CK) 
triangle  -  Pike  (PK) 
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Figure  A. 33.   A  plot  showing  the  average  phase  of  the 

photoacoustic  response  vs  log  modulation  frequency 
for  four  of  the  eight  hard  red  winter  wheat  varie- 
ties tested.   The  photoacoustic  responses  from  a 
total  of  five  separate  kernels  were  used  to  obtain 
the  average  results  for  each  variety. 

Legend: 

square  -  Arkan  (AK) 
plus  -  Bounty-310  (BY) 
diamond  -  Hawk  (HK) 
triangle  -  Mustang  (MG) 
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Figure  4.34.   A  plot  showing  the  average  phase  of  the 

photoacous tic  response  vs  log  modulation  frequency 
for  the  remaining  four  hard  red  winter  wheat  varie- 
ties tested.   The  photoacous tic  responses  from  a 
total  of  five  separate  kernels  were  used  to  obtain 
the  average  results  for  each  variety. 

Legend: 

square  -  Newton  (NW) 
plus  -  Scout-66  (SC) 
diamond  -  Tam-105  (TM) 
triangle  -  Vona  (VN) 
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variation  within  a  given  variety  was  generally  at  least  as  great 
as  this  separation.  Two  varieties  which  seemed  to  consistently 
produce  photoacous tic  signals  significantly  stronger  than  the 
other  varieties  tested  were  Arkan  and  Coker-916.  In  addition, 
the  phase  of  the  response  from  Coker-916  seemed  to  vary  with 
modulation  frequency  somewhat  differently  than  typically  observed 
with  the  other  varieties.  Based  on  these  observations,  the 
identification  of  a  kernel  from  the  twelve  varieties  tested  as 
being  either  of  these  two  varieties  appears  to  be  immediately 
feasible.  More  detailed  analysis  of  data  obtained  by  the  above- 
described  technique  may  permit  further  identification  as  well. 

One  obvious  variability  factor  in  individual  wheat  kernels 
is  the  kernel  size.  In  an  attempt  to  help  reduce  the  variation 
that  this  might  cause  in  the  photoacous tic  response,  a  normaliza- 
tion on  the  basis  of  kernel  cross-sectional  area  was  attempted. 
The  normalization  procedure  tried  assumed  el ipsoid-shaped  kernels 
of  constant  density,  using  the  functional  dependence  of  cross- 
sectional  area  on  the  kernel  mass  to  normalize  the  data  on  the 
basis  of  the  area  of  the  kernel  being  illuminated.  The  normal- 
ized data,  however,  did  not  appear  to  exhibit  any  less  kernel  to 
kernel  variation  within  a  given  variety,  and  was  often  found  to 
actually  increase  the  scatter.  Such  a  normalization  procedure 
assumes  that  the  spot  size  of  the  incident  e-m  radiation  is 
larger  than  the  kernel's  cross  sectional  area.  Since  for  these 
experiments  the  chamber  depth  was  fairly  shallow  (0.47  cm),  and 
the  previously  described  aperture  restricting  the  cone  of  radia- 
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tlon  entering  the  chamber  was  used,  it  is  possible  that  the  spot 
size  of  significantly  intense  radiation  at  the  kernel  depth  was 
smaller  than  the  size  of  the  wheat  kernels  sampled.  At  the  cham- 
ber depth  which  was  used,  the  3  mm  diameter  aperture  should 
restrict  the  spot  size  to  a  diameter  of  approximately  2.8  mm  at 
the  chamber  floor,  which  is  slightly  larger  than  the  width  of  a 
typical  wheat  kernel.  However,  the  radiation  having  intensity 
greater  than  or  equal  to  half  the  intensity  occurring  at  the 
center  of  the  spot  was  confined  to  a  diameter  of  only  1.1  mm, 
assuming  the  manufacturer's  specifications  for  the  diode  to  be 
accurate.  Thus,  it  is  probable  that  all  significantly  intense 
direct  irradiation  impinged  upon  the  kernel.  Assuming  this  to  be 
the  case,  such  a  normalization  procedure  would  be  counterproduc- 
tive, as  we  found  it  to  be  with  our  data. 

Another  possible  factor  increasing  the  observed  variation  in 
the  photoacoustic  response  among  kernels  of  the  same  variety  is 
the  presence  of  moisture  remaining  in  the  kernel  after  the  drying 
process  which  we  used  was  completed.  As  previously  stated, 
standardization  of  the  moisture  content  is  very  important,  due  to 
its  large  effect  on  the  photoacoustic  signal.  Although  our 
drying  studies  seemed  to  indicate  that  after  fifty  hours  drying, 
the  majority  of  moisture  is  lost  by  the  kernel,  it  is  possible 
that  a  sufficiently  large  amount  of  moisture  is  still  present  in 
the  kernel  to  influence  the  response.  If  this  amount  varies  from 
kernel  to  kernel,  It  should  contribute  to  the  variation  seen  in 
the  photoacoustic  response.   To  reduce  this  problem,  it  might  be 
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necessary  to  use  a  much  longer  drying  time  or  higher  temperature 
to  drive  any  remaining  moisture  from  the  kernel.  Alternatively, 
a  standardization  procedure  bringing  all  kernels  to  the  same 
moisture  content  by  moisture  saturation  might  be  effective,  al- 
though this  would  decrease  the  amplitude  of  the  overall  response. 
Figures  4.35  and  4.36  show  the  average  pho toacous t ic  re- 
sponse for  all  hard  wheat  kernels  and  all  soft  wheat  kernels 
tested.  It  is  seen  that  the  average  amplitude  of  the  response  is 
slightly  greater  on  the  average  for  hard  wheats,  except  for  in 
the  intermediate  portion  of  the  modulation  frequency  range.  The 
greater  amplitude  for  hard  wheat  is  consistent  with  results 
obtained  in  an  earlier  study,  comparing  the  pho toacous tic  re- 
sponse of  a  typical  hard  and  a  typical  soft  variety  (27).  The 
average  phase  of  the  photoacoustic  response  for  hard  wheat  varie- 
ties tends  to  be  slightly  less  than  for  soft  varieties  at  modula- 
tion frequencies  above  approximately  40  Hz,  although  at  frequen- 
cies  below   this    value   there   seems    to  be   little   difference. 
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Figure  4.35.   A  log-log  plot  showing  the  average  amp- 
litude of  the  photoacoustic  response  vs  modulation 
frequency  of  all  of  the  hard  red  winter  wheat  ker- 
nels and  all  of  the  soft  red  winter  wheat  kernels 
tested . 

Legend: 

square  -  average  response,  soft  red  winter 
plus  -  average  response,  hard  red  winter 
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Figure  4.36.   A  plot  showing  the  average  phase  of  the 
photoacoustic  response  vs  log  modulation  frequen- 
cy of  all  of  the  hard  red  winter  wheat  kernels  and 
all  of  the  soft  red  winter  wheat  kernels  tested. 

Legend : 

square  -  average  response,  soft  red  winter 
plus  -  average  response,  hard  red  winter 
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CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

The  results  obtained  from  these  studies  clearly  show  that 
different  kernels  of  the  same  variety  may  produce  photoacoustic 
responses  that  are  significantly  different  from  one  another, 
which  implies  that  there  is  a  certain  amount  of  variation  in 
thermal  properties  between  Individual  kernels,  even  when  the 
kernels  are  of  the  same  variety.  If  moisture  content  is  the  main 
contributor  to  the  observed  variations,  then  more  extensive  sam- 
ple preparation,  either  by  increased  drying  or  some  other  method, 
might  sufficiently  eliminate  differences  In  moisture  content,  and 
allow  the  intrinsic  thermal  properties  of  the  kernels  themselves 
to  dominate  the  photoacoustic  response. 

Despite  the  kernel  to  kernel  variation  which  was  found  to 
occur  In  our  results,  we  did  find  significant  differences  in  the 
photoacoustic  responses  of  Arkan  and  Coker-916,  both  of  which 
provided  signals  relatively  strong  in  comparison  to  the  other 
varieties  tested.  This  seems  to  suggest  that  positive  identifi- 
cation on  the  basis  of  variety  may  be  possible  using  photoacous- 
tic analysis,  if  reduction  in  the  variation  between  kernels  of 
the  same  variety  can  be  achieved.  Although  Coker-916  is  a  soft 
wheat  variety,  and  it  exhibited  a  relatively  strong  photoacoustic 
amplitude,  the  average  amplitude  of  soft  wheats  was  found  to  be 
slightly  lower  than  for  hard  wheats,  and  the  average  phase  seemed 
to  be  slightly  greater  In  the  higher  modulation  frequency  range. 
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This  implies  that  even  with  the  kernel  to  kernel  variations, 
there  seem  to  be  differences  in  the  thermal  properties  of  hard 
and  soft  varieties  which  can  be  detected  using  pho toacous t ic 
analysis  techniques. 

The  kernel  drying  studies  carried  out  in  this  project  seemed 
to  yield  interesting  results,  even  though  they  were  performed 
prior  to  improvements  made  in  both  the  extent  of  the  modulation 
frequency  range  which  was  able  to  be  monitored  and  the  ratio  of 
the  sample  chamber  contribution  to  the  pho toacous tic  signal. 
Since  moisture  has  such  a  pronounced  effect  on  the  photoacous tic 
response,  and  since  frequency  modulation  of  the  incident  e-m 
radiation  permits  profiling  of  a  sample  as  a  function  of  depth, 
detection  of  the  relative  distribution  of  moisture  in  a  kernel  as 
it  is  either  dried  or  allowed  to  absorb  water  seems  to  be  feasi- 
ble. This  may  prove  valuable  for  more  detailed  studies  of  either 
moisture  evaporation  or  absorption  in  grains.  Alternatively, 
photoacoustic  analysis  may  provide  a  rapid  method  for  determining 
the  relative  moisture  content  of  a  single  kernel. 

With  the  addition  of  the  aperture  to  the  photoacoustic  cell, 
the  ratio  of  the  empty  chamber  response  to  the  response  with  a 
wheat  kernel  sample  present  was  greatly  improved.  However,  it  is 
still  possible  that  the  chamber  itself  was  interfering  signifi- 
cantly with  the  sample  response,  especially  at  higher  modulation 
frequencies,  where  the  empty  chamber  and  wheat  kernel  responses 
were  comparable.  A  further  reduction  of  the  empty  chamber  re- 
sponse might  decrease  the  similarity  of  the  responses  obtained 


168 


for    kernels    of    different    varieties,     thus    increasing    the    ability 
to    make    more    clear-cut    distinctions    between    them.        Improvements 
to    the   apparatus   which   could   be   made   include   the  use  of  a   sharply 
focused,    higher   intensity   e-m   radiation   source.      This   would   serve 
to   increase   the   amplitude   of    the   pho toacous tic   response,    thus 
improving    the    signal    to    noise    ratio.       If    the    spot    size    were 
focused    so    that    it   were    smaller    than   the   cross-sectional    area   of 
a    typical    wheat   kernel,     the   chamber's    contribution    to    the    re- 
sponse   should   not    increase   nearly   as   much,    as   only   an   increase   in 
indirect    illumination    caused    by    radiation    scattered    from    the 
sample   would   occur.      Use   of   a    sufficiently    intense   near-Infrared 
laser   could    be    used    to    achieve    this,    although    it    would    most 
likely    require   a   redesign   of    the   pho toacous tic   cell    used    in    these 
experiments.      Variation   of    the   material    used    to    construct    the 
photoacoustic   cell    might   be   advantageous    with   regard    to   decreas- 
ing   the   unwanted    contribution    to    the    photoacoustic    response    from 
the  chamber   itself.     Using  stainless   steel,   for  example,  would 
provide  a   thermal   diffusion   length   approximately   five    times    lower 
than   aluminum   at   a   given   modulation   frequency   (21),    due   primarily 
to    its    much    larger    density.       Since    the    photoacoustic    response 
typically    increases    with    increasing    thermal    diffusion    length,     a 
much   smaller   empty  chamber   response    should    be    expected.      Use    of    a 
photoacoustic   cell    with    transparent  walls    would   essentially   eli- 
minate   absorption    of    radiation    by    the    sample    chamber    Itself,     and 
thus    would   greatly   decrease    the   contribution    to    the   photoacoustic 
response    from    the    chamber.       However,     acoustic    Isolation    and    the 
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ability   of   heat    to   diffuse   from    the   sample  chamber  may  be   compro- 
mised. 

In  photoacoustlc  sample  analysis,  it  Is  the  purpose  of  the 
gas-filled  photoacoustlc  cell  to  allow  for  indirect  detection  of 
temperature  variations  at  the  surface  of  the  sample  being  anal- 
yzed, by  monitoring  the  resultant  pressure  fluctuations  occurring 
in  the  cell.  If  these  temperature  fluctuations  could  be  measured 
more  directly,  It  might  be  possible  to  eliminate  many  of  the 
interfering  effects,  restrictions,  and  variables  associated  with 
the  gas-microphone  configuration,  while  the  same  information 
could  be  obtained.  Using  an  infrared  radiation  detection  device, 
thermal  radiation  emitted  from  the  sample  surface  could  be  moni- 
tored as  a  function  of  the  modulation  frequency,  instead  of 
monitoring  the  acoustic  signal,  allowing  for  direct  detection  of 
temperature  fluctuations  at  the  sample  surface.  Such  an  analysis 
technique  would  offer  several  advantageous  features.  Since  the 
sample  would  not  have  to  be  confined  to  a  gas-filled  cell,  the 
need  for  maintaining  a  helium  gas  atmosphere  would  be  eliminated. 
Also,  the  sample  could  be  easily  positioned  for  illumination,  and 
high-Intensity  and  well-focused  light  sources  such  as  lasers 
could  thus  more  easily  used.  The  problem  with  contributions  to 
the  signal  from  the  sample  chamber  In  our  gas-microphone  photoa- 
coustlc measurements  would  be  eliminated,  as  well  as  would  be 
problems  with  background  noise.  Again,  since  the  temperature 
fluctuations  at  the  sample  surface  provides  the  essential  infor- 
mation   about    the    thermal    properties    of    the    sample,     more    precise 
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results  might  be  obtainable  from  experiments  similar  to  those 
carried  out  in  this  project,  utilizing  a  method  which  incorpor- 
ates some  of  these  improvements. 
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APPENDIX  I 

A  general  flowchart  of  the  data  acquisition  routine  Incor- 
porated into  the  programs  used  by  the  microcomputer  operated  data 
acquisition  system.  Complete  listings  of  the  actual  programs 
used  are  provided  in  Appendices  II  and  III. 
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Canpute  voltage  being 
sait  to  function  gener- 
ator VCO  Input 


Compute  function  generator 
modulation  frequency  as  a 
function  of  VCO  input 


Correct  canjuted  modulation 
frequency  value  for  high 
function  generator  range 


Correct  computed  modulation 
frequency  value  for  high 
function  generator  range 


*© 
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Increase  data  point  counter 
by  I 


Save  modulation  frequency 
in  data  matrix 


Send  new  voltage 
to  function  gen- 
erator VCD  Input 


Execute  delay  loop  to 
allow  lock-In  amplifier 
to  stabilize 


Obtain  50  : 
ate  readings  of 
PIN  diode  output 


Compute  average  of  PIN  dl- 
ode  output  signal  readings 


Convert  average  value  into 
microvolts  and  save  In  data 
matrix 
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Obtain  nunerous  amp-  7 
litude  and  in-phase  / 
camp,  readings  of  / 
microphone  signal       / 


Compute  average  amplitude 
and  in-phase  component  of 
microphone  signal 


Convert  average  values  Into 
microvolts  and  save  in  data 

matrix 


Calculate  phase  angle  from 
amplitude  and  In-phase  com- 
ponent values 
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Increase  value  by  which 
output  value  to  function 
generator  Is  Increased 
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Increase  output  value  to 
be  sent  to  function  gen- 
erator by  appropriate 
Increment 


Send  new  voltage  / 

function  generator       / 
VXD  Input  / 

I 


Increase  sensitiv- 
ity setting  of  lock- 
In  anplif ier  j 


Accept  acknowledgement 
to  proceed  with  acqui-  j 
sltlon  process  at  new  / 
sensitivity  setting    J 


& 


Redefine  sensitivity  of 
lock-in  amplifier  to  accom- 
modate setting  change 


O 
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[  Redefine  appropriate  para- 
I  meters  for  second  portion 
,  of  acquisition  process 


Increase  range  set- 
ting of  function 


Set  output  value  to  func- 
tion generator's  VCO  Input 
at  lowest  value  for  the 
new  range  setting 


Send  output  voltage 
to  function  generator 
VC0  input 


■'  Accept  acknowledgensit 
to  proceed  with  por- 
tion *2  of  acquisition  i 
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APPENDIX  II 

The  acquisition  program  used  in  acquiring  data  for  the  1-50 
Hz  modulation  frequency  range.  The  program  was  written  to  be 
used  with  the  Microsoft  BASICA  Interpreter.  A  discussion  of  the 
program  is  provided  in  Section  3.4,  and  Appendix  I  gives  a  gen- 
eral flowchart  of  the  portion  of  the  program  in  which  the  data  Is 
actually  acquired. 
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10  ADC.SErUP=&H8B4:ALTOUT12.B=&H72E:ALTODT12.BT=&H78A:ALTOUT12.R=&H7E9 

15  ALTODT12.RT=&H847:ALTODT12.S=&H6E5:ALTODT12.ST=&H708:ALTOUT8.B=&H2CF 

20  ALTOUT8 .  BT=&H359 : ALTOUT8 .  R=&H3E9 : ALTOUT8 .  RT=479 : ALTOUT8 .  S=&H24A 

25  ALT0UT8.ST=&H28A:DAC.SETUP=SH1A:INADC12.B=&HCDF:INADC12.B'I^&HD4B 

30  INAIX12.EXTaxDCK^HIX;7:INAI)C12.S^HC70:imDC12.ST=SHC9C:INADC8.B==SH97B 

35  INADC8 .  BT=&H9E2 :  INADC8 .  EXTCLOCK=&HA4C:  INADC8 .  S=&H904 :  INADC8 .  ST=SH933 

40  OUTDAC12.B=&H564:OUTDAC12.B'P=&H5C4:OUTDAC12.R=&H627:OUTDAC12.RT=SH67B 

45  O0TDAC12.S=&H50F:OUTDAC12.ST=SiH538:OOTDAC8.B=4HBC:O0TDAC8.B'P=S,HllD 

50  OUTDAC8 .  R=&H181 : 00TDAC8 . RT=SH1E4 : OUTDAC8 .  S=&H6D: OUTDAC8 .  ST=&H93 

55  SCAN12.B=SHEE5:SCAN12.BT=SHF7A:SCAN12.S=SHE36:SCAN12.ST=&HE82 

60  SCAN8.B=&HB4B:SCAN8.BT=&HBDA:SCAN8.S=SHA9E:SCAN8.ST=&HAE9 

65  DEF  SEG=0 

70  CSEG=(256*PEEK(&H3CB))+PEEK(SH3CA) 

75  DEF  SEG=CSEG 

91  ' 

92  ' 

93  'THE  FIRST  100  LINES  ARE  RESERVED  FOR  LABSTAR  SOFTWARE 

94  "THE  APPLICATION  PROGRAM  STARTS  AT  LINE  100 

95  ' 

96  ' 

97  ' 

98  ' 

99  ' 

100  CLS: PRINT: PRINT: PRINT 

110  PRINT'THE  FOLLOWING  DATA  ACQUISITION  PROGRAM  IS  DESIGNED  TO  OBTAIN  AND" 

120  PRINT'PROCESS  THE  PHOTOACOUSTIC  SIGNAL  FROM  A  LOCK  IN  AMPLIFIER  AS  A" 

130  PRINT" FUNCTION  OF  THE  MODULATION  FREQUENCY  OF  THE  LIGHT  SOURCE  AS" 

140  PRINT"DRIVEN  BY  A  VARIABLE  FREQUENCY  WAVEFORM  GENERATOR." 

150  : 

160  REM  —THE  FOLLOWING  SEVERAL  LINES  DEFINE  CONSTANTS  NECESSARY  FOR  LMPLEMENTAT 

ION  OF  PROG. 

170  : 

175  ZERO$="0000000000" 

180  C0=100:C1=163.654:C2=3.627054:C3=-. 4178307 

190  0VER=256:INCR=15:P0RTI0N=1 

210  CORRECT=10".5:LG=LOG(10) 

215  TMG%=3994:LNG%=50:ADJ=2 

220  DIM  DATMAT(  1000,5)  :DIM  ARRAY%(200)  :DIM  AAYZ(50) 

230  ADDRESS%=&H304 

240  CALL  DAC.SETOP(ADDRESS%) 

250  CHANNEL%=0 

260  VALUE2=128  'ENABLES  INITIAL  MODULATION  FREQUENCY  OF  10.0  HERTZ. 

270  CALL  OUTDAC8.S(CHANNEL%,VALUE%)  'SENDS  INITIAL  VOLTAGE  TO  WFG. 

280  PRINT:  PRINT:  PRINT :  PRINT'TURN  ON  POWER  TO  PHOTOACOUSTIC  SAMPLING  EQUIPTMENT." 

290  PRINT:  PRINT:  PRINT 

300  PRINT"SET  WAVEFORM  GENERATOR  RANGE  AT  100  HZ.  AND  ADJUST  UNTIL  COUNTER  CONSI 

STANTLY" 

310  PRINT"REGISTERS  A  100  HZ.  READING." 
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320  PRINT: INPUT'WHEN  ADJUSTMENT  IS  COMPLETE,  SET  WFG  RANGE  TO  10  HZ.  AND  THEN  PR 

ESS  RETURN.  ";RTN 

330  CLS:PRINT:PRINT:PRINT:PRINT"FROM  THE  METERS  ON  POWER  SUPPLY  DRIVING  LED:":PR 

INT:  PRINT 

340  INPUT'ENTER  DC  VOLTAGE  TO  LED  (IN  VOLTS) ";LEDVLT 

350  PRINT: INPUT'ENTER  CURRENT  TO  LED  (IN  AMPS)";CURRENT 

360  PRINT: PRINT : INPUT'ENTER  RANGE  SETTING  OF  MICROVOLTMETER  (IN  MICROVOLTS)" ;MIC 

RO 

362  CLS:rciNT:PRINT:PRINT:PRINT:PRINT:PRINT:PRINT''SET  WAVEFORM  GENERATOR  TO  1  HZ 

.  RANGE  AND  MAKE  NECESSARY  ADJUSTMENTS  ON  LIA." 

365  PRINT:PRINT:INPUT"ONCE  FINISHED,  PRESS  RETURN";RRN$ 

370  CLS: PRINT:  PRINT: PRINT 

380  PRINT'ENTER  THE  FOLLOWING  LOCK  IN  AMPLIFIER  SETTINGS :":PRINT: PRINT 

390  INPUTTIME  CONSTANT  (IN  SECONDS)" ;TICON 

400  INPUTSENSITIVITY  (IN  MICROVOLTS)" ;SENSI 

410  SENSI$=STR$(SENSI) 

420  IF  LEFT$(SENSI$,2)="  3"  THEN  SENSI=(SENSI/3)*CORRECT 

430  INPUT'PHASE  SETTING  (IN  DEGREES)" ;PHSET$ 

440  CLS:PRINT:PRrNT:PRINT"IS  A  DELAY  TIME  OF  "TIC0N*4"  SECONDS  ADEQUATE  TO  ALLOW 

LOCK  IN  AMPLIFIER" 
450  INPUT'TO  STABILIZE"  ;SEY$ 

460  IF  LEFT$(SEY$,1)="Y"  OR  LEFT$(SEY$,D="y"  THEN  DELAY=TICON*4:GOTO  510 
470  CLS:PRINT:PRINT:PRINT 

480  PRINT'ENTER  DESIRED  DELAY  TIME  (IN  SECONDS)  BETWEEN  MODULATION  FREQUENCY" 
490  PRINT" INCREMENT  AND  SUBSEQUENT  ACQUISITION  OF  SIGNAL  FROM  LOCK  IN  AMP." 
500  PRINT:  INPUT  DELAY 

510  DELAY=DELAY*840  'CONVERTS  LOOP  DELAY  PARAMETER  TO  SECONDS. 
520  as :  PRINT:  PRINT: PRINT 
530  PRINT" IS  50  DATA  POINTS  TAKEN  PER  MODULATION  FREQ.  INCREMENT  SUFFICIENT,  YES 

OR  NO" ; :  INPUT  ESY$ 

540  IF  LEFT$(ESY$,1)="Y"  OR  LEFT$(ESY$,l)="y"  THEN  POINTNO=50:GOTO  570 
550  PRINT: PRINT" HOW  MANY  DATA  POINTS  ARE  TO  BE  TAKEN  AT  EACH  MODULATION  FREQUENC 
Y" 

560  INPUT  POINTNO:IF  INT(POINTNO)OPOINTNO  OR  POINTNCK1  THEN  GOTO  550 
570  LENGTH%=POINTN0 
580  RATE=nCON/POT.NTNO 

590  PRINT: PRINT: PRINT: PRINT'IS  A  SAMPLING  RATE  OF  ";INT(RATE*10000)/10000;"  SECO 
NDS  BETWEEN  ACQUISITION  OF  DATA" 

600  PRINTPOINTS  AT  A  GIVEN  MODULATION  FREQUENCY  SUFFICIENT" ;:  INPUT  AN$ 
610  IF  LEFT$(AN$,1)="Y"  OR  LEFT$(AN$,l)="y"  THEN  GOTO  630 
620  PRINT: INPUT'ENTER  DESIRED  SAMPLING  RATE  (IN  SECONDS)" ;RATE 
630  BETWEEN=RATE/. 000005 

640  IF  BETWEEN>65536!  THEN  PRINT:  PRINT'DESIRED  SAMPLING  RATE  IS  TOO  LARGE,  IT  CA 
NNOT  EXCEED  0.327  SECONDS":PRINT:PRINT:GOTO  620 
650  IF  BETWEEN>32767  THEN  BETWEEN=BETWEEN-65536 ! 
660  TTMING%=INT(BETWEEN) 
670  as : PRINT:  PRINT:  PRINT 
680  INPUT"  ENTER  SAMPLE  IDENTIFICATION  NAME" ;  SAMPLE$ 
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690  PRINT:PRINT:PRINT"ENTER  PRESENT  DATE:  MONTH/DAY/YEAR.";  :INPUT  TODAY$ 

700  CLS: PRINT: PRINT: PRINT" SAMPLE  IDENTIFICATION  NAME  IS:  ";SAMPLE$ 

710  PRINT :PRINT"TODAY  IS  ";TODAY$: PRINT: PRINT" YES  OR  NO";: INPUT  YYY$ 

720  IF  LEFT$(YYY$,1)<>"Y"  AND  LEFT$(YYY$, !)<>"y"  THEN  GOTO  670 

730  CLS:PRINT:PRINT:PRINT:INPOT"WHEN  READY  TO  BEGIN,  PRESS  RETURN. ";START$ 

740  IF  START$="STOP"  THEN  GOTO  1910  'SENDS  TO  END  OF  PROGRAM. 

750  LPRINT:LPRINT:LPRINT"DATA  FOR  ";SAMPLE$ 

760  LPRINT:LPRINT"DATE  DATA  ACQUIRED:  ";TODAY$ 

770  CLS: PRINT :PRINT"DATA  FOR  ";SAMPLE$ 

780  PRINT: PRINT" CHOP.  FREQ.";TAB(15) ;"SIG.  AMP.";TAB(30) ;"A  STN(p.a.)";TAB(45);" 

PHASE  ANGLE" ;TAB( 60); "PIN  DIODE  VLT.": PRINT 

790  LPRINT:LPRINT:LPRINT"CHOP.  FREQ.";TAB(15) ;"SIG.  AMP.";TAB(30) ;"A  SIN(p.a.)"; 

TAB(45);  "PHASE  ANGLE" ;TAB(60) ; "PIN  DIODE  VLT.":LPRINT 

800  : 

810  REM  —DATA  ACQUISITION  PORTION  OF  PROGRAM  BEGINS  HERE. 

820  : 

830  WFGINP=INT(((VALUE%-128)*.03906)*100)/100  'COMPUTES  VOLTAGE  SENT  TO  WFG. 

840  CHFREQ=INT( (C0+OlWGINIH<;2*WFGINP-2+C3<OTGINP~3)*10)/10 

850  IF  PORTION=l  THEN  CHFREQ=INT(CHFREQ)/100:GOTO  890 

855  IF  PORTION=2  THEN  CHFREQ=INT(CHFREQ)/10 

860  XYZ=INT(CHFREQ/HICUT) 

870  IF  CHFREQ<XYZ*60+LOCUT  OR  CHFREQ>XYZ*60+HICUT  THEN  GOTO  890 

880  GOTO  1440  'SENDS  TO  POINT  WHERE  FREQ.  IS  ALLOWED  TO  INCREASE. 

890  C0UNT=C0UNT+1 

900  DATMAT(COUNT,l)=CHFREQ 

910  CHANNEL%=0 

920  ADDRESS%=&H304 

930  CALL  DAC.SETUP(ADDRESS%) 

940  CALL  OUTDAC8.S(CHANNEL%,VALUE%)  'SENDS  DESIRED  VOLTAGE  TO  WFG. 

950  FOR  1=1  TO  DELAY:  NEXT  I  'DELAY  LOOP  TO  ALLOW  LOCK  IN  AMP.  TO  STABILIZE. 

960  ADDRESS%=&H300 

970  CALL  ADC.SETUP(ADDRESS%) 

971  CHANNEL%-1 

973  CALL  INADC8.B(CHANNEL%,TMG%,LNG%,AAY%(1)) 

974  PNVLT=0 

975  FOR  1=1  TO  LNG% 

976  PNVLT=PNVLT+AAY%(I) 

977  NEXT  I 

978  PNVLT=PNVLT/LNG% 

979  DATMAT(COUNT,5)=PNVLT*ADJ*.0019608*MICRO 

980  BCHANNEL%=4:ECHANNEL%=5 

1000  CALL  SCM8.B(IiCHANNEL%,ECHANNEL%,TIMING%,LENGTH%,ARRAn(l)) 

1010  FOR  1=1  TO  POINTNO*2  STEP  2 

1020  ASUM=ASUM+ARRAY%(I) 

1030  SSUH=SSUM+ARRAY%(I+1) 

1035  NEXT  I 

1040  AV=(ASUM/POINTOO)*. 0019608  'CONVERTS  BINARY  SIGNAL  TO  LIA  OUTPUT  VOLTAGE. 

1050  SV=(SSUM/POINTNO)*. 0019608 


184 


1060  DATMAT(COUNT,2)=AV*SENSI  'CONVERTS  VOLT.  TO  LIA  METER  READING  (IN  MICROV.). 
1070  DATMAT(COUNT,3)=SV*SENSI 

1080  X=DATMAT(COUNT,3)/DATMAT(COUNT,2)  'FINDS  SINE  OF  PHASE  ANGLE 
1085  IF  X>1  THEN  DATMAT( COUNT, 4)=90.99999:GOTO  1100 

1090  DATMAT(COUNT,4)=ATN(X/SQR(-X*X+l))*57. 29578  'FINDS  PHASE  ANGLE  (IN  DEG.) 
1100  ASUM=0:SSUM=0 
1110  FOR  1=1  TO  5 

1120  PRINT  DATMAT(C0UNT,I);TAB(15*I); 
1130  NEXT  I 
1140  PRINT 
1150  FOR  1=1  TO  5 

1160  LPRINT  DATMAT(COUNT,I);TAB(15*I); 
1170  NEXT  I 
1180  LPRINT 

1230  IF  AV/.0019608>80  OR  ABC$="C"  THEN  GOTO  1440  'CHECKS  IF  SENSITIVITY  SHOULD 
BE  INCREASED. 

1240  FOR  1=1  TO  4:BEEP:FOR  J=l  TO  90:NEXT  J:NEXT  I 
1250  IF  MES=1  THEN  GOTO  1310 

1260  PRINT:PRINT:PRINT:PRINT"SIGNAL  AMP.  HAS  DECREASED  TO  A  POINT  WHERE  THE  SENS 
ITTVITY  CAN  BE  INCREASED." :MES=1 

1270  PRINT"IF  YOU  WISH  TO  MAINTAIN  THE  SAME  SENSITIVITY  THROUGHOUT  THE  EXPERIMEN 
T,  ENTER" 

1280  PRINT'C  FOR  CONTINUE.  OTHERWISE  SIMPLY  ENTER  A  BLANK  LINE.": INPUT  ABC$ 
1290  IF  ABC$="C"  THEN  GOTO  1440 

1300  IF  VALUE%+INCR>=OVER  AND  PORTION=2  THEN  GOTO  1650 
1310  VALUE*=VALUE%+INCR 
1320  ADDRESS%=&H304 
1330  CALL  DAC.SETUP(ADDRESS%) 
1340  CHANNEL%=0 

1350  CALL  OUTDAC8.S(CHANNEL%, VALUE!) 

1360  PRINT: PRINT: PRINT" SENSITIVITY  MAY  NOW  BE  INCREASED.": PRINT: PRINT 
1370  SENSI=SENSI/CORRECT 

1380  PRINT"AFTER  SENSITIVITY  SETTING  IS  INCREASED,  AND  LOCK  IN  AMP.  HAS  STABILIZ 
ED," 

1390  INPUT"PRESS  RETURN.  ";RTN$ 
1400  JKL$=STR$(INT(SENSI+.5)) 

1405  IF  LEFT$(JKL$,2)="  3"  THEN  JKL$=LEFT$(JKL$,2)+LEFT$(ZER0$,LEN(JKL$)-2) 
1410  PRINT: PRINT"NEW  SENSITIVITY  IS";JKL$;"  MICROVOLTS.": PRINT: PRINT 
1420  LPRINT: LPRINT; "NEW  SENSITIVITY  IS";JKL$;"  MICROVOLTS.": LPRINT 
1430  IF  VALUEKOVER  THEN  GOTO  830 

1440  VALUE%=VALUE%+INCR  'ALLOWS  INCREMENT  IN  VOLTAGE  SENT  TO  WFG. 
1450  IF  VALUEKOVER  THEN  GOTO  830 
1460  IF  P0RTI0N=2  THEN  GOTO  1650 

1470  FOR  1=1  TO  5:BEEP:FOR  J=l  TO  130:NEXT  J:NEXT  I 
1480  PORTION-2 
1490  INCR=6:OVER=200 
1500  L0CUT=57.5:HICUT=62.5 
1510  PRINT: PRINT: PRINT'PORTION  #1  OF  DATA  ACQUISITION  IS  COMPLETE." 
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1520  PRINT: PRINT" CHANGE  APPROPRIATE  SETTINGS  ON  UFG  AND  LIA." 

1530  PRINT: PRINT'WHEN  READY  TO  CONTINUE  SAMPLING,  HIT  RETURN:" 

1540  PRINT'OR,  IF  YOU  WISH  TO  DISCONTINUE  SAMPLING,  ENTER  STOP."; 

1550  CHANNEL%=0 

1560  VALUE%=128  'ENABLES  MODULATION  FREQUENCY  OF  10  HZ. 

1570  CALL  OUTDAC8.S(CHANNEL%,VALUE%)  'SENDS  NEW  VOLT.  SIG.  TO  WFG  FOR  PORTION  2. 

1580  INPUT  GO$ 

1590  IF  GO$=""  THEN  GOTO  830 

1600  IF  GO$="STOP"  THEN  GOTO  1650 

1610  GOTO  1580 

1620  : 

1630  REM  —DATA  ACQUISITION  PORTION  OF  PROGRAM  ENDS  HERE. 

1640  : 

1650  LPRLNT:LPRINT:LPRINT 

1660  LPRLNT"NUMBER  OF  DATA  POINTS  TAKEN:  ";COUNT 

1670  LPRINT" INITIAL  LIA  SENSITIVITY  SETTING:  ";SENSI$;"  MICROVOLTS." 

1680  LPRINT'LIA  TIME  CONSTANT  SETTING:  ";TICON; "SECONDS." 

1690  LPRINTDELAY  TIME  USED:  ";DELAY/840;"SECONDS." 

1695  LPRINT'NUMBER  OF  SAMPLINGS  PER  DATA  POINT:  ";POINTNO 

1700  LPRINT" SAMPLING  RATE:  ";RATE+.000033;"SECONDS." 

1710  LPRINTVOLTAGE  TO  LED  (AT  10  HZ.):  ";LEDVLT;"  VOLTS." 

1720  LPRINT"CURRENT  TO  LED  (AT  10  HZ.):  ";CURRENT;"  AMPS." 

1730  LPRINT'RANGE  SETTING  OF  MICROVOLTMETER:  ";MICRO;"  MICROVOLTS." 

1740  LPRINT: LPRINT: LPRINT 

1750  PRINT: PRINT: PRINT'SAMPLING  IS  COMPLETED." 

1760  : 

1770  REM  —THE  FOLLOWING  LINES  ALLOW  FOR  FILING  OF  ACQUIRED  DATA. 

1780  : 

1790  PRINT:PRINT:PRINT:INPUT"DO  YOU  WISH  TO  FILE  ACQUIRED  DATA";Y$ 

1800  IF  LEFT$(Y$,1)="N"  OR  LEFT$(Y$,l)="n"  THEN  GOTO  1910 

1810  CLS:PRINT:PRINT:PRINT"DO  YOU  WISH  TO  FILE  ACQUIRED  DATA  UNDER  FILENAME  ";SA 

MPLE$;: INPUT  YE$: PRINT: PRINT 

1820  IF  LEFT$(YE$,1)="Y"  OR  LEFT$(YE$,l)="y"  THEN  NME$=SAMPLE$:GOTO  1850 

1830  PRINT: PRINT: PRINT" INPUT  DESIRED  FILENAME  FOR  ";SAMPLE$ 

1840  PRINT: INPUT  NME$ 

1850  OPEN  "0",#1,NME$+".PRN" 

1860  FOR  1=1  TO  COUNT 

1870  PRINT#1,DATMAT(I,1);LOG(DATMAT(I,1))/LG;DATMAT(I,2);LOG(DATMAT(I,2))/LG;DAT 

MAT(I,4);DATMAT(I,5) 

1880  NEXT  I 

1890  CLOSE 

1900  CLS:PRrNT:PRINT:PRINT:PRINT"FILING  IS  COMPLETE." 

1910  END 
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APPENDIX  III 

The  acquisition  program  used  In  acquiring  data  for  the  10- 
1000  Hz  modulation  frequency  range.  The  program  was  written  to 
be  used  with  the  Microsoft  BASICA  Interpreter.  A  discussion  of 
the  program  Is  provided  in  Section  3.4,  and  Appendix  I  gives  a 
general  flowchart  of  the  portion  of  the  program  in  which  the  data 
is  actually  acquired. 
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10  ADC.SETnP=4H8B4:ALTOUT12.B=&H72E:ALTOUT12.BT=&H78A:ALTOUT12.R=&H7E9 

15  ALTOUT12.RT=SH847:ALTOnT12.S=&H6E5:ALTOUT12.S'I>&H708:ALTOUT8.B=&H2CF 

20  ALTOUT8 .  BT=&H359 :  ALT0UT8 .  R=&H3E9 : ALTOUT8 .  RT=479 : ALTOUT8 .  S=&H24A 

25  ALTOUT8.ST=&H28A:DAC.SETUP=SH1A:INADC12.B=&HCDF:INADC12.BT=&HD4B 

30  INAK;i2.EXTCMCK==4HDC7:INADC12.S=&HC70:INADC12.S'P=&HC9C:INADC8.B==&H97B 

35  INADC8 .  BT=&H9E2 :  INADC8 .  EXTCL0CK=«HA4C :  INADC8 .  S=&H904 :  INADC8 .  ST=&H933 

40  O0TDAC12.B=&H564:OUTDAC12.BT=iH5C4:OUTDAC12.R=&H627:OUTDAC12.RT=&H67B 

45  OOTDAC12.S=&H50F:OUTDAC12.ST=&H538:OUTDAC8.B=&HBC:OOTDAC8.B'P=&H11D 

50  00TDAC8.R=&H181  :OUTDAC8.RT^HlE4:OUTDAC8.S=&H6D:OUTDAC8.ST=&H93 

55  SCAN12.B=&HEE5:SCAN12.BI^SHF7A:SCAN12.S=&HE36:SCAN12.ST=&HE82 

60  SCAN8.B=4HB4B:SCAN8.BT=SHBDA:SCAN8.S=SHA9E:SCAN8.S'f=&HAE9 

65  DEF  SEG=0 

70  CSEG=(256*PEEK(SH3CB))+PEEK(5.H3CA) 

75  DEF  SEG=CSEG 

91  ' 

92  ' 

93  "THE  FIRST  100  LINES  ARE  RESERVED  FOR  LABSTAR  SOFTWARE 

94  'THE  APPLICATION  PROGRAM  STARTS  AT  LINE  100 

95  ' 

96  ' 

97  ' 

98  ' 

99  ' 

100  as : PRINT: PRINT: PRINT 

110  PRINTTHE  FOLLOWING  DATA  ACQUISITION  PROGRAM  IS  DESIGNED  TO  OBTAIN  AND" 

120  PRINT"PROCESS  THE  PHOTOACCOUSTIC  SIGNAL  FROM  A  LOCK  IN  AMPLIFIER  AS  A" 

130  PRINT"FUNCTION  OF  THE  MODULATION  FREQUENCY  OF  THE  LIGHT  SOURCE  AS" 

140  PRINT"DRIVEN  BY  A  VARIABLE  FREQUENCY  WAVEFORM  GENERATOR." 

150  : 

160  REM  —THE  FOLLOWING  SEVERAL  LINES  DEFINE  CONSTANTS  NECESSARY  FOR  IMPLEMENTAT 

ION  OF  PROG. 

170  : 

175  ZERO$="0000000000" 

180  C0=100:C1=163.654:C2=3.627054:C3=-.4178307 

190  OVER=256:INCR=6:PORTION=l 

195  FINC=1:IFREQ=250:STFREQ=250 

200  LOCUT=57.5:HICUT=62.5 

210  CORRECT=10~.5:LG=LOG(10) 

215  TMG%=3994:LNG%=50:ADJ=2 

220  DIM  DATMAT(1000,5):DIM  ARRAY%(200)  :DIM  AAY%(50) 

230  ADDRESS%=SH304 

240  CALL  DAC.SETUP(ADDRESS%) 

250  CHANNEL%=0 

260  VALUK-128  'ENABLES  INITIAL  MODULATION  FREQUENCY  OF  10.0  HERTZ. 

270  CALL  OUTDAC8.S(CHANNEL%,VALUE%)  'SENDS  INITIAL  VOLTAGE  TO  WFG. 

280  PRINT:PRINT:PRINT:PRINT"TURN  ON  POWER  TO  PHOTOACOUSTIC  SAMPLING  EQUIPTMENT." 

290  PRINT: PRINT: PRINT 

300  PRINT'SET  WAVEFORM  GENERATOR  RANGE  AT  100  HZ.  AND  ADJUST  UNTIL  COUNTER  CONSI 
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STANTLY" 

310  PRINT'REGISTERS  A  100  HZ.  READING." 

320  PRINT :INPOT"WHEN  ADJUSTMENT  IS  COMPLETE,  SET  WFG  RANGE  TO  10  HZ.  AND  THEN  PR 

ESS  RETURN." ;RTN 

330  CLS:PRINT:PRINT:PRINT:PRINT"FROM  THE  METERS  ON  POWER  SUPPLY  DRIVING  LED:":PR 

INT:  PRINT 

340  INPUT" ENTER  DC  VOLTAGE  TO  LED  (IN  VOLTS )";LEDVLT 

350  PRINT: INPUT'ENTER  CURRENT  TO  LED  (IN  AMPS)"; CURRENT 

360  PRINT : PRINT : INPUT'ENTER  RANGE  SETTING  OF  MICROVOLTMETER  (IN  MICROVOLTS)" ;MIC 

RO 

370  as : PRINT: PRINT: PRINT 

380  PRINTENTER  THE  FOLLOWING  LOCK  IN  AMPLIFIER  SETTINGS:": PRINT: PRINT 

390  INPUT'TIME  CONSTANT  (IN  SECONDS)" ;TICON 

400  INPUT'SENSITIVITY  (IN  MICROVOLTS)";SENSI 

410  SENSI$=STR$(SENSI) 

420  IF  LEFT$(SENSI$,2)="  3"  THEN  SENSI=(SENSI/3)*C0RRECT 

430  INPUT'PHASE  SETTING  (IN  DEGREES) ";PHSET$ 

440  as : PRINT: PRINT: PRINT" IS  A  DELAY  TIME  OF  "TIC0N*4"  SECONDS  ADEQUATE  TO  ALLOW 

LOCK  IN  AMPLIFIER" 
450  INPUTTO  STABILIZE"  ;SEY$ 

460  IF  LEFT$(SEY$,1)="Y"  OR  LEFT$(SEY$,l)="y"  THEN  DELAY=TICON*4:GOTO  510 
470  as :  PRINT:  PRINT:  PRINT 

480  PRINTENTER  DESIRED  DELAY  TIME  (IN  SECONDS)  BETWEEN  MODULATION  FREQUENCY" 
490  PRINT" INCREMENT  AND  SUBSEQUENT  ACQUISITION  OF  SIGNAL  FROM  LOCK  IN  AMP." 
500  PRINT:  INPUT  DELAY 

510  DELAY=DELAY*840  'CONVERTS  LOOP  DELAY  PARAMETER  TO  SECONDS. 
520  as : PRINT: PRINT: PRINT 
530  PRINT"IS  50  DATA  POINTS  TAKEN  PER  MODULATION  FREQ.  INCREMENT  SUFFICIENT,  YES 

OR  NO";: INPUT  ESY$ 

540  IF  LEFT$(ESY$,1)="Y"  OR  LEFT$(ESY$,l)="y"  THEN  POINTNO=50:GOTO  570 
550  PRINT:  PRINT'HOW  MANY  DATA  POINTS  ARE  TO  BE  TAKEN  AT  EACH  MODULATION  FREQUENC 
Y" 

560  INPUT  POINTNO:IF  INT(POINTNO)OPOINTNO  OR  POINTNO<1  THEN  GOTO  550 
570  LENGTH%=POINTNO 
580  RATE=TICON/POINTNO 

590  PRINT: PRINT: PRINT: PRINT" IS  A  SAMPLING  RATE  OF  ";INT(RATE*10000)/10000;"  SECO 
NDS  BETWEEN  ACQUISITION  OF  DATA" 

600  PRINT"  POINTS  AT  A  GIVEN  MODULATION  FREQUENCY  SUFFICIENT" ;:  INPUT  AN$ 
610  IF  LEFT$(AN$,1)="Y"  OR  LEFT$(AN$, l)="y"  THEN  GOTO  630 
620  PRINT:  INPUT" ENTER  DESIRED  SAMPLING  RATE  (IN  SECONDS)"  ; RATE 
630  BETUEEN=RATE/. 000005 

640  IF  BETWEEN>65536!  THEN  PRINT:  PRINTDESIRED  SAMPLING  RATE  IS  TOO  LARGE,  IT  CA 
NNOT  EXCEED  0.327  SECONDS": PRINT: PRINT:GOTO  620 
650  IF  BETWEEN>32767  THEN  BETWEEN=BETWEEN-65536! 
660  TIMING%=INT(BETWEEN) 
670  aS:PRINT:PRINT:PRINT 

680  INPUT'ENTER  SAMPLE  IDENTIFICATION  NAME" ;  SAMPLE$ 
690  PRINT:PRINT:PRINT"ENTER  PRESENT  DATE:  MONTH/DAY/ YEAR." ;  :INPUT  T0DAY$ 
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700  CLS:PRINT:PRINT:PRINT"SAMPLE  IDENTIFICATION  NAME  IS:  ";SAMPLE$ 

710  PRINT :PRINT"TODAY  IS  ";TODAY$: PRINT: PRINT" YES  OR  NO";: INPUT  YYY$ 

720  IF  LEFT$(YYY$,1)<>"Y"  AND  LEFT$(YYY$,l)<>"y"  THEN  GOTO  670 

730  CLS: PRINT: PRINT: PRINT :INPUT"WHEN  READY  TO  BEGIN,  PRESS  RETURN. "; START$ 

740  IF  START$="STOP"  THEN  GOTO  1910  'SENDS  TO  END  OF  PROGRAM. 

750  LPRINT:LPRINT:LPRINT"DATA  FOR  ";SAMPLE$ 

760  LPRINT:LPRINT"DATE  DATA  ACQUIRED:  ";TODAY$ 

770  CLS: PRINT :PRINT"DATA  FOR  ";SAMPLE$ 

730  PRINT:PRINT"CHOP.  FREQ.";TAB(15)  ;"SIG.  AMP.";TAB(30);"A  SIN(p.a.)";TAB(45);' 

PHASE  ANGLE";TAB(60);"PIN  DIODE  VLT.":PRINT 

790  LPRINT:LPRINT:LPRrNT"CHOP.  FREQ.";TAB(15) ;"SIG.  AMP.";TAB(30);"A  SIN(p.a.)" 

TAB(45);  "PHASE  ANGLE" ;TAB(60); "PIN  DIODE  VLT.":LPRINT 

800  : 

310  REM  —DATA  ACQUISITION  PORTION  OF  PROGRAM  BEGINS  HERE. 

820  : 

830  WFGINP=INT(((VALUE%-128)*.03906)*100)/100  'COMPUTES  VOLTAGE  SENT  TO  WFG. 

840  OlFREQ=INT((CCH<:in^GINPH<;2nTOINP"2+C3*WFGrNP"3)*10)/10 

850  IF  P0RTI0N=1  THEN  CHFREQ=INT(CHFREQ)/10 

860  XYZ=INT(CHFREQ/HICUT) 

870  IF  CHFREQ<XYZ*60+L0CUT  OR  CHFREQ>XYZ*60+HICUT  THEN  GOTO  890 

880  GOTO  1440  'SENDS  TO  POINT  WHERE  FREQ  IS  ALLOWED  TO  INCREASE. 

890  COUNT=COUNT+l 

900  DATMAT(COUNT,l)=CHFREQ 

910  CHANNEU=0 

920  ADDRESS%=&H304 

930  CALL  DAC.SETUP(ADDRESS%) 

940  CALL  OUTDAC8.S(CHANNEL%,VALUE%)  'SENDS  DESIRED  VOLTAGE  TO  WFG. 

950  FOR  1=1  TO  DELAY:  NEXT  I  'DELAY  LOOP  TO  ALLOW  LOCK  IN  AMP.  TO  STABILIZE. 

960  ADDRESS%=&H300 

970  CALL  ADC.SETUP(ADDRESS%) 

971  CHANNEL%=1 

973  CALL  INADC8.B(CHANNEU,TMG%,LNG%,AAY%(1)) 

974  PNVLT=0 

975  FOR  1=1  TO  LNG% 

976  PNVLT=PNVLT+AAY%(I) 

977  NEXT  I 

978  PNVLT=PNVLT/LNG% 

979  DATMAT(COUNT,5)=PNVLT*ADJ*.0019608*MICRO 

980  BCHANNEL%=4:ECHANNEL%=5 

990  CALL  SCAN8.B(BCHANNEL%,ECHANNEL%,TIMINa,LENGTH%,ARRAY%(l)) 

1000  FOR  1=1  TO  P0INTN0*2  STEP  2 

1010  ASUM=ASUM+ARRAY%(l) 

1020  SSUM=SSUM+ARRAY%(I+1) 

1030  NEXT  I 

1040  AV=(ASUM/POINTNO)*. 0019608  'CONVERTS  BINARY  SIGNAL  TO  LIA  OUTPUT  VOLTAGE. 

1050  SV=(SSUM/POINTNO)*. 0019608 

1060  DATMAT(C0UNT,2)=AV*SENSI  'CONVERTS  VOLT.  TO  LIA  METER  READING  (IN  MICROV.). 

1070  DATMAT(COUNT,3)=SV*SENSI 
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1080  X=DATMAT(C0UNT,3)/DATMAT(C0UNT,2)  'FINDS  SINE  OF  PHASE  ANGLE 
1085  IF  X>1  THEN  DATMAT(COUNT,4)=90.99999:GOTO  1100 

1090  DATMAT(COUNT,4)=ATN(X/SQR(-X*X+l))*57. 29578  'FINDS  PHASE  ANGLE  (IN  DEG.) 
1100  ASUM=0:SSUM=0 
1110  FOR  1=1  TO  5 

1120  PRINT  DATMAT(C0UNT,I);TAB(15*I); 
1130  NEXT  I 
1140  PRINT 
1150  FOR  1=1  TO  5 

1160  LPRINT  DATMAT(COUNT,I);TAB(15*I); 
1170  NEXT  I 
1180  LPRINT 

1190  IF  CHFREOXSTFREQ  THEN  GOTO  1220  'THESE  3  LINES  ALLOW  FOR  INCR.  IN  FREQ.  INC 
REMENT  IN  INTERVALS  OF  IFREQ. 
1200  STFREQ=STFREO+IFREQ 
1210  INCR=INCR+FINC 

1220  IF  VALUE%+INCR>=OVER  AND  PORTION=2  THEN  GOTO  1650 

1230  IF  AV/.0019608>80  OR  ABC$="C"  THEN  GOTO  1440  'CHECKS  IF  SENSITIVITY  SHOULD 
BE  INCREASED. 

1240  FOR  1=1  TO  4:BEEP:F0R  J=l  TO  90:NEXT  J:NEXT  I 
1250  IF  MES=1  THEN  GOTO  1310 

1260  PRINT:PRINT:PRINT:PRINT"SIGNAL  AMP.  HAS  DECREASED  TO  A  POINT  WHERE  THE  SENS 
ITIVITY  CAN  BE  INCREASED. ":MES=1 

1270  PRINTIF  YOU  WISH  TO  MAINTAIN  THE  SAME  SENSITIVITY  THROUGHOUT  THE  EXPERIMEN 
T,  ENTER" 

1280  PRINT'C  FOR  CONTINUE.  OTHERWISE  SIMPLY  ENTER  A  BLANK  LINE.": INPUT  ABC$ 
1290  IF  ABC$="C"  THEN  GOTO  1440 

1300  IF  VALUE%+INCR>=OVER  AND  PORTION=2  THEN  GOTO  1650 
1310  VALUE%=VALUE%+INCR 
1320  ADDRESS%=&H304 
1330  CALL  DAC.SETUP(ADDRESS%) 
1340  CHANNEL%=0 

1350  CALL  OUTDAC8.S( CHANNELS, VALUES) 

1360  PRINT:PRINT:PRINT"SENSITIVITY  MAY  NOW  BE  INCREASED.": PRINT: PRINT 
1370  SENSI=SENSI/CORRECT 

1380  PRINT"AFTER  SENSITIVITY  SETTING  IS  INCREASED,  AND  LOCK  IN  AMP.  HAS  STABILIZ 
ED," 

1390  TNPUT'PRESS  RETURN. ";RTN$ 
1400  JKL$=STR$(INT(SENSI+.5)) 

1405  IF  LEFT$(JKL$,2)="  3"  THEN  JKL$=LEFT$(JKL$,2)+LEFT$(ZER0$,LEN(JKL$)-2) 
1410  PRINT: PRINT"NEW  SENSITIVITY  IS"jJKL$;"  MICROVOLTS.": PRINT: PRINT 
1420  LPRINT:LPRINT;"NEW  SENSITIVITY  IS";JKL$;"  MICROVOLTS.": LPRINT 
1430  IF  VALUEKOVER  THEN  GOTO  830 

1440  VALUE%=VALUE%+INCR  'ALLOWS  INCREMENT  IN  VOLTAGE  SENT  TO  WFG. 
1450  IF  VALUE%<OVER  THEN  GOTO  830 
1460  IF  P0RTI0N=2  THEN  GOTO  1650 

1470  FOR  1=1  TO  5:BEEP:F0R  J=l  TO  130:NEXT  J:NEXT  I 
1480  P0RTI0N=2 
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1490  INCR=2:OVER=256 

1500  LOCUT=55:HICUT=65 

1510  PRINT: PRINT :PRINT"PORTION  111   OF  DATA  ACQUISITION  IS  COMPLETE." 

1520  PRINT:PRINT"CHANGE  APPROPRIATE  SETTINGS  ON  WFG  AND  LIA." 

1530  PRINT :PRINT"WHEN  READY  TO  CONTINUE  SAMPLING,  HIT  RETDRN:" 

1540  PRINT"0R,  IF  YOU  WISH  TO  DISCONTINUE  SAMPLING,  ENTER  STOP."; 

1550  CHANNEL%=0 

1560  VALUE%=128  'ENABLES  MODULATION  FREQUENCY  OF  100  HZ. 

1570  CALL  0UTDAC8.S(  CHANNELS,  VALUE%)  'SENDS  NEW  VOLT.  SIG.  TO  WFG  FOR  PORTION  2. 

1580  INPUT  GO$ 

1590  IF  GO$=""  THEN  GOTO  830 

1600  IF  GO$="STOP"  THEN  GOTO  1650 

1610  GOTO  1580 

1620  : 

1630  REM  —DATA  ACQUISITION  PORTION  OF  PROGRAM  ENDS  HERE. 

1640  : 

1650  LPRINT:LPRINT:LPRINT 

1660  LPRINT"NUMBER  OF  DATA  POINTS  TAKEN:  ";COUNT 

1670  LPRINT" INITIAL  LIA  SENSITIVITY  SETTING:  ";SENSI$;"  MICROVOLTS." 

1680  LPRINT'LIA  TIME  CONSTANT  SETTING:  ";TICON; "SECONDS." 

1690  LPRINT'DELAY  TIME  USED:  ";DELAY/840; "SECONDS." 

1695  LPRINT'NUMBER  OF  SAMPLINGS  PER  DATA  POINT:  ";POINTNO 

1700  LPRINT" SAMPLING  RATE:  ";RATE+.000033;"SECONDS." 

1710  LPRBTTVOLTAGE  TO  LED  (AT  10  HZ.):  ";LEDVLT;"  VOLTS." 

1720  LPRINT" CURRENT  TO  LED  (AT  10  HZ.):  ";CURRENT;"  AMPS." 

1730  LPRINT"RANGE  SETTING  OF  MICROVOLTMETER:  ";MICRO;"  MICROVOLTS." 

1740  LPRINT: LPRINT: LPRINT 

1750  PRINT:PRINT:PRINT"SAMPLING  IS  COMPLETED." 

1760  : 

1770  REM  —THE  FOLLOWING  LINES  ALLOW  FOR  FILING  OF  ACQUIRED  DATA. 

1780  : 

1790  PRINT: PRINT: PRINT: INPUT'DO  YOU  WISH  TO  FILE  ACQUIRED  DATA";Y$ 

1800  IF  LEFT$(Y$,1)="N"  OR  LEFT$(Y$,l)-"n"  THEN  GOTO  1910 

1810  CLS:PRINT:PRINT:PRINT"DO  YOU  WISH  TO  FILE  ACQUIRED  DATA  UNDER  FILENAME  ";SA 

MPLE$ ; :  INPUT  YE$ :  PRINT :  PRINT 

1820  IF  LEFT$(YE$,1)="Y"  OR  LEFT$(YE$,l)="y"  THEN  NME$=SAMPLE$:GOTO  1850 

1830  PRINT: PRINT: PRINT" INPUT  DESIRED  FILENAME  FOR  ";SAMPLE$ 

1840  PRINT: INPUT  NME$ 

1850  OPEN  "0",#1,NME$+".PRN" 

1860  FOR  1=1  TO  COUNT 

1870  PRINT)l'l,DATMAT(I,l);LOG(DATMAT(I,l))/LG;DATMAT(I,2);LOG(DATMAT(I,2))/LG;DAT 

MAT(I,4);DATMAT(I,5) 

1880  NEXT  I 

1890  CLOSE 

1900  CLS:PRINT:PRINT:PRINT:PRINT"FrLING  IS  COMPLETE." 

1910  END 
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ABSTRACT 

A  phenomenon  commonly  referred  to  as  the  photoacous tic  ef- 
fect is  known  to  occur  when  an  intensity  modulated  source  of 
electromagnetic  radiation  impinges  on  a  solid  sample  enclosed  in 
a  gas-filled  chamber.  Upon  absorption  of  the  incident  radiation, 
an  acoustic  signal  is  produced  in  the  gas  with  a  frequency  equal 
to  the  modulation  frequency.  The  amplitude  and  phase  of  the 
signal  depends  on  the  thermal,  and  in  many  cases,  the  optical 
properties  of  the  sample,  thus  allowing  the  effect  to  be  used  as 
a    sample   analysis    technique. 

A  series  of  experiments  was  performed  in  this  research 
project  to  study  the  photoacous tic  responses  of  individual  wheat 
kernels.  The  technique  used  involved  varying  the  modulation  fre- 
quency of  the  incident  radiation  over  a  1-1000  Hz  frequency 
range,  and  monitoring  the  photoacous  tic  response.  The  acoustic 
signal  was  detected  with  a  sensitive  microphone  and  amplified 
using  a  lock-in  amplifier.  A  microcomputer-driven  data  acquisi- 
tion system  was  used  to  step  the  modulation  frequency  and  to 
record  the  data.  By  systematically  varying  the  modulation  fre- 
quency of  the  e-m  radiation  illuminating  the  sample,  it  is  possi- 
ble to  obtain  a  profile  of  the  sample's  thermal  properties  as  a 
function  of  depth  beneath  its  surface.  The  purpose  of  the  inves- 
tigation was  to  determine  whether  the  photoacous tic  response  of 
wheat  kernels  can  be  related  to  other  kernel  properties,  such  as 
hardness    and    moisture    content. 


In  this  project,  we  studied  a  number  of  different  wheat 
varieties,  including  both  hard  and  soft  red  winter  wheats.  Our 
studies  indicate  that  some  varieties  of  wheat  exhibit  photoacous- 
tic  responses  measurably  different  from  others,  although  substan- 
tial scatter  amongst  individual  kernels  of  a  given  variety  is 
often  present.  Improvements  in  the  standardization  procedure  for 
the  individual  kernels  and/or  modifications  in  the  apparatus 
which  increase  the  magnitude  of  the  kernel's  contribution  to  the 
photoacous tic  response  may  help  to  improve  the  separation  between 
varieties,  and  decrease  the  kernel  to  kernel  scatter  found  in  our 
results . 
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